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CXVIII. The Four Possible Waves in a Magneto-Ionic Medium 


By J. H. Prppineton 
Division of Radiophysics, C.S.1.R.0., Sydney, Australia* 


[Received June 20, 1955] 


SUMMARY 

It is well known that several different types of waves may be transmitted 
through a magneto-ionic medium ; that is, ionized gas permeated by a 
steady magnetic. field. Foremost is the pair of radio waves whose 
properties are so well described by the ‘ magneto-ionic theory ’ of Appleton 
and others. More recently several hydromagnetic (or magneto- 
hydrodynamic) waves have been studied. A ‘plasma’ or electron- 
pressure wave has received little attention except in the simple case of 
zero magnetic field. Finally there are sound waves which may propagate 
in any gas. 

The relationship between these various waves is disclosed : the magneto- 
ionic medium is quadruply refracting, transmitting four, and only four, 
different waves. The hydromagnetic ordinary and extraordinary waves 
are identical with the radio O and E waves, the latter being the simple 
limiting case of the former when heavy ion motion is negligible. The 
other two waves comprise a ‘magnetic sound ’ wave (a type of hydro- 
magnetic wave) and a ‘ magnetic plasma ’ (P) wave. 

The P wave is discussed in some detail, using transport equations. Its 
velocity and polarization are determined within certain ranges and its 
spontaneous growth and decay are considered. Examples of P waves 
in the solar atmosphere are given. If, as seems probable, these waves 
can grow, then they may cause the O and E waves to grow and so provide 
an explanation of non-thermal solar radio noise. 

The effect of electron pressure on the O and E waves is determined 
quantitatively. 


§ 1. INTRODUCTION 


A MAGNETO-I0NIC medium, comprising ionized gas permeated by a 
steady magnetic field, may transmit magneto-ionic waves which are 
electromagnetic waves associated with oscillations of the ions. Radio 
physicists are well acquainted with the ‘magneto-ionic theory’ of Appleton 
(1932) and others which has proved so satisfactory in explaining radio 
wave transmission in the ionosphere and the solar atmosphere. According 
to this theory the magneto-ionic medium is doubly refracting, transmitting 
an ‘ordinary’ (O) and ‘extraordinary’ (E) wave. In the present 


r. . 
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paper it is shown that the magneto-ionic medium is in fact not doubly 
but quadruply refracting. 

The original magneto-ionic theory was developed to deal with waves 
of such high (radio) frequencies that perturbations of the heavy ions of 
the gas may be neglected.* Furthermore, both electron and heavy ion 
thermal motions and the resulting pressures of the electron gas and heavy 
ion gas are neglected.} If we continue to neglect gas pressure but decrease 
the wave frequency until heavy ion motion is important the two waves 
change gradually into hydromagnetic (or magneto-hydrodynamic) waves 
(Astrém 1950, Piddington 1955a). The O and E waves retain their 
identity throughout the frequency change,} but become more complex 
when heavy ion motion is considered. In fact the radio O and E waves 
are really simple limiting cases of hydromagnetic O and E waves when 
heavy ion motion is neglected. 

In addition to the O and E waves, two others may propagate in the 
magneto-ionic medium provided the gas is not at absolute zero tempera- 
ture. These waves depend in part on the elastic forces resulting from 
compression of the gas (either as a whole or the electron gas alone). 
The waves vanish when these forces are neglected. Nevertheless they 
are true electromagnetic waves, having associated electric and magnetic 
fields and a Poynting flux indicating a flow of electromagnetic energy. 
One of the waves appears at hydromagnetic frequencies and probably 
does not propagate in the higher frequency part of the spectrum. The 
other occurs at radio frequencies ; it is not yet clear how far into the 
low-frequency spectrum this wave may extend. 

The low-frequency wave (which might be called a “ magnetic sound ’ 
or S wave) makes its appearance when thermal motions of the heavy 
gas particles, both ions and atoms, are taken into account ; that is when 
the whole gas has a finite pressure (Piddington 1955a). It may be 
considered either as a third hydromagnetic wave or as a sound wave 
disturbed by the steady magnetic field. It differs considerably from the 
other two hydromagnetic waves: the gas movement tends to be longi- 
tudinal rather than transverse and the velocity is of the order of the 
velocity of sound rather than the familiar ‘ hydromagnetic velocity ’. 
When propagation is along the magnetic field the wave is almost a pure 
sound wave, having an associated electric field only. 

The S wave has already been discussed in some detail and the (fourth. 
degree) dispersion equation for a magneto-ionic medium with all four | 
waves present has been given and discussed (Piddington 1955 b). In 


* The effects of heavy ions on wave polarization has been investigated by 
Baker and Green (1933). 

+ Electron motions are considered insofar as they cause collisions with heavy 
ions and atoms and so absorption of the waves by a scattering process. 

+ Although there may be portions of the spectrum where the refractive index 
of a particular wave is an imaginary quantity so that the wave cannot propa gate. 
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the present paper we will, for simplicity, confine ourselves to (radio) 
frequencies where heavy ion motion may be neglected. This eliminates 
the S wave and reduces the O and E waves to the relatively simple form 
of radio waves instead of hydromagnetic waves. 

The fourth (high-frequency) magneto-ionic wave appears when electron 
thermal motions and electron gas pressure are taken into account. At 
the frequencies concerned the medium may be regarded as triply 
refracting. ; 

It is the object of the present paper to 


(a) present the overall picture of the four magneto-ionic waves ; 
(6) indicate the characteristics of the fourth magneto-ionic wave ; 


(c) show the effects of electron pressure on the O and E waves; it is 
not certain that this effect is negligible in regions of high temperature such 
as the solar corona ; 


(d) consider the possibilities of spontaneous growth of the fourth 
magneto-lonic wave and its effectiveness in generating solar radio 
emission. 


§ 2. Puasma OSCILLATIONS AND WAVES 


The reason for the neglect of the third magneto-ionic wave may be 
traced historically. When the magneto-ionic theory was developed, 
electron pressure terms were omitted for simplicity and the great success 
_of the theory seemed to justify this step, at least as far as propagation in 
the ionosphere was concerned. ‘This is so as far as the O and E waves are 
concerned, but what was not realized was that the omission entirely 
removed one root from the dispersion equation and one possible wave 
from the medium. 

In the meantime a new type of electrical oscillation was studied, not 
as a wave travelling through the ionosphere, but as a stationary oscillation 
in a discharge tube. This is usually called a plasma or space-charge 
oscillation (Tonks and Langmuir 1929). Expressed mathematically in 
the form exp i(wt—kx) for variation in time t and space x it has the 
restriction that k is zero or imaginary. Seeking standing waves in a 
magneto-ionic medium Martyn (1947) found two as opposed to Tonks 
and Langmuir’s one. However, neither of these two corresponds to the 
earlier wave; they are the O and E waves in a particular part of the 
spectrum where the real part of & is zero. A change in frequency or 
electron density will cause them to propagate through the gas. To 
these two must be added the third ‘ plasma’ (P) wave which continues 
to exist when a steady magnetic field is added. The presence of three 
stationary oscillatory modes was noted by Westfold (1949) and the 
presence of three travelling waves may be inferred from Bailey’s (1948) 
equations, although these have been complicated by the introduction of 
a steady electron drift so that there seem to be twelve different waves 
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present. The nature of all these ‘waves’ has been discussed briefly 
elsewhere* (Piddington 1955 b). 

When electron pressure is considered, the P wave (which might be 
called a ‘magnetic plasma’ wave) becomes a travelling electromagnetic 
wave and the medium triply refracting. The dispersion equation is now 
of the third degree and considerably more difficult to discuss as it does not 
factorize. 

The P wave has been studied for propagation along the magnetic field, 
when it becomes a pure electric wave and quite independent of the two 
electromagnetic waves 0 and E. Thomson and Thomson (1933) found the 
dispersion equation k2 Bm : 

a ary 1—X) 
WwW K 


where m, 7’, e and N, are the electron mass, temperature, charge and 
density, « is Boltzmann’s constant, X=47N,e?/mw* and f is a constant, 
equal to unity in this case. Electrostatic units are used. Bohm and 
Gross (1949 a) derived a similar equation starting from Boltzmann’s 
equations instead of using transport methods; they found B=}. The 
P wave has also been studied by Gross (1951) for the case of propagation 
perpendicular to the magnetic field. He finds some surprising results 
and incidentally demonstrates that the use of Boltzmann’s equations in 
the extended (3 wave) magneto-ionic theory is likely to prove too difficult 
mathematically. Whatever the inaccuracy entailed, the use of transport 
methods seems justified at least in a preliminary survey. 


§ 3. THE DISPERSION EQUATION 


The equations necessary to describe the three magneto-ionic waves 
are Maxwell’s field equations and Maxwell’s momentum transfer equation 
for an electron gas. For the time being collisions of electrons with heavy 
ions are neglected. The equations are ; 


- ere 
curl H=47c~1j-+-c¢ re .. 4 Glen 
oH 
curl E= —c~! rae dae! he ieoy ON 
div E= dr Nes, 5) 2 ni ee ee 
d 
m . +N~1 grad p,t+e(E+c-1v x H)=0. «ee ops Naren ee 


* The general dispersion equation, derived for a stationary observer in a 
magneto-ionic medium, is of the fourth degree in k? or w?. Each type of wave 
is represented by two waves moving in opposite directions but otherwise 
identical. When the observer moves relative to the gas he sees waves which 
have undergone a Lorentz transformation and the pairs of waves are no longer 
identical. In addition the drift of the gas or observer appears to introduce 
further waves ; these are not real, however: they have velocity equal to the 
drift velocity and mean, physically, that any stationary irregularity in the gas 
appears as a travelling wave when seen by a moving observer. 
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Here H, E and j are respectively the magnetic and electric fields and the 
current density and v, N and p, the electron gas velocity, density and 
pressure. For small oscillations about a steady state, products of two 
variables may be neglected so that eqn. (4) becomes 


0 , 
5p tNo ir grad n+ “E+vxQ=0 arg ise ee eT) 


where n=N—WNp, Hy is the steady part of H, and Q=H,e/mc. Also 
grad p, has been replaced by mz grad n where 7 is one-third the mean 
square random velocity of the electrons, and electron temperature 
fluctuations are neglected. 


fee eV. weet eee. <3) AO ee) 

Consider a right-handed system of axes (x, y, z) and a plane wave of the 

form exp i(wt—kx) propagated along Oxv. Let @ have components 

Q, and 2, along Ox and Oz respectively. The vector equations may be 

written in terms of the perturbation components v,.,,, H,,, and H 
of N and of the wave frequency and wave number. They become : 


tw,+ Qpv,te/mH,—tkN,-1™m=0 | 
sarge sen 8) | Cee en ea) 
1wv,— 27v,+e/mH,=0 
4rN ev, twH,=0 
4 N ev,—twH,+ickH,—0 ae aS) 
4 N ev,—twH,—ickH ,=0 


“LY 2, 


ckE,+wH,=9 (9) 
ckhE ,—wH,=0 
(hE anes 0. Yen 2s) 2. (10) 


Using eqns. (9) and introducing the function X, defined above, eqns. 
(8) reduce to 
wmXv,=1eH , 
2),2 
wmXv,=eL, (- =) 


ww? 


; 24,2 
womXv,=1eH, (1 Z iy 


w 


(11) 


These equations and (10) may be used to eliminate v and n from (7); 
the new variables Y,= 2,/o, Yp= 2 7/o are introduced. 


2 af 
(\-x- =) Bay, e “5 )B,=0 
(62) 
ays 2 
(1—-x- By i¥ (0-S) E,A1Y pH,=0 aL) 
(60) 


c2k2 : Crk 
(1-x- =) B,AtY, (1 ==) E,=09. 
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Eliminating E we have the dispersion equation 
rk2 2h? 2 2k? 2 
(-- 3) (3) -ro(-S} 
2 ri) 7k? 
—Y>p (1-<F)(1—x— >) == 0. (13) 

The basic equations used above are similar to those employed by Bailey 
(1948) with his steady drift velocity omitted. Equation (13) may be 
derived from Bailey’s full (twelfth degree) dispersion equation by reducing 
the drift velocity to zero. The preceding derivation was given to avoid the 
confusion due to the ion drift and because eqns. (7) to (12) are needed 
later to determine the form of the P wave. 

Equation (13) is of the third degree in /? indicating three pairs of electro- 
magnetic waves ; these are the O, E and P waves. When 7=0 one pair 
of waves disappears and, as Bailey has shown, the equation reduces to the 
well-known magneto-ionic equation. 

To introduce the remaining (S) magneto-ionic wave, eqn. (6) is replaced 
by j=N,e(v’—v) where v’ is the positive ion perturbation velocity and 
is given by another equation similar to (4) describing heavy ion momentum 
transfer. The dispersion equation is then of the fourth degree indicating 
four pairs of waves. This equation is difficult to use and if the S wave 
is to be studied some simplifying assumptions are desirable. A convenient 
simplification is effected by introducing the electrical conductivity tensor 
(Piddington 1955 a) in place of the two equations of momentum transfer. 
At the low frequencies concerned the displacement current may also be 


neglected and the properties of the three waves then found without exces- 
sive difficulty. 


§ 4. LONGITUDINAL AND TRANSVERSE PROPAGATION 


In the special cases of propagation along and perpendicular to the 
magnetic field, eqn. (13) factorizes. In each case two of the modes 
remain linked with a quadratic dispersion equation, while the third 
separates out. 

For longitudinal propagation (Yp~—0) the P wave separates out; its 
dispersion equation, with the effect of electron collisions with heavy 
particles included (frequency v sec~?), is 


k? 
as! Se he 
Fa peek (1—X —1Z) 


where Z=v/w. This shows the way in which electron collisions cause 
wave attenuation, just as found earlier by Bohm and Gross (1949 b) 
using a microscopic treatment. The wave velocity is (r/1—X)¥* which 
differs by a fixed factor from that found by Bohm and Gross whose 
dispersion equation is 


ing | 


w 37 


(lesz): 
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The O and E waves are exactly as in the simple magneto-ionic theory 
being unaffected by electron pressure. 

When propagation is perpendicular to the magnetic field (Y,;=0) the 
O wave separates out, once again unaffected by the electron pressure 
term. The E and P waves are linked in a quadratic equation whose 
roots are (writing Y for Y,) : 


ke 1 [eee ee V2 We) inl \ eee 2 
ae (=X) (+ 3) = = - 5{0-") te > 7 a= 
4Y2X ik 


«| ea TF 
The positive sign corresponds to the P wave and the negative to the E 
wave. ‘The equation simplifies for very small values of the second term 
in square brackets. When c?>7 this condition is 
47Y2X 
Ce eeara 4 
and the two waves are given approximately by 


Vhs 
=r \(1-X—Y?) 


Paes, “ X(1—X)(1—7/c?) ) 
ot { ~ (Xx) — 78) — 2 
The first is the P wave having velocity (7/1—X— Y*)!?, provided this 
velocity is much less than c; otherwise the full equation must be used. 
The second is the EK wave and reduces, when 7—0, to the well-known form 
(Booker 1934) of the original magneto-ionic theory. Thus electron 
pressure only changes the E wave appreciably when 7 is comparable 
with c¢ (temperatures of order 108 °K or more) in which case a relativistic 
treatment should be used ; or when the refractive index is nearly zero. 
Doubt has been cast on the above results by Gross who compared the 
results of the transport treatment with his own derived from the Boltz- 
mann equation. Under certain conditions his dispersion equation for 
the P wave differs considerably from that found above and assumes most 
unexpected forms. In particular the medium does not transmit any 
wave when 1/Y is an integer; thus the transmitted spectrum has a series 
of narrow gaps at frequencies which are multiples of the gyro-frequency. 
Also Gross’s dispersion equation takes a more complex form : in the case 
when rk?/w? < Y? it has three different forms, as follows : 
Y(1—X — Y?) 
ee 
AY? Wf? 27-"X + ¥*7—1), 
Y?(X+ Y?—1) 
pf ; 


and 


when Se iio Sr 


X<¥? Plw®=271 


The last two waves only propagate when X-+Y?>1, that is when, 
according to the equation derived above, the medium is overdense. 
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Gross remarks that “ the same type of difference can be found if one does 
not restrict the discussion to the plasma-type solutions, but considers 
the strongly coupled waves when ck/w~Y, X4¥?”’ (that is, the E wave). 
' Whatever inaccuracies may eventually be found to result from the use 
of transport methods, their use seems justified in a preliminary survey 
of obliquely propagated waves. The use of Boltzmann’s equations would 
appear to offer excessive mathematical difficulties. 


$5. OBLIQUE PROPAGATION 


In the general case, where propagation is neither along nor perpendicular 
to the magnetic field, eqn. (13) does not factorize and a full discussion 
becomes lengthy, beyond the scope of this paper. However, with the 
restriction that 7+<c? the waves may be separated by using our prior 
(approximate) knowledge of their characteristics. 

Let us assume that the values of the refractive index (ck/w) of the O 
and E waves are not substantially changed by electron pressure. In the 
regions of main interest ck/w lies between zero and unity and so 

tk? rok 

we ae <P 
Replace rk?/w? in eqn. (13) by a small constant (6); the two roots corre- 
sponding to the O and E waves are then given by 


rend, ee Y,2 ( Y,2 2 ol v2 ]-2 
sot Dg resp as YP : sual 


\ 

This reduces to the well-known magneto-ionic equation (Booker 1934) 
when 6=0 and, in general, is affected only slightly by the presence of 6. 
However, in regions of nearly zero or very high refractive index the value 
of ck/w may be more critically dependent on 5. Also when 7~c? it is 
evident that the effects of electron pressure may be of major importance. 
It may be noted that 5 only affects the Y, terms and so increases in 
importance as Y/Y, increases. 

To separate the P wave we assume its phase velocity is much less than 
c, that is c?k?/w?>1, and also that X is of order unity or less. Equation 
(13) then reduces approximately to 


which may be written 

Na ? YY, 

wot (I-X— 55). er ee aS) 
It may be noted that this equation does not hold in the interesting range 
where the term in brackets approaches zero, that is, near the natural 
plasma frequency. Here the phase velocity approaches and exceeds 
¢ and our simplifying assumption no longer holds. At yet higher phase 
velocities when c?k?/w? <1 the dispersion equation takes the form 


k2 ie = Y,°(1—X) 
co) aman Hinx- eh, é : (16) 
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: Some interesting properties of the P wave may be displayed by plotting 
p*(=ck?/w?, the square of the complex refractive index) against X. This 
method was used by Ratcliffe (1933) to show the characteristics of the 
O and E waves. The curves are drawn for longitudinal propagation and 
transverse propagation ; the region between the curves corresponds to 
the infinite number of oblique waves. One of Ratcliffe’s curves is drawn 


Biel 


nN 


i 
(@) 
a) 


Ws 


or 


\ 
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O&E Waves \\\\\\\\ 
PWave “2 


Areas on the 2X diagram corresponding to the O and P waves. Full lines 
correspond to transmission along the magnetic field (Yp~=0), dashed 
lines to propagation across the magnetic field (Y;=0). The hatched 
areas between these lines correspond to oblique propagation. 


in fig. 1 for the O and E waves ; the full lines correspond to longitudinal 
propagation and the dashed lines to transverse propagation. The 
corresponding curves for the P wave are also shown and the two areas of 
oblique propagation are cross hatched. The equations used were those 
of §4, using the full transverse equation. The value of 
Y(=(Y,?+Y,°)12) chosen was 0-5 and real values of » are seen to exist 
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for values of X up to 1+ Y. The value of c?/7 used was 4 which is very 
small, but was chosen to make the diagram a reasonable shape. It 
corresponds to an electron temperature of ~1-4 x 10° °K. 

The polarization of the P wave may now be determined, unfortunately 
only in the relatively uninteresting region away from the resonance 
frequency. Assuming c?k?/w?> 1 eqn. (12) yields, approximately, 


E, =—iY, 
EB, --ictk® /1—Y,? 
oe i, wat, \GReeys 


Thus, provided Y, is of order unity, #, and H, are of the same order but 
out of phase (elliptically polarized) and provided 1— Y,?/Y > is of order 
unity, H,>H, (or H,). The greater part of the electric field associated 
with the wave is a potential field (see eqns. (3) and (10)) and lies in the 
direction of propagation. 

The components of H may be found from eqn. (9) 


“hk 
a= ~iY,— He 


k 
and i= ie 
(62) 
Jal . 
so that Tie Fy 


the waves showing elliptical polarization of the magnetic field with 
Hoa. 
The components of v are found from eqn. (11), remembering that 


; Pes gx : 

Ps 5 L Hae 
= = —4 ( 7 and —=—iY,. 
() 1 VU 

1 


Thus the electron velocity components tend in general to be of the same 
order of magnitude. 


§6. PLasmMA Waves IN THE SOLAR ATMOSPHERE 


Perhaps the region in which plasma waves may have the most interesting 
observable effects is the solar atmosphere. They may, for example, be 
a link in the mechanism of the non-thermal radio emission which requires 
ordered motion of the electrons (Piddington 1953). Magnetic fields exist 
throughout a large part of the solar atmosphere but particularly above 
sunspots where their distribution of strength and direction may be 
inferred with reasonable accuracy from visual observations. The 
characteristics of some possible plasma waves are given for such a region, 
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The phase velocity of the P wave is given by 
y— = =71/2 F(X.Y) 


where F is given by eqns. (15) and (16) for the restricted conditions 
concerned but has not been determined for the interesting velocity range 
neare. Itis likely that F has some roughly similar form within this range 


Mognetic fickd strength(gauss)for Frag. 100 Mcf. 
356 ‘356 356 3 


10 


Wave velocity = V em, sec, 7! 
Wavelength (cm) for Frequency 100 Me/s. 


} - 
ee ne) 10 
Gyro to wave frequency ratio-Y 
Plasma (P) wave velocities plotted against the ratio Y of electron gyro-frequency 
cae wave Rncecas for several different values of plasma frequency X. 
Also shown are the corresponding wavelengths and magnetic field 
strengths when the wave frequency is 100 Mc/s. 
and so, for illustrative purposes, its value is taken from eqn. (15) and 
assumed to apply when V=c. It is assumed that the actual values of 
V are in error when V ~c but that the general form of the velocity function 
so found is representative of the true function. 
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In fig. 2 the values of V are plotted for P waves in a gas at temperature 
107 °x which corresponds to the ‘hot regions’ which occur normally 
above sunspots (Piddington and Minnett 1951). The waves move obliquely 
across the (uniform) magnetic field at an angle of 7/4 and V is plotted for 
Y(=\/2Y,=v/2Y 7) ranging between 0-1 and 10. The corresponding 
values of magnetic field-strength when the wave frequency is 100 Mc/s 
are plotted, ranging from 3-56 gauss to 356 gauss. Five values of X 
are chosen ranging from 0-6 to 3-0; the electron density corresponding 
to X=1 for 100 Mc/s waves is approximately 108 cm~*. The wavelength 
of 100 Mc/s waves is also plotted, ranging from 1-0 to 10? cm. 

For values of X<1-0 the V curves are discontinuous, being split into 
the ranges Y>1 and Y<1. The former has values of V<7!/?, falling 
to zero when asymptotic to the line Y= 1/2. The latter has values of 
V>7!/, rising to infinity at a resonance value of Y. When X>l, 
each V curve is continuous and complete in the region Y>1. They are 
all asymptotic to the line Y= 1/2, at which value V0; they are also 
asymptotic to larger values of Y at which the plasma is resonant and 
V+oo. The curves for any other value of 7 may be found by simply 
changing the vertical scale, remembering V ocr”. 

The minimum wavelength of a plasma wave in the absence of a magnetic 
field is substantially greater than the Debye wavelength (Bohm and Gross 
1949 a) which is («7'/47ne*)1/*. For the frequency 100 Me/s and n~108 em-* 
(ile. X~1) this is about 2-2 cm so that the waves having V=71!/? cannot 
propagate. However, the presence of the magnetic field may assist 
propagation at lower values of V because of the restriction on electron 
motion. The electron gyro-radius when H=40 gauss is about 0-7 cm. 

The curves of fig. 2 suggest interesting possibilities of wave growth and 
ion acceleration. However, in discussing these possibilities it should be 
borne in mind that in the region V Sc the curves are at best inexact. 


§7. Growing P Waves anp LINKED WAVES 


Above sunspots there will be marked gradients in the strength and 
direction of the magnetic field and hence in V. The effects of gradients 
in X are likely to be smaller except when hydrodynamic shock waves are 
present, or P waves which are particularly sensitive to changes in X. 

Gradients in V may lead to growth or decay of the P wave as shown by 
Bohm and Gross (1949b) when a magnetic field is absent. Growth 
results when the wave velocity is decreasing and when a proportion of the 
ions present have velocities equal to the wave velocity. The same 
mechanism will operate in the presence of a steady magnetic field, but 
as the ions may have steady velocities only along the field, their velocity 

must now be V sec % where tan J=Y,/Y,. The ions are decelerated 
by the component of the wave electrostatic field parallel to the steady 
magnetic field, that is to H, sin j (for the wave velocity range discussed 
in §5). This type of wave growth is likely to occur wherever V<c 
and is decreasing away from the resonance region, owing to decrease in Y, 
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The wave will continue to grow at the expense of the fast ions which are 
retarded with the wave. This continues until V—71/2 and the wave 
either disintegrates or decays, its energy resulting in local heating. 

This mechanism, operating in reverse (with V increasing), may cause 
ion acceleration and could perhaps be responsible for solar cosmic rays. 
Provided V increases slowly enough the ions will be accelerated until 
Vc when they must leave the wave and start moving independently. 

A numerical example of a growing P wave may be useful. Consider 
the wave corresponding to the X=0-8, Y<1 curve of fig. 2. In moving 
away from its resonance boundary (Y=0-577) it decelerates from V—= oo 
through ¢ down to V=2-8x 10° when Y=0-1. Suppose there are a few 
electrons having a component of velocity along 0x lying between 8-1 x 109 
and 1:75 10° cm sec”. The wave decelerates through this range as 
Y falls from 0-57 to 0-55; over a medium sunspot such a change occurs 
in a distance of about 400 km for waves of frequency ~100 Mc/s (Pidding- 
ton and Minnett 1951). The time of travel of a wave crest through this 
region is about 3 10-3 sec so that it would suffer negligible absorption 
due to electron collisions with heavy particles. If there were, say, 10 
electrons per cm? having an average velocity of 101° cm sec~! and they 
lost half of their kinetic energy to the wave, the increase in energy would 
be about 210-7 erg cm-?. If this energy is turned into energy of the 
electrostatic field of a P wave it corresponds to a peak value of about 
1 volt cm. 

It is shown above that P waves may never emerge from the solar 
atmosphere and so cannot provide a direct explanation of the intense 
solar radio emission which reaches the earth. However, an indirect 
explanation may result in terms of coupling between the P wave and the 
O or E waves. The simple theory of weak waves shows no such coupling, 
but this theory only applies rigorously for infinitely weak disturbances. 
For finite disturbances and particularly for waves which have grown to 
considerable strength as indicated above, there is coupling between the 
waves. This coupling becomes effective in transferring and equalizing 
wave energies when the two waves have equal or nearly equal velocities 
and so remain in phase over many oscillations. We have seen that the 
P wave, within certain regions, has the same velocity as the O and E 
wave, so that strong coupling should result. The transfer of energy from 
the strong (P) to the weak (O or E) wave is inevitable since they comprise 
two systems resonating at the same frequency and linked by a common 
component of H or E. The four magneto-ionic waves may be considered 
as mechanical oscillations of ions with which are associated electromagnetic 
fields, the latter providing the mechanical forces on the ions. When 
two such oscillations resonate, energy must be transferred from the 
stronger to the weaker because they will share, in general, a common 
part of the electromagnetic field. 

A partial theory of intense solar radio noise is now clear. The P wave 
grows by one of the well-known mechanisms, and having become strong 
enough to introduce non-linear effects it transfers part of its energy to an 
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O or E wave having the same frequency and wavelength. The P wave 
then disintegrates, leaving the O or E wave to emerge from the solar 
atmosphere. The electric field calculated above and extending over a few 
hundred km would, when transferred to the O or E waves, be sufficient 
to explain solar radio bursts of moderate intensity. 


§ 8. OrHER PossipLE Maaneto-lonic WAVES 


The four magneto-ionic waves discussed in § 1 may occur in a highly 
ionized gas in which electrons and heavy ions predominate. In a lightly 
ionized gas, in which neutral atoms or molecules predominate, yet further 
types of waves may propagate. First consider the ions in this gas and 
assume that they suffer no collisions with neutral particles; they may 
then oscillate in the four modes already discussed. If a few collisions 
occur between ions and neutral particles these do not essentially alter the 
modes but cause dissipation of the waves by scattering. When collisions 
are sufficiently frequent, enough energy is transferred to cause the heavier 
and more sluggish mass of neutral particles to oscillate. Each neutral 
particle must suffer many collisions in a wave period which means that the 
frequency range of these oscillations is some orders of magnitude below 
that of the hydromagnetic waves in the ion plasma. 

This process may be extended further by including other heavier 
particles (for example smoke) in a sufficiently tenuous gas which, in turn, 
is lightly ionized. It hardly seems desirable, however, to describe such 
waves as magneto-ionic waves of which either three (at radio frequencies) 
or four (over a wider frequency range) may be taken as fundamental. 
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ABSTRACT 


Thin gold films deposited on a substrate of bismuth oxide or, better, 
between two films of bismuth oxids, are found to have exceptionally 
high electrical conductivity and optical transmission. Heating to 
several hundred degrees centigrade further enhances these properties, 
irreversibly. Typical results, for gold thickness 100 A and heat-treatment 
at 450°C, are film-resistance 3} ohms ‘per square’ and white-light 
transmission 75° (or 82°% using controlled interference effects). 

The measured electrical and optical constants are remarkably close to 
those of gold in bulk. This appears to be due to a very fine-grain, con- 
tinuous structure. Electron microscopy and diffraction show that even 
at 60 A thickness no holes of resolvable size exist in the film, that the 
lateral dimension of the gold particles is some ten times the film thickness, 
and that the enclosing bismuth oxide films are amorphous, at least 
initially. The effect of the heating appears to include a considerable 
reduction in internal surface and also a marked reduction in surface- 
roughness. Surface reflection of conduction-electrons is found to change 
from fully diffuse to almost fully specular, as a result of the heating. 

The function of the bismuth oxide seems to be to preserve, up to quite 
high temperatures, surface energy conditions controlling mobility of 
grains, and opposing increase in gold surface area. This would account 
for the absence of aggregation during deposition, and of grain-growth 
perpendicular to the surface during heating which would otherwise 
disrupt the films at much lower temperatures. In this respect. and 
possibly also in promoting a smooth gold surface, the oxide enclosure 
seems to simulate one of gold itself. 

The films have a possible use as transparent conducting coatings, for 
the purpose of de-icing aircraft windows by direct electrical heating. 


§1. INTRODUCTION 


Durtxe work on transparent, electrically conducting film-coatings for 
glass, two of the authors (E. J. G. and J.S. P.) were asked to develop a 
coating suitable for de-icing aircraft windows by direct electrical heating. 
A power dissipation of about 1 w/cm? was envisaged. With a low voltage 
of about 100 volts, applied between bus-bars along the longer edges of 
the usual rectangular window, such a high dissipation is possible only 
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with a low film resistance of about 10 ohms ‘ per square ’ (i.e. between 
line-contacts along opposite edges of a square specimen). At the same 
time, the light-transmission had to exceed 70%. As a coating-material, 
cadmium oxide (Preston 1950) in a thickness giving 10 ohms per square 
was found too yellow and absorbent. Stannic oxide (U.K. Patent 
Specifications Nos. 632, 227 and 632, 256), though very transparent, was 
found difficult to prepare with adequate surface uniformity. Other semi- 
conductors appeared still less favourable. 

So, attention was turned to metal films—particularly gold films because 
of the high specific electrical conductivity and chemical stability of this 
metal. Simple calculation, from the accepted bulk-metal constants, 
gives a white-light transmission of over 70°, and resistance well below 
10 ohms per square, for an ideal lamina of gold of thickness 100 4. These 
are much more favourable values than any experimental ones revealed 
in a search through the literature on. gold films. However, the simple 
calculation neglects two factors. First, on the assumption that the 
surfaces of a conductor scatter or reflect conduction-electrons diffusely, 
the surfaces give rise to an excess component of resistance which, while 
generally entirely negligible, becomes practically significant in films whose 
boundary surfaces are closer together than, say, 100 or 200 A (depending, 
of course, on the electronic mean-free-path for the material and tempera- 
ture in question). At the start of the work, such an effect was taken to be 
inherent and unalterable. However, results were obtained later which 
indicate that electron scattering at the film-surfaces was in fact almost 
eliminated by the chosen modes of preparation and treatment of the 
films. The second point neglected in the simple calculation is the more or 
less granular structure of real gold films, contrasting with the ideal laminar 
model. This is a matter of common observation (Sennett and Scott 
1950), and is held generally to be the chief reason why observed values 
of optical absorption and electrical resistance for thin gold films are 
excessive. 

The present work started, then, with the object of finding a way to 
control granulation, and to hold the transmission and resistance of the 
resulting films consistently within satisfactory practical limits, close to 
those predicted by the simple calculation. The guiding principle was that 
the degree of granulation is controlled by the balance between the adhesive 
forces between the gold and its substrate (favouring ‘ wetting ’, and con- 
tinuity of the film), and surface tension in the film (favouring disruption 
and granulation). Search was then made for a substrate-film material 
which could be relied on to exercize the desired influence on the structure 
of the superposed gold, would be itself highly transparent and, of course, 
would adhere well to glass. The search was successful, and preliminary 
observations have been published (Gillham and Preston 1952, article in 
The Chemical Age, July 5, 1952). Bismuth oxide was selected for the 
substrate material. This paper describes more fully the preparation of 
gold films upon it, and their distinctive electrical and optical properties. 
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The electrical conductivity of the bismuth oxide is entirely negligible 
compared with that of the gold, and so is left out of account. 


§ 2. PREPARATION AND GENERAL DESCRIPTION OF FILMS 


The films were all sputtered, in a plant and by a method already 
described (Preston 1950). The essential conditions were :—current- 
density at cathode 1-5 ma/em?, cathode dark-space 5 mm, specimen to 
be coated situated just outside the dark-space. For the gold films, a 
continuous flow of argon (99-8°%, commercial) at suitable pressure was used. 
For the bismuth oxide films, the cathode was of bismuth, and air was 
added to the inflowing argon in proportion just sufficient to ensure deposi- 
tion of clear oxide films obviously free from excess metal. The exact 
composition of such films is not certain, but excess oxygen in their 
preparation is to be avoided. Both cathodes, of fine gold and of bismuth, 
were first etched by dummy sputtering-runs until clean of all foreign 
matter. 

The thickness of bismuth oxide films so made was estimated by optical 
methods, high accuracy not being necessary. The rate of deposition was 
found to be 6-7 A per second. For the gold films an accurate estimate of 
thickness was desired, also preferably in terms of sputtering-time since 
direct measurement of each sample would be neither easy nor convenient. 
By relative colorimetric estimation of colloidal gold, taken into solution 
from a series of films, it was shown that the mass of gold deposited was 
proportional to sputtering-time, between 0 and 20 seconds. Then, the 
absolute rate was found by weighing several 20-second specimens, and 
calculating thickness on the basis of the bulk-density of 19-3. The rate 
was found to be 1114 per second. All gold-thicknesses quoted here- 
after are on this basis, which may not of course give the exact physical 
or optical thickness of films, in so far as their density might not be 19-3. 
It is, however, the only systematic basis for discussion. As an extra 
precaution, in making a series of gold films of different thicknesses, an 
additive method was used. For instance, specimens A and B were 
sputtered simultaneously for X seconds, then A and C simultaneously 
for X seconds, so that the thickness of A was accurately the sum of the 
thicknesses of B and C. This step-wise doubling of thickness afforded a 
check on, and correction for, slight variations in the rate of deposition 
over the whole series. 

Tt was found that a 1004 film of bismuth oxide was sufficient to obliterate 
all influence of the necessary glass supporting base. On this film, a 60 A 
gold film had initial resistance 20 ohms per square, overall light transmission 
(tungsten 2700°x) 73%, and colour pale straw. Heat-treatment for a few 
minutes at 200°c was found to improve these properties, irreversibly, to 10 
ohms per square and 75%. At higher temperatures the gold apparently 
disrupted, the film changing colour through the more familiar pale green to 
pale blue, with considerable decrease of conductivity and transmission. It 
seemed that up to 200°cincrease in mobility ofthe gold permitted an increase 
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in continuity or compactness of structure, under forces of adhesion, but that 
above 200° surface tension preponderated. To diminish still further the 
surface tension effect, therefore, specimens were prepared in which a second ~ 
100 A bismuth oxide film was deposited on top of the gold. These were 
found to withstand heat-treatment or annealing* up to 350°C or more 
without disruption. In this way, with gold again of thickness 60 A, 
a resistance of 7 ohms per square and transmission 76°, were obtained. 

The best overall technical performance was obtained with 100 A gold 
sandwiched between two bismuth oxide films each of thickness about 
450 A, annealed for a few minutes at 350°c. Such a multiple coating is 
effectively a high-index half-wave layer, so that light-reflection is at a 
minimum. Also, because the relatively thin gold is situated approximately 
in the plane where the resultant electric vector of the interfermg wave- 
trains is at a minimum, its optical absorption is at a minimum. In this 
way a combination of resistance 4 ohms per square and transmission 
82°, was achieved. The colour was then rather dependent on control 
of the thickness of the bismuth oxide, but for the correct thickness was 
like the slightly greenish colour of ordinary thick plate glass. Reproduci- 
bility of both types of sandwich film was found to be good, however, and 
it was also found that higher annealing temperatures than those quoted 
could be used safely, particularly for greater gold thicknesses. For 
instance, a 100 A sandwich film can safely be annealed at 450°c, and then 
has a resistance of only 3} ohms per square. 

To sum up, it was found that a substrate of bismuth oxide favoured 
consistent and low values of electrical resistance and optical absorption 
in the thin gold films. Also, there was the purely empirical observation 
that a final heat treatment was then possible and led to further improve- 
ment—which was very considerable in the case of gold films enclosed 
between two films of bismuth oxide, when the treatment could be carried 
up to a temperature of 350—450°c without break-up of the gold. The 
fact that in this way the resistance of a 100 A gold film could be brought 
down to 3} ohms per square—not much more than the value of 2-45 
deduced directly from the normal bulk-resistivity and the thickness— 
prompted a closer investigation into the resistive properties of such films, 
and measurements of their optical constants, in comparison with those of 
gold in bulk. 


§3. Errects oF ANNEALING, AND VARIATION OF ELECTRICAL 
RESISTIVITY WITH THICKNESS AND AMBIENT TEMPERATURE 


For the electrical study a range of sandwich films was made, all with 
equal upper and lower bismuth oxide films 120 A thick, but with gold 
thicknesses between 55 and 880 A. They were made in triplicate according 
to the step-wise procedure already mentioned. From the measured 


* Hereafter the word annealing is used as the equivalent of heat-treatment in 
the broad sense, for brevity. 
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resistances and linear dimensions of the specimens, the specific resistivities 
were calculated in ohm-cm. Then the mean value for the three specimens 
at each thickness was found. Before annealing, resistance measurements 
were made at ambient temperatures of 20°c and —183°c. After annealing 
at 350°C, measurements were made at 200°c, 20°c and —183°c. The 
mean values of specific resistivity so obtained are given in table 1. 


Table 1. Specific Resistivities of Gold Films in ohm cm x 10-® 


=a 
Thickness of gold in Unannealed Annealed films 
A films at at 
20° + —183°c 200° 20°c- —188°e | 
55 16-8, 15-6 6-6, 48 3-1; 
110 11-9, 10-5 5:6; 40, 2:3 
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Variation of specific resistivity of gold films with thickness, annealing, and 
ambient temperature. 


These results are plotted as crosses in fig. 1. At the right-hand edge 
of fig. 1 three short lines are also given, showing for comparison the 
accepted values for the specific resistivity of bulk gold at the three 
temperatures in question. Each column of resistivities can be fitted, 
within probable limits of error, to an equation of the form 


Film resistivity ep=pp+hy/t eee a (E) 
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where p, and & are constants, and ¢ is the film thickness. Appropriate 
values of p,; and k are given in table 2. 

The five full curves in fig. 1 are those for py, obtained by inserting the 
five pairs of values from table 2 into eqn. (1). They are seen to represent 
the observations quite closely. 


Table 2 
State of Film ppin ohm cm x 10-6 k in ohm cm?x10-¥ 
Diff. 

Unannealed, at 20°C 6:7 1-6 5:7 

A », —183°C 5-1 5:7 
Annealed, at  200°o SL eihee key 0-66 

2”? ” 20°C 3°6 0-66 

i eros each 1:85} 1-75 0-66 


Consider now the interpretation of the results, on the basis of this 
systematized representation. First, pg is the value to which the film- 
resistivity tends for large values of thickness, i.e. it is the specific resistivity 
of the material of the films, in bulk. We see that this is markedly and 
irreversibly diminished by the annealing process. The reduction is 
like the one observed when bulk metal specimens are annealed, but is 
much larger. It may well be due to expulsion of occluded gas, as well as 
to elimination of dislocations, or even to some major change in structure 
from a more or less amorphous to a crystalline form (Kramer 1937). 
It can hardly be due to elimination of impurity in the ordinary chemical 
sense. It is interesting to note, however, that the values of p; in table 2 
conform well to Matthieson’s rule, which often applies to metal specimens 
whose resistivity is excessive because of impurity or of dislocations in 
considerable number. According to this rule one would write 


PB=PatPeo 


where pg is the temperature-dependent specific resistivity of normal pure 
gold, while pg is a temperature-invariant, ‘ residual’, component due 
to the impurity or structural imperfection of the actual specimen. In 
such a case, the change in p, produced by a change in ambient tempera- 
ture would be found equal to the corresponding change in p, Now, 
the reference tables give 1-69 x 10-6 and 1-46 x 10-6 for the changes in pg 
for the steps —183 to 20°c, and 20 to 200°c, respectively. The 
experimental changes in pz (see table 2) are 1:75 10-6 and 1:6 10-8 
for annealed and unannealed films, for the first step, and 1-5 x 10-° for 
unannealed films for the second step. These show quite reasonable 


agreement with the accepted characteristics of pz, in confirmation of the 
rule. 
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Thus, the film material behaves like pure bulk gold in a state of struc- 
tural disorder, in respect of both the irreversible reduction of its bulk 
resistivity by annealing, and the size of the reversible changes in resistivity 
which accompany changes in ambient temperature, before as well as 
after annealing. This supports the view that low resistance has been 
reached in the present films by a process which secures an unusual degree 
of order or continuity of structure in so far, at any rate, as the 
substance of the gold in thin films is concerned. 

Consider, now, surface effects. We may suppose the term k/t of eqn. (1) 
to represent the additional resistivity-component due to diffuse electron- 
reflection, or electron-scattering, at the film surfaces. Qualitatively, 
the observed values conform to this supposition, for from table 2 we see 
that k, or k/t for a given thickness, is not dependent on ambient tempera- 
ture, so that 4/t behaves as another residual-resistivity component, as one 
would expect for such an electron-scattering mechanism. Also, this 
residual component increases as thickness diminishes. This again would 
be expected, for as the thickness diminishes, the surface extent remains 
unchanged while the volume of film diminishes with thickness. Surface 
scattering thus becomes relatively more and more predominant, with a 
correspondingly increasing specific resistivity. It is not, of course, to be 
assumed that the surface component is represented precisely by a term 
varying as the reciprocal of thickness. This simple function is chosen as 
giving a satisfactory fit with measured values over the thickness-range 
covered in the experiments. The important point is that it is seen to be 
unaffected by changes in ambient temperature. 

The really surprising feature of the figures given in table 2, however, is 
the very marked reduction in & produced by the annealing. The con- 
clusion is inevitable that annealing produces some change in the surface 
topography of the films which results in a marked drop in diffuse electron 
reflection at the surfaces. It is worth examining this conclusion in terms 
of electron mean-free-path theory which Sondheimer (1952) has lately 
reviewed. It is sufficiently exact to use his approximate relationship 
for films having thicknesses greater than the mean-free-path. This can 
be written 


oud ; 
pr=po| 1+ 5°; 1-2) | 6 ‘ . 5 G a (2) 


where / is mean-free-path, and p is the specular component of electron 
reflection at the film surfaces. If we compare eqns. (1) and (2) we 
find 1 (1—p)=(8/3)(k/p,). Using values of pz (for 20°C) and & from 
table 2 this gives /,(1—p,)=2-27x<10-*§ cm for unannealed films, and 
1,(1—p)=0-49 x 10 ® cm for annealed films, at 20°c. ne 
For the present purpose, we may assume that mean-free-path is in- 
versely proportional to the basic specific resistivity pz, so that 1, /t,=3-6/ 6:7 
and we find (1—p,)/(l—p,)=0-12. Now (1—p,) cannot exceed unity, 
therefore (1—p,) cannot exceed 0-12, and p, must therefore be at 
least 0-88. The values p,—0, and p,=0-88 are further confirmed because 
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they lead to /,=227 4 for the mean-free-path in unannealed films, and 
1,=408 A in annealed films, both quite reasonable values of mean-free- 
path in metal specimens having the resistivities just quoted. 

There is room for some degree of inexactitude in the above treatment 
and assumptions, but there is clearly strong ground for concluding that 
electron reflection at the surfaces cf the unannealed films is completely 
diffuse or almost so, while for annealed films the diffuse component is 
only about 0-1 and the specular component is about 0-9. This effect of 
annealing is of little practical importance in thick films. In films less 
than about 200 A thick, however, it is more important than the drop in 
py, in lowering the actual resistance of the films. It is believed that this 
is the first investigation in which it has been possible to demonstrate the 
practical possibility of almost completely specular electron reflection at 
the boundary surfaces of a metal specimen. 

The precise structural changes produced by annealing remain for 
speculation. Annealing might well expel occluded gas, and so diminish 
the extent of internal surface in the film and lead to closer packing. It 
would reduce the number of lattice imperfections and dislocations, per- 
haps transforming an almost amorphous layer into a microcrystalline 
one. As already mentioned these effects are quite consistent with the 
observed drop in pz. They neither require nor suggest, however, an 
increase in surface smoothness of the films, while the sort of grain-growth 
which usually accompanies annealing would suggest, rather, a decrease 
in smoothness. In contrast, however, the change from zero to 0-9 for 
the specular electron reflection coefficient at the surfaces seems to require 
development of a very high degree of surface smoothness, with excrescences 
which are either only quite rare or else are of dimensions generally not 
exceeding the effective wavelength ofthe conduction electrons, viz. about 5 A 
(see Miiser 1954, on the diffraction of electrons by certain types of surface- 
structure). A possible structural change would be a preferred growth or 
extension of single-crystal domains along the run of the film but not 
perpendicular to it, perhaps also with a preferred orientation of certain 
crystal planes also parallel to the run of the film. Such a preferred growth 
or orientation need not be assumed to extend far from the surfaces into 
the body of a thick film. In a very thin film however, one would expect 
it to penetrate the whole thickness, and perhaps then be demonstrable by 
electron microscopy or diffraction. 

Thus far, in using eqns. (1) and (2), we have followed the orthodox 
analysis. py, is treated as the bulk-resistivity of the film-material, 
constant over a given set of films all in the same given condition (in 
practice, it may vary somewhat—but there is no experimental check on 
this). The variation of film-resistivity with thickness, over the set, is 
assigned quite separately to properties of the film surfaces. The assump- 
tion is that these properties are in no way related to, or determined by, 
the resistivity or structure of the body of the specimen. For a thick 
specimen, this seems reasonable. In really thin films however, or films 
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built up like the present ones by continued deposition of similar material, 
one might expect to find some systematic relationship between surface 
properties and internal properties, dependent on some single parameter of 
the structure. For example, in wholly granular films, both surface 
roughness and extent of internal surface (or holes) would depend on grain 
size. Again, in the present films, annealing produces two simultaneous 
effects, namely a reduction in roughness and a decrease in p,, the 
temperature-independent component of the bulk resistivity. These 
might both arise from a change in one single feature of the structure. 
We shall now look briefly at this kind of proposition, for what light it may 
throw on the physical nature of the films and of the annealing process. 

Suppose a film to resemble, initially, a slice cut from sponge-rubber, 
full of holes and with surfaces which are rough because covered with 
sections through holes. Given a suitable hole-size, the surfaces will 
cause diffuse scattering of conduction-electrons, giving rise to a 
temperature-independent ‘surface-resistivity’ like the term h/t of 
eqn. (1). Inside the film the holes, being of size similar to that of the 
surface asperities, must also cause similar electron scattering, and give 
rise to another temperature-independent component of resistivity, 
which we may identify with p, since it will not vary with ¢t. If, then, 
annealing be supposed to reduce the number of holes in the bulk as well 
as in the surfaces, both these resistivity terms will decrease together, 
as in fact they are observed to do. This picture is not altered if we 
replace the concept of holes by one of grains, except that we must then 
suppose annealing to provoke a general coalescence of grains, and not a 
growth in the size of the grains at the expense of their number. We may, 
then, suppose our films to be initially granular, but with a smaller grain 
size than usual because of low mobility of gold over the bismuth oxide 
substrate. During annealing, we may suppose the oxide layers to preserve 
large surface tension forces tending to flatten the surfaces of the gold, and 
permitting the escape of ‘ bubbles ’ of gas, occluded in the inter-granular 
spaces, without permanent break-up of the gold surfaces. In this picture, 
the effect of annealing appears as a reduction in internal surface rather than 
in the number of dislocations, along with a reduction in the real extent of 
the external surface. It is probably as good an overall picture of the 
physical nature of the films, and of the effects of annealing, as can be 
drawn from the results of the electrical measurements. 

In this picture, one envisages a grain size somewhere in the range 
between the electronic wavelength and the mean-free-path, in the films in 
their initial state, at least. It is not supposed that, in the unannealed films, 
collisions of the electrons with the film surfaces are in themselves inelastic, 
but rather that the redistribution of electron velocities caused by specular 
reflection at the individual grain facets results, effectively, in a reduction 
in the resultant drift velocity in the neighbourhood, and therefore in an 
excessive measured resistance. This seems to be the usual assumption in 
theoretical treatments of the subject, though it is often implicit rather than 
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explicit. The grain size of about 40 4 revealed by the electron microscope 
examination described later does not conflict with this general picture.* 


§4. Mrruop oF MBASURING THE OPTICAL CONSTANTS 


We turn next to an optical examination. The difference in colour and 
transmission, between unannealed and annealed films, was less striking 
than the differences in electrical properties, and an examination of both 
types of film did not seem likely to give any useful additional information 
on the structure. The optical measurements were therefore confined to 
annealed films, as contrasting with the greener films usually obtained on a 
glass base, and likely to give the closest comparison with the accepted 
constants for bulk gold. 

Sandwich films were used, similar to those for the electrical measure- 
ments. There was a choice of methods, depending on what three 
quantities were selected from which to calculate the three unknowns, 
viz. the optically effective thickness, and the real and imaginary parts of 
the complex refractive index. Transmission, and back and front reflection 
factors are three such quantities, ideally. For very thin films, however, 
the relation between them is such that, unless measured with extreme 
accuracy, they amount in practice to only two independent quantities. 
Thickness, then, was taken as the third. Only that on the mass-per-unit- 
area basis was available, not necessarily exactly equal to the effective 
thickness. In addition, therefore, measurements of phase change on 
transmission and reflection were made. Phase change can be shown to 
be particularly sensitive to thickness (Ishiguro and Kuwahara 1951), 
and so these observations were used to check this, rather than to replace 
transmission and reflection factors as primary measured quantities. They 
were also of considerable use as a guide in certain methods of successive 
approximation used in the computations. 

In the present case, the measurements of transmission and reflection 
factors were necessarily made on the whole composite film, situated 
on a glass base which could contribute an attenuation due to its own 
absorption and the reflection !oss at its free surface. (Troublesome 
inter-reflections were avoided by using a wedge-shaped glass.) Corrections 
were thus necessary for the effects of the base. Further, to compute 
constants for the gold alone from the corrected measurements, certain 
optical data on the actual bismuth oxide present in the composite films 
were needed. The procedure was thus to examine the glass base plates 
first when uncoated. Then the whole area of the plates was coated with 
the bismuth oxide films, but the gold was confined, by masking, to a 
central strip. Thus the composite film was flanked on either side by 
areas having no gold, but an identical total thickness of oxide upon them. 
Then, by making measurements of transmission and reflection, on the two 


side areas as well as on the central strip, all the additional data required 
were obtained. 


* But see also §7, ‘‘ Note added in proof”. 
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Phase changes were measured with the aid of a Twyman—Green 
interferometer, combining the light reflected or transmitted by the films 
with that reflected from a suitable optical flat. Equal-thickness fringes, 
running at right angles to the lines dividing the gold sandwich from the 
adjacent clear oxide layers, were observed. Their shift, across these 
lines, was a measure of the phase change due to the gold. A necessary 
second-order correction for the actual thickness-step present was made, 
using the assumed thickness of gold as the height of the step. 

All the measurements were made in collimated monochromatic light, 
incident practically normally. The form of the computations is given 
in §8. 

§5. Optica Resuirs 


Three specimens were made for optical examination, with gold of 
thicknesses 39, 71, and 103 A, sandwiched between equal films of oxide 
each 100 A thick. Annealing was limited to 300°c to avoid distorting the 
optically flat base plates. The electrical conductivities of the two thicker 
specimens were normal, for an annealing temperature of 300°c. That of 
the thinnest was low by a factor of two, and its structure must therefore 
be assumed to be in some way not quite typical. The results, namely 
v and k, the real and imaginary parts of the refractive index of the gold, 
are given in table 3. In addition, as a matter of practical interest, 7’, 
R, and R,, the transmission and front and back reflection factors, are also 
given. They are the measured values after elimination of attenuation due 
to absorption and free-surface reflection loss due to the glass support. 
Measured phase changes, not so interesting, ranged in angular measure 
up to + 25° for transmission, and between —87° and + 182° for reflection. 


Table 3. Optical Properties of three High-Conductivity Gold Films 


Film 


E Film Transmission ant 

ee resistivity | Wavelength and reflection eben 

ae PF factors : 
are 

(A) (ohm cm) (mp) 2 R, Ry v k 

436 0-643 0-192 0-134 | 2-13 1-46 

aoe 481 0-705 0-155 0-099 | 1-78 1:55 

39 17-2 10™ 546 0-776 0-104 0-057 | 1-29 210 

636 0-803 0-075 0-021 | 0-93 3-01 

436 0-583 0-177 0-091 | 1:80 1-66 

iA 481 0-663 0-138 0-060 | 1-44 1-70 

a 6-7 X10 546 0-785 0-075 0-017 | 068 2:41 

636 0-829 0-062 0:014 | 0-43 3-06 

436 0:524 0-174 0-074 | 1-79 1-61 

481 0-617 0-133 0-039 | 1-42 1-58 

103 5-6 x 10-° 546 0-763 0-071 0-002 | 0-64 2-14 

636 0-770 0-113 0-060 | 0:30 3:38 


ee ee 
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The check-thicknesses derived solely from the optical measurements are 
not given. 

They were about 20°, greater than the assumed, tabulated thickness. 
This difference could not with certainty be regarded as significant, being 
just on the limit of the range corresponding to the estimated probable 
error of the phase measurements. It is, however, in the direction which 
might be expected. 

Table 3 shows, for the 71 and 103 4 films, values of vy which agree almost 
within likely limits of error. At 39 A v is appreciably greater, by a factor 
which increases with wavelength. It is natural to associate this with the 
poorer conductivity of this film. The values of & have a moderate spread, 
but show no regular trend with thickness. 

The general agreement of the values of v and k for the two thicker films, 
coupled with the normal high conductivity, suggests that they are 
typical of the films here dealt with, in thickness of the order 100 A. It is 
therefore of interest to compare them with recorded values for bulk gold 
on the one hand, and films deposited on glass on the other. For this 
purpose, the authors have chosen from the literature data which seem 
reliable, on internal evidence. The comparison is made in terms of the 
product vk, proportional to the optical absorption (Wolter 1937), and is 
shown in table 4. 


Table 4. Values of the Product vk 


. oO —— 
Nature of gold Authority Re bd bes op, me: 
Bulk metal Meier 1910 
(interpolated) 2°20 -1:93° "1-67" e116 
High-conductivity 
film, ca. 100 A Present paper 288 2:24 L37  LOl 
100 A film on glass Goos 1937 
(interpolated) 3°22 22 2:16 5-0 


The present results are, overall, lower and more like the bulk metal 
values than Goos’. In particular, toward longer wavelengths, they show 
the steady downward trend characteristic of practically all published bulk 
values, and not a sharp upward tendency. This is consistent with the 
pale yellowish colour, as compared with the bluish-green reported by 
Goos and commonly seen. It supports the conclusion, from the electrical 
measurements, that the present films approach closely the ideal of a plane 
lamina having * bulk ’ properties. 

It seems possible, in fact, that the optical constants we have obtained 
for the present gold films are more nearly typical of bulk gold than any 
hitherto obtained by the use of thick or massive specimens. Ina massive 
specimen, only a surface layer some 200 A thick is involved in the measure- 
ments. The properties of this layer are known to be very dependent on 
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the method of preparing and polishing the specimen, but they obviously 
cannot be checked by a direct d.c. measurement of resistivity, which is 
a sensitive test of normality. For thin films, such a measurement is 
possible and affords a check on the normality of the whole of the trans- 
parent, operative, metal layer. The surface-resistivity of massive 
specimens can indeed be measured approximately by high-frequency 
methods. Where these methods have been used on massive gold, they 
have usually given results indicating an appreciably higher resistivity than 
is observed for the present films (Reuter and Sondheimer 1948), so far 
confirming, in this respect, the superiority of the films as suitable 
specimens. 

With regard to the trend of the product vk toward the longer wave- 
lengths, an upward sweep is characteristic of much reported work on thin 
gold films on glass. Corresponding film-resistances are too seldom given, 
but the present authors have failed to produce greenish 100 4 film on glass, 
with less than about four times the resistance of the yellowish high- 
conductivity annealed films of the same thickness. In resistance, as well 
as in the product vk, the 39 A film of table 3 lies (through some fault in 
preparation) about half way between the two types. A close relationship 
between optical properties and resistivity is thus confirmed, but not 
elucidated in detail. 


§ 6. OBSERVATIONS BY ELECTRON MIcROSCOPY AND DIFFRACTION 


During the electrical and optical work, attempts were made to study 
the films by transmission in the electron microscope. It was not found 
possible to sputter the films on to formvar pellicles without destroying 
the pellicle, unless this was backed by a support. Wet-stripping of the 
already coated pellicle was thus necessary, for transfer to the specimen 
grids of the microscope. That this did not alter the properties of the film 
was confirmed by electrical measurements on a large sample before and. 
after stripping. It was found very easy to overheat the films in the 
microscope, by bombardment with normal conditions. Therefore, the 
magnification, beam current, and exposure time were kept as low as 
possible. Moreover, only non-sandwich films (annealed only up to about 
200°c) could be examined because of the large electron-scatter from the 
greater total mass of material in sandwich films. With these restrictions, 
it was unfortunately found that no distinction could be observed between 
films which had, and had not, been annealed before examination in the 
microscope. The photographs obtained were thus limited to the following 
formvar-based specimens :—Figure 2 (figs. 2, 3, Plate 20, fig. 4, Plate 21), 
Gold, thickness 60 A, as sputtered ; fig. 3, Bismuth oxide, thickness 100 A, 
as sputtered; fig. 4, Bismuth oxide, 100 4 plus gold 60 A, annealed at 
200°c for 5 min; fig. 5 (fig. 5, Plate 21, fig. 6, Plate 22), Bismuth oxide, 
100 4 plus gold 60 A, heated just sufficiently to discolour it, and raise the 
resistance; and fig. 6, the same specimen as in fig. 5 after severe electron 
bombardment in the microscope. 
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In all cages the magnification in the microscope was 3750, while that of 
the photographic reproductions is 10 000 overall. 

Points to note are the similarity of fig. 2 to those for similar films 
obtained by other workers (Sennett and Scott 1950) ; the complete absence 
of observable structure in the bismuth oxide; the shaded, rather than 
granular, appearance of the gold structure of the typical film in fig. 4; 
the change of this to something much more like fig. 2 after overheating ; 
and the final separation of the gold into quite large aggregates, fig. 6, 
due to severe heating. 

Next, an attempt was made to study quantitatively the variation in 
gold thickness over the area of the gold+bismuth oxide film. Micro- 
graphs of this, and of the bismuth oxide alone, were taken under identical 
conditions of exposure and photographic processing, so as to permit of 
densitometric comparison of the photographs. The results showed 
conclusively the absence of any holes of resolvable size extending right 
through the gold layer. It was concluded, therefore, that the gold was 
continuous, but somewhat pitted or rough. 

Finally, transmission—diffraction patterns were observed. For the 
bismuth oxide only two broad haloes indicative of amorphous structure 
were found. For the annealed gold on bismuth oxide, broadened rings 
for gold were seen, together with a few faint additional broad rings 
believed to result from slight crystallization of the bismuth oxide by the 
heat-treatment. Measurements of the half-breadths of the gold rmgs gave 
an approximate value of 40 A for the dimension of the gold particles in the 
plane of the film. The relative intensities of the various (/ k 1) reflections 
were not inconsistent with random orientation of the gold crystallites. 

The chief information given by the microscope study is first that the 
yellowish high-conductivity films have definitely a more continuous 
structure than the bluish-green ones without oxide substrate. From the 
appearance of the photographs, this is thought to be due to a lesser degree 
of aggregation of the elementary particles or grains. Secondly, from the 
appearance of the overheated film, it may be concluded that a temperature 
above 200°C is required to produce aggregation of the same order as 
occurs naturally in films on a glass base at about room temperature. 
Thirdly, there is the observation that the gold crystallites are of size about 
40 A in the plane of the film and have no certain preferred orientation. 
This observation cannot be correlated closely with the inferences drawn 
from the electrical and optical properties. The particle size is, neverthe- 
less, unusually small for gold films. However, the non-sandwich film 
annealed only to 200°c is not a very suitable type of specimen, being a 
necessary compromise. The technical difficulties in getting good electron 
micrographs of the gold in sandwich films, of avoiding bombardment 
heating so as to study the true unannealed state, and of preparing suitable 
pellicle-based specimens annealed at a temperature as high as 350 or 400°c 
are being studied, with the intention of further investigations by electron 
microscopy and diffraction, 
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Lastly, the observation that the bismuth oxide is amorphous may be a 
clue to its apparent ability to control the gold structure, though it seems 
to rule out simple epitaxial growth as the relevant mechanism. The 
microscope evidence does, however, lend considerable support to the basic 
idea of the whole investigation, namely, that the bismuth oxide in some 
way limits or imposes orderliness on the mobility of the gold, and so 
prevents the kind of free grain growth—especially perpendicular to the 
film—which must not proceed far if break-up is to be avoided, and 
continuity assured, in films as thin as 100 A. 


§7. SUMMARY AND CONCLUSIONS 


It has been found, then, that sputtered films of gold as thin as 100 4 
enclosed between sputtered films of bismuth oxide, have a consistent 
specific electrical conductivity about a quarter that of gold in bulk. 
Unlike unprotected films on glass, they withstand heat treatment at 
400-500°c, after which their conductivity is three-quarters that of bulk 
gold. Their colour is pale yellowish. The optical transmission is raised 
by the heat treatment and is, throughout, higher than that of bluish- 
green films, on glass, of the same mass per unit area. The optical constants 
_ of the gold in the heat treated films are in good agreement with authorita- 
tive results of measurements on bulk specimens. 

The chief interest lies in the unusually close approach, with films as thin 
as 100 A or less, to the ideal of a plane lamina having the same constants as 
the bulk metal. Also of importance, however, are some of the conclusions 
we draw, in these circumstances, directly from observation. One is that 
electron reflection at the surfaces of the annealed films must be almost 
entirely elastic or specular, because effects related to the ratio of electronic 
mean-free-path to film thickness are only measurable with certainty in 
conditions when this ratio is upwards of about 20. Another is that the 
variation of specific conductivity with both thickness and temperature 
might find a simple alternative physical explanation in terms of granu- 
larity. affecting at one and the same time the surface roughness and the 
resistivity of the film material in bulk so that the electrical effects of these 
two factors, usually considered separately, might be related to a single 
variable, perhaps grain size. We note that the figure of 40 A found by 
electron-diffraction for the grain size in the present films, is probably the 
smallest hitherto recorded for thin gold films. For gold on formvar, 
Sennett and Scott (1950) find a rather larger grain size, with good con- 
tinuity only in films of thickness greater than about 180 A. 

Of almost equal interest are the means by which the close approach to 
bulk properties is obtained. It stands to reason that the films in the 
final state must be highly continuous, and this is confirmed in a general 
way by electron microscopy. The heat treatment seems to be annealing 
in the broadest metallurgical sense, that is to say, it probably expels 
occluded gaseous impurity, reduces lattice dislocations, and increases the 
areas of ordered structure in the gold. The function of the bismuth oxide 
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seems to be to preserve, up to usefully high annealing temperatures, 
conditions in which increase in the area of the gold would be accompanied 
by an increase, rather than a decrease of potential energy—conditions 
disfavouring the break-up of the gold into agglomerates of * droplets ’, 
both during deposition and during the later annealing (Cabrera an 
Terrien 1949). There seems here to be a field of study to which further 
work with thin films could make its own unique contribution of funda- 
mental value. 


NOTE ADDED IN PROOF 


Since this paper was prepared, electron microscope observations have 
been successfully made on a gold film 100 4 thick, between two bismuth 
oxide films and supported on a carbon pellicle, after annealing for half an 
hour at 450°c. A micrograph is shown in fig. 7 (Plate 22) for comparison 
with fig. 4 (Plate 21). A rather similar, but larger-scale, structure is 
evident. Electron diffraction clearly indicated that the size of the gold 
particles in the plane of the film is about 1000 A, and that there is some 
preferred orientation. The bismuth oxide is no longer entirely structure- 
less, but still exhibits only a comparatively small crystallite size. 

Thus, in a typical fully annealed high-conductivity film, the gold is — 
probably composed of a continuous sheet of monocrystalline lamellae, 
each of lateral dimensions some 10 times the mean thickness of the film, 
connected possibly by smaller regions in which the gold is in a greater 
state of disorder. The type of structural growth suggested in § 3 is 
thus largely confirmed. 

By the contrast it affords, the present work throws a little more light, 
perhaps, on the significance of the abnormal conductivity and optical 
constants of the more familiar green gold films, in terms of resistance and 
film-structure, but makes no serious contribution to existing theories 
(for example, Garnett 1904). It does however raise the question whether, 
particularly if the technique described could be extended to other metals, 
high-conductivity films might not be as good a medium for the study of the 
optical properties of metals as the surface layers of bulk specimens, seeing 
that the conductivity and structure of the former are open to direct 
observation while those of the latter, throughout the relevant depth of 
about 200 A or so, are much less easy to examine for ‘normality’. 


§8. APPENDIX ON OPTICAL CALCULATIONS 


Mainly algebraic calculation was used to derive the refractive index of 
the gold from the observed optical results, using the matrix notation of 
4-terminal network theory (Herpin 1947). Complex #, and H, repre- 
senting the electric and magnetic fields to the left of the optical system 
are related to those on the right by the square matrix A 


Lin ]=4 Lat} 
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For a uniform slab, of refractive index n (possibly complex, n=v—ik) 
and optical thickness 9=27nl/A 
2 sin 0 


| cos 4 
2 n 


| imsin@ cos 0 


provided the units of H and H are chosen so that, for a forward travelling 
wave in the medium, H=n#H. 

Denoting by A, the A-matrix for each oxide layer, and by A, that for 
the gold film, we have for the whole system 


vets 2 ee ee 2? @) 


To determine A ,, and thence A,, from the experimental results. it is con- 
venient to introduce the amplitude transmission factor t= \/T exp [id(T)], 
where 7’ is the measured factor and ¢(7') the corresponding 
measured phase change, and similar amplitude reflection factors 7, and 
r>. ‘The matrix # which relates the amplitudes (actually electric ampli- 
tudes multiplied by the square root of the appropriate refractive index) 
of the forward and backward travelling waves on the left of the optical 
system to those on the right may then be written 


We 1/t —Talt 
Put be yPolt 


As might be expected, the relation between H and A, involves the 
refractive indices n, and n, of the media to the left and right of the system, 


thus n 1 1 1 In 
=i {| (OG Ee Veena ey 3 
Ay 1/2 B Beetles ) 


Combining (2) and (3) we obtain 


Meat I 1 I 1/no 
iy nae | oe ood bee Pome . eae 
oe Sp “| Ny ee oad i E a4 b (4) 


A, is known from the measured thickness and index of the oxide layers, 
n,—1, and ny is the known index of the glass base, so that 4 ,, which describes 
completely the optical properties of the gold film, may be calculated 
from (4) in terms of measured t,7,; andr,. If the gold film is in fact homo- 
geneous, A, will be of the form given by (1) and it will be possible to 
calculate from its elements the index n of the gold, the optical thickness 8, 
and thence the geometric thickness 1. 

The dependence of index on thickness in (4) is only in implicit form. 
Hence, to calculate index from thickness and intensity ratios, a process of 
successive approximation was used. In this, ¢(7’), the argument of 
complex ¢, was used as an adjustable parameter. The arguments of 7, 
and r, may be obtained graphically in terms of any selected ¢(7’) from 
the following equation, which holds for a symmetrical system on a glass 
base of index 75. 


(1) 


No— 1 
A = el eee . . . . . (5) 
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A trial value of 4(7'), together with ¢(R,) and ¢(R,) found from (5), 
and also the measured intensity ratios, were inserted in (4), giving a 
calculated thickness. This was then brought into sufficiently close 
agreement with the assumed (mass per unit area) thickness by step-by- 
step adjustment of 4(7') and repetition of the calculations. In this way, 
values of 7 and k for the gold were derived from (4) in conformity with the 
measured transmission and reflection factors and the assumed thickness. 

Finally, the observed factors and phases were inserted directly into (4) 
to derive a thickness based solely on the optical measurements, providing 
the cross-check on thickness for which the phase measurements were 
mainly intended. 
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CXX. Measurement of the Elasticity of Solid Argon with an 
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ABSTRACT 


The adiabatic compressibility and adiabatic Poisson’s ratio for solid 
_ argon have been measured by determining the velocities of longitudinal 
and transverse ultrasonic waves in fine-grained polycrystalline specimens 
between 65° and 80°K. The cryostat, the ultrasonic interferometer and 
the preparation of the specimens are described. The compressibility 
is estimated to rise from 0-42 to 0-62 10-19 cm?/dyn between 0°K and 
the melting-point ; Poisson’s ratio is 0-30. Existing data are used to 
estimate the expansivity at various temperatures. The compressibilities 
are compared with other measurements and with various theoretical 
estimates. The results are consistent with a Gruneisen parameter 
y=2V/K,C,=1-2. 


§ 1. INTRODUCTION 


THE inert gases have spherically symmetrical interatomic forces 
and crystallize in close-packed structures. They have therefore been 
the subject of many theoretical investigations on the solid state which 
ideally require a knowledge of the specific heat, expansivity and com- 
pressibility as functions of temperature and pressure to verify their 
predictions. Such information has been obtained by Dugdale and 
- Simon (1953) for solid helium, but in this case the large influence of zero 

point energy makes helium unsuitable for checking theories based on 
a simple interatomic potential function (see, for example, Dugdale and 
MacDonald 1954). The heavier inert gases (argon, krypton, xenon) 
are much less influenced by zero point energy and, of these, argon is 
available to a high degree of purity at small cost, whereas krypton and 
xenon are much rarer and expensive. Argon is therefore, in many ways, 
a suitable substance for a detailed investigation of the properties of the 
solid state. 

The specific heat has already been measured and, as will be shown, 
an estimate of the thermal expansion can be made from existing data. 
This paper describes the measurement of the velocities of longitudinal 
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and transverse ultrasonic waves in solid argon in the temperature range 
65° to 80°x. A knowledge of these velocities, together with the appro- 
priate densities, enables the adiabatic compressibility, and adiabatic 
Poisson’s ratio, at zero pressure to be directly calculated. These elastic 
coefficients can then be used with the other data on solid argon, to check 
the deductions that have been made from various models of the solid 


state. 


§ 2. GENERAL EXPERIMENTAL METHOD 


Elasticity measurements, particularly of metals, have shown that only 
measurements with strain-free, homogeneous specimens yield results 
that are characteristic of the substance rather than its crystalline state. 
Ideally, for solid argon, as with other cubic crystals, the three independent 
elastic constants (C11, C12, C44) Should be measured using suitable single 
crystals. Because of the difficulty of producing and handling such 
crystals at low temperatures, we have made our measurements on 
fine-grained, polycrystalline specimens which are elastically isotropic. 

The usual static methods of measuring elasticity are not easily adapted 
to low temperatures and, in the case of the low boiling point elements | 
which have very small elastic limits, are limited in sensitivity by the 
necessity of avoiding plastic deformation. For argon, Stansfield (1954) 
has found that the tensile strength is about 50 g wt/mm? for coarse- 
grained, polycrystalline specimens at 75°K (cf. 20 kg wt/mm? for copper 
at room temperature). An ultrasonic method that does not involve 
large elastic displacements is therefore desirable. The composite resonator 
and pulse methods have both been extensively used with single crystals 
at low temperatures, but are only suitable for polycrystalline material 
when the grain-size can be made very small compared with the wave- 
length of the sound. 

A preliminary investigation (Barker, Dobbs and Jones 1953) showed 
that a compact, transparent, polycrystalline solid of grain-size of the 
order of 0-1mm (as estimated from x-ray diffraction photographs) 
could be formed by depositing the solid from the vapour on to a cooled — 
surface at a rate of growth of about 0-1mm/min. The velocity of 
longitudinal waves in such a specimen was found to be 1300 m/s at 78°K, 
so that the desired acoustic wavelength (A) of 1mm or more would 
only be produced by sound at frequencies below 1:3 Mc/s. This is too 
low a frequency for accurate measurements by the pulse method in solid — 
argon and so we have adapted the quartz crystal acoustic interferometer, 
commonly used for gases and liquids, to be suitable for the low boiling- 
point elements in the solid state. 

Kssentially, a plane wave of sound (the active area of the transmitting 
crystal is about 20 4) is propagated through an aggregate of argon crystals 
about 0-1) in diameter. This sound is reflected at a prepared, flat, 
solid-vapour interface producing the usual resonance when stationary 
waves are set up in the solid argon for each A/2 increase in thickness. 
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The velocity of sound is thus found from the known frequency of the 
quartz oscillations and the measured thicknesses of the specimen at each 
resonance. 


§ 3. THE APPARATUS 


A general outline of the apparatus, which is a development of that 
used in the earlier work, is shown in fig. 1. The solid argon is formed in 
the cylindrical Perspex cell a mounted between the liquid oxygen reservoir 
d and the quartz crystal transducer q, shown in detail in the inset figure. 


Fig. 1 


LIQUID OXYGEN 
INLET VALVES 
NOT SHOWN 


Perspex (ZQuartz Crystal S88] Tufnol 


Cryostat and ultrasonic interferometer 
(inset : quartz crystal transducer in its Perspex cell). 


The quartz crystals used were thick discs, 1-000+-0-001 in. diameter, 
and operated in thickness vibration at their resonant frequency. They 
were manufactured and gold-plated by the Quartz Crystal Co. on all 
faces, except for an annular ring (? in. int. dia.) on the back face 
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(fig. 2). This back face was mounted on a ground brass block b (fig. 1), 
surrounded by an insulating Tufnol ring. The front face was pressed 
at its rim against the square shoulder of a brass tube h by a steel spring 
held between the lower Tufnol ring r and the brass block b, both of which 
could slide along the brass rod c. In this way electrical connections were 
made to the upper face through the tube h, and to the lower face through 
the rod ¢, a phosphor bronze spring and the ground block b. At the same ~ 
time the contact between b and the quartz was sufficiently rough to _ 
prevent the sound being transmitted downwards. The springs not 
only maintained these contacts at the low temperatures of the experi- 
ments, but also enabled crystals of thicknesses from } to 6 mm to be used 
in the same holder. 
Fig. 2 


Argon gas 
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Solid Argon} Perspex 
Flattening of solid argon: (a) and (b) at second minimum ; (c) and (d) at eighth. 


The crystal holder was rigidly connected by the Perspex tube a to the 
brass cylinder d through which passed a ground brass rod, ending in the 
flat steel disc g. In this way this disc could be displaced along an axis 
perpendicular to the quartz crystal, by the micrometer M (at room 
temperature) acting through the thin stainless steel tube t. Connections 
were made using ball-bearings so that the coupling was flexible. The 
tube t, and hence the dise g, could also be rotated about a vertical axis, 
using the head R. The lower side of the disc g was acoustically flat, 
but was covered with fine radial grooves for cutting a flat face on the 
solid argon. For reflection of the sound, the flatness and parallelism 
of the sides of the argon block should be within A/4 (see, for example, 
Gold 1951), which in our case is about 0-25 mm—they were, in fact, true 
to better than 0-1 mm as checked with an optical flat. 

The complete interferometer was mounted in a vacuum enclosure — 
formed by the brass cone B, the Pyrex tube G and the glass-to-metal 
seal 8. It could be evacuated through the tube e by the high vacuum 
system V,. The Pyrex tube G was sealed with an O-ring into the brass 
top-plate D mounted similarly on the glass Dewar vessel H, which had 
two }in. wide vertical slits left on opposite sides in its silvering as 
observation windows. The connections from the tuned circuit O were 
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an earth lead E via the tube e to the tube h, and an R/F input lead through 
a pin seal in D to the rod c. Also going through D were the inlet tubes 
for the liquid oxygen valve, its control rod, and the wide tube W through 
which the Dewar vessel could be evacuated by the high-speed pumping 
system V, at a controlled rate. 


§ 4. EXPERIMENTAL PROCEDURE 


Starting at room temperature, the high vacuum system V, was used 
to exhaust to below 10-4 mm mercury the 3 1. flask F, the argon inlet tubes 
up to the needle valve N, and the vacuum enclosure BGS. An Avo 
signal generator, loosely coupled to the tuned circuit O, was then adjusted 
to resonance at the frequency of the quartz crystal, as indicated by an 
anode-bend detector joined to the tuned circuit O. (Typical circuit 
diagrams are given by Vigoureux, 1950.) The Dewar vessel D was then 
half-filled with liquid oxygen and the pumping system V, used, together 
with the liquid oxygen inlet valve, to attain the temperature required 
for the experiment. 

This was roughly indicated by the vacuum gauges py and/or P, acting 
as oxygen vapour pressure thermometers and was measured by the 
thermocouples T, and T,. These thermocouples were made from specially 
selected 33 s.w.g. eureka/40 s.w.g. copper wires supplied by the London 
Electric Wire Co. and were calibrated against a standard oxygen vapour 
pressure thermometer (Scott 1941) in the temperature range 60°-90°K 
to an absolute accuracy of 0-1°k. The thermocouple T, showed that this 
simple cryostat took up to 1 hr. to attain equilibrium and then main- 
tained its temperature within 0-2°K at the quartz crystal during the 
definitive runs of about 6 hr., while the liquid oxygen level was kept 
within } in. of that shown in fig. | by a trickle from the inlet valve. 

Before admitting the argon gas, the reservoir d was filled with liquid 
oxygen (the tubes for this also attached d to the brass cone B and are 
omitted from fig. 1 for clarity). The interferometer was then isolated 
from V, and argon gas from the cylinder A admitted through a two-stage 
regulator and the needle valve N at a very slow rate. The gas filled the 
flask F and the interferometer at a pressure, indicated by the vacuum 
gauge P,, that rose until it was a little above the vapour pressure of 
argon at the temperature of the experiment, but below the triple point 
pressure (52cm mercury). Under these conditions solid argon sublimed 
from the supersaturated vapour on the colder parts of the interferometer, 
but not on the tool g, nor in the copper inlet spiral s, since both were 
at 90°K, i.e. above the triple-point (83-8°K). After a few minutes, a 
’ transparent layer of solid argon could be seen on the quartz crystal. 
It was formed from the cooled argon gas (the reservoir d was kept full 
of liquid oxygen) and its thickness increased by about 0-1 mm/min. 

As the argon covered the quartz crystal it damped its voltage resonance, 
shown by a minimum in the anode current of the detector as read on 
a backed-off microammeter. The growth of the argon could be observed 


1074 J. R. Barker and E. R. Dobbs on 


through slits in the Dewar and its surface remained flat enough to show 
maxima and minima on the microammeter for each 4/4 increase in 
thickness. It was allowed to grow until the second minimum had just 
been passed (fig. 2 (a)), when the argon supply and the flask F (fig. 1) 
were disconnected from the interferometer, leaving only the vacuum 
gauge P, and the narrow tube e connected to Perspex cell. The argon 
quickly reached equilibrium, as shown by the gauge P, (acting as a vapour 
pressure thermometer) and the microammeter. 

The tool g was then lowered, using the micrometer M, until it reached 
the solid argon when it was alternately rotated and further depressed 
into the solid argon, subliming away first part of the ridge at its edge 
and then cutting out a cylinder 3 in. diameter in the block (fig. 2 (6)). 
Initially, the stainless steel tool was at about 90°K, but its lower surface 
was rapidly cooled in the sublimation process to be only a little warmer 
than the solid argon. When the surface was seen from the appearance 
of the turning marks of the tool to be flat, the readings of the microam- 
meter and of micrometer became consistent and were recorded. 
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Resonance in solid argon. 


As the tool slowly cut into the argon, the ultrasound, radiated as a 
plane wave from an active area of } in. in diameter of the quartz crystal, 
was reflected at the moving solid argon—vapour interface. The 
microammeter readings, plotted against the position of this surface 
from an arbitrary zero (fig. 3), showed a pronounced minimum when the 
thickness of the block was one ultrasonic wavelength. During these 
readings the frequency of the signal generator was kept within 100 c/s 
of the resonant frequency of the quartz crystal, as detected by a standard 
frequency meter. 


or 
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Having followed through the resonance curve corresponding to a 
thickness of A (the second minimum), the tool g was removed from the 
vicinity of the solid argon and argon gas was again slowly admitted to the 
interferometer. The microammeter showed the solid argon growing 
through the second minimum and then on until the third minimum had 
just been passed, when the process of subliming away the surplus argon 
and observing the resonance curve on the microammeter and micrometer 
was repeated. A number of minima were detected in this way, figs. 2 (c) 
and (d) showing the appearance of the block before and after cutting 
through the eighth minimum. The resonance curves became less 
pronounced as the attenuation of the ultrasound increased in the thicker 
argon blocks: for example fig. 4 shows a resonance curve corresponding 
to the eighth minimum of a run with a 1 Mc/s X-cut quartz crystal at 
67°K. 

This attenuation was mainly due to the scattering of the ultrasound 
at the grain-boundaries and was kept as small as possible by choosing 
a quartz crystal that produced a wavelength of 1 mm or more in the solid 
argon. It was found, however, that above 75°K the attenuation rapidly 
increased with thickness and prevented more than four consecutive 
minima being accurately observed, whereas up to ten could be recorded 
at lower temperatures. This increase of attenuation can probably be 
accounted for by recrystallization and grain-growth of the solid argon, 
for this type of attenuation is very sensitive to the size of the grains—as 
shown, for example, by Mason (1954). 


§ 5. Resutts or MEASUREMENTS 


(i) Ultrasonic Velocities 


The interferometer was used as an absolute instrument, the velocities 
being given by the product of the frequency of the oscillations and their 
mean wavelength in the solid argon. -As a check, however, it was used to 
find the longitudinal velocity in two liquids (ether and carbon 
tetrachloride) having similar velocities at room temperature to those 
observed in the argon. In each case the measurements agreed with the 
published values to within 1%. 

The argon gas was supplied by the British Oxygen Company and 
was stated to be at least 99-8% argon by volume; the main impurity 
was nitrogen with less than 2 in 10° of other gases. The solid argon 
formed from it was found to melt sharply, however, at the triple point, 
showing that the impurity level in the solid was very low. 

The accuracy of the velocity measurements in solid argon was essentially 
limited by the grain-size attenuation, which prevented a large number of 
minima being observed at the high temperatures. The wavelength was 
given by the mean of successive half-wavelength measurements from an 
arbitrary zero, thus eliminating any errors due to end-effects or thermal 
contractions of the various links between the micrometer and the argon 
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surface. The longitudinal velocities (v,), measured with a 1 Me/s X-cut 
crystal, had probable errors ranging from }% at 67°K to 2% at 79°K. 
The transverse velocities (v,) measured with a 600 ke/s Y-cut crystal, 
were associated with greater grain-size attenuation and had probable 
errors of 3°% at 67°K and 4% at 74°K. 


Fig. 5 ; Fig. 6 
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Both velocities are shown plotted against temperature in fig. 5, with 
a reasonable extrapolation to the absolute zero in each case. Smoothed 
values taken from the graphs are given in table 1. 


Table 1 
temp. (°K) i) 30 60 70 80 84 
v, (m/s) 1510 1485 1420 1375 1320 1290 
v, (m/s) 810 785 750 735 715 705 
v,/0,: 0-54 0-53 0-53 0-53 0-54 0-55 


(ii) Hlastic Coefficients 

The elasticity of an isotropic medium can be specified in terms of any 
two elastic coefficients ; those most suitable for our purpose are the 
adiabatic compressibility («,) and adiabatic Poisson’s ratio (o,). In our 
case, they are given by «,=1/p(v,,—3v7) and o,=(40,2—v7?)/(v2—v?), 
where p is the density of the solid argon. 

The density of solid argon has been found, using x-rays, to be 
1-68+-0-06 g/cm? at 20°K by de Smedt and Keesom (1925) and to be 
1-65+-0-02 at 40°K by Simon and von Simson (1924). At the triple 
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point the density of liquid argon is 1-42+-0-01 (Baly and Donnan 1902) 
and the expansion on melting (cm*/mole) has been found to be 2:96 
(Simon, Ruhemann and Edwards 1930), 3-15 (Bridgman 1934) and 
3-53 (Clusius and Weigand 1940). These give a mean density at the 
triple point of 1-608+0-015. Using the fact that the expansivity is 
zero at O°K, a smooth curve of density against temperature was drawn 
from which the three experimental points differed by less than 0-01 g/em3, 
Bes of the density and expansivity («) from this curve are given in 
table 2. 


Table 2 
temp. p 101 k, 104 107°C 10% er 
(°K) (g/em?) (em2/dyn) °5 (/°c) (*) (em2/dyn) °2 y 
() — ilexaipey/ 0-42 0:30 0 0 0:42 0:30 = 
30 ~=-:1-664 0-44. 0:31 4:3 18-4 0-44 () 3 lee 7 
60 1-636 0:49 0:30 6:8 26-8 0-51 0:30 1-26 
70 31-625 ' 0-53 0:30 75 29-3 0-56 0:29 1-19 
80 1-612 0-59 0:29 8-1 32-2 0-63 0-28 1-06 
84 1-607 0-62 0:29 8-4 35:5 0-67 0-27 1-00 
Liquid argon 


84 1-415 0-92 — 45 41-6 2-00 — — 


* (erg/mole/°c) 


Using the smoothed values of the velocities and density, a graph of 
compressibility («,) and temperature has been drawn (fig. 6). This 
compressibility and the adiabatic Poisson’s ratio (o,) are given for several 
temperatures in table 2; their experimental error does not exceed 10%. 
The specific heat (C,,) has been measured from 10°K to 80°K by Clusius 
(1936) and can be used to calculate the isothermal compressibility, 
Kp=k,+277/pC,, and the isothermal Poisson’s ratio, 


Oa Or ae L( Loe 20)(1 +o)(2?7’)]/(3xpC;,). 


These results are listed in table 2 and «7 is plotted on fig. 6. 


§ 6. DIscussION 


In comparing our measurements of the compressibilities with those of 
other workers, it is important to remember that we have measured directly 
only the ultrasonic velocities and in calculating the compressibilities 
have used data on the density (for «,) and specific heat (for «p) that have 
been measured independently. The experimental error in these data is 
about 1%, but since the measurements have been made on solid specimens 
formed in rather different ways, they may differ from specimen to specimen 
in an unknown way. Measurements of the density and specific heat on 
essentially similar specimens to those used in the present work are, 
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however, in progress at Queen Mary College and may lead to a revision 
of the following conclusions on solid argon. 

With these reservations, our results for «, and «p are accurate to about 
10°/, and may be compared (fig. 6) with the measurements of «7 made by 
Stewart (1955) using the piston-displacement. method. His condensed 
argon was a loosely packed aggregate of crystals in a cylindrical steel 
cavity of 0-8 cm3, which he compacted irreversibly at the temperature of 
condensation by rapidly applying his maximum pressure, 4000 kg/cm?. 
His curves of piston-displacement versus pressure then became reproducible 
and by fitting a curve to these measurements of 4V/V he was able to 
calculate «, as a function of the applied pressure. As shown in fig. 6, 
his zero-pressure values of x, are about 40°, higher than our measure- 
ments. Moreover, when his values are converted to adiabatic 
compressibilities («,), they are seen to lie in a line which coincides with 
«, for the liquid at the melting-point, as measured ultrasonically by 
Itterbeek and Verhaegen (1949) and Galt (1948). It therefore seems 
possible that his maximum pressure did not completely compact the 
solid, but left it less dense than a normal polycrystalline specimen. No 
measurements of the compressibility at different compacting pressures 
nor any values of the density are given in his paper. We checked our 
density at 65°K, using a volumenometer similar to that used by Baly 
and Donnan (1902) for liquid argon, and obtained 1-61-+-0-05 compared 
with 1-63 from the x-ray data (table 2). Our values of «, do indeed show 
an increase of 50°% on melting, which is not unreasonable when compared 
with the other increases on melting of 15° in volume and 25°, in specific 
heat. 

In a previous paper (Barker et al. 1953) we showed that our preliminary 
estimates of the compressibility lay within the range covered by the 
theoretical estimates of Kane (1939) and Rice (1946). Since then, 
Barron and Domb (1954) and Salter (1954 a) have extended the dynamical 
theory of the face-centred cubic lattice so that a more direct comparison 
can be made with theory. Barron and Domb show that for polycrystalline 
material where the intermolecular forces are of the Lennard—Jones type, 
Poisson’s ratio is 0-27 and v,/v, is 0-56, independent of the force constant 
when only nearest-neighbour interactions are considered. These values 
are substantially in agreement with the values of Poisson’s ratio (table 2) 
of 0-30 and v,/v, of 0-54 (table 1) derived from our measurements. Using 
an interatomic potential of the form 4(r)=4e[(o/7)!2—(a/r)*], and values 
of «, o derived from the second virial coefficient of argon gas, Barron and 
Domb obtain values at 0°K of v,=1470 m/s and «=0-38 x 10-10 em2/dyn, 
which they expect to be increased to about 1620 m/s and 0-42 10-1 
em?/dyn when distant neighbour interactions and zero-point energy are 
taken into account. We obtain 1510 m/s and 0-42 10-10 em2/dyn from 
our measurements of v, and « when they are extrapolated to 0°xK. 

Salter (1954 a), allowing for distant neighbour interactions in his 
static energy and also incorporating an approximate expression for the 
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zero-point energy, has calculated expressions for the elastic constants in 
terms of de Boer’s reduced coordinates. We have evaluated the com- 
pressibility of solid argon at 0°K from Salter’s equations and obtain 
0-52 10-1 cm?/dyn, assuming that p=1-68 (table 2). Since Salter’s 
equations involve differences between terms of the order p® and_ p%, 
they are very sensitive to the value of p adopted and the agreement 
between theory and experiment must be considered satisfactory until 
more accurate density data are available. 

Our measurements enable Gruneisen’s parameter, y=« V/x,C,, to be 
calculated, yielding the values shown in the last column of table 2. It 
will be seen that y is falling slightly as the melting-point is approached. 
The total change is, however, hardly outside our experimental error. 
It is of interest to note that the compressibility data of Stewart, although 
very different from ours, also lead to a y falling (from 1-07 at 65°K to 
0-77 at 77°K) near the melting-point. Previous estimates of y for argon 
are those of Simon and Kippert (1928), who obtained 4:5, and Rice 
(1946) who used 4-3 in his calculations. But Simon has said (quoted by 
Salter 1954 b) that he considers y=4-5 is suspect as being too high. 
Certainly our value of y=1-2 is much nearer the values for, say, the 
alkali metals (1-17 to 1-48) or the alkali halides (1-25 to 2-00). 

Our measurements are continuing and are being extended to lower 
temperatures. A detailed comparison with theory will be made later. 
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SUMMARY 


Measurements have been made of the dielectric properties of a number 
of electroluminescent phosphors in the form of powders and polycrystalline 
lumps. There appear to be two quite separate mechanisms giving rise 
to dispersion: one which is present in the dark and the other which is 
stimulated by irradiating the phosphor with ultra-violet light. The 
latter is very similar to the dispersion found in ordinary non-electro- 
luminescent powders and has not been studied further. The former, 
which seems to be characteristic of electroluminescent phosphors, produces 
dispersion at very low frequencies and has been investigated in detail. 
Measurements have also been made of the variation of dielectric properties 
with voltage and with temperature. 

To explain the effects observed it is assumed that electronic barriers 
are formed in the specimen and that motion of some of the ions in the 
barrier region takes place. On the basis of this model it is found possible 
to account in a semi-quantitative manner for the experimental results. 
Values are deduced for the height of the barrier and for the number and 
diffusion coefficient of the mobile ions. 


§ 1. [INTRODUCTION 


Since Destriau (1936, 1947) first showed that certain phosphors emit 
light when placed in a sufficiently strong alternating electric field many 
investigators have studied this phenomenon, but it must be admitted 
that the underlying mechanism of electroluminescence is still somewhat 
obscure. A number of authors have drawn attention to the interesting 
electrical properties of these phosphors. Thus, in papers by Roberts 
(1952, 1953), Zalm, Diemer and Klasens (1954), and by Butler and 
Waymouth (1954), attention has been drawn to the fact that the d.c. 
conductivity varies greatly with the applied voltage and curves have been 
given showing how the permittivity and power-factor depend on the 
applied alternating voltage at various frequencies. Some of these 
authors also noticed that the electrical characteristics of a specimen were 
different when the phosphor was kept in the dark from what they were 
when it was irradiated with visible or ultra-violet light. 


* Communicated by the Authors. 
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Many of the effects reported in these papers had been observed in this 
laboratory also, and it was therefore decided to make a more detailed 
study of the electrical behaviour of some electroluminescent phosphors. 
The results which have been obtained are reported below in two separate 
groups. The first of these is concerned with the variation of the dielectric 
properties of the phosphor with frequency and with temperature while, 
in the second, the variation with applied voltage is described. 


§ 2. MATERIALS AND APPARATUS 


Most of the measurements reported below have been carried out on an 
electroluminescent polycrystalline lump of copper-activated zine sulphide 
(sample A). X-ray examination suggested that the total number of 
crystallites in this specimen was of the order of 104. For the d.c. 
measurements, which required a different arrangement of electrodes, 
a similar polycrystalline piece (sample B) was used. The deeper colour 
of sample A suggested that it contained a higher proportion of activator 
than sample B. Some a.c. measurements also were made on sample B, 
which showed the same general behaviour as sample A. The rough 
dimensions of these samples were 0-4 0-5 x 0-14 cm. 

A few measurements were made on a single crystal of electroluminescent 
copper-activated zine sulphide of much smaller dimensions (0-2 x 0-2 « 0-05 
em), but this was difficult to handle and was fractured before it could be 
thoroughly investigated. 

Measurements have also been made on a number of electroluminescent 
powder phosphors. These were of the zine sulphide type with suitable 
activators to produce light of blue, green or yellow colour in the different 
cases. The exact composition of these commercial phosphors was not 
known. 

When powders were used, a known quantity of the phosphor was 
suspended in a solution of polystyrene in benzene and allowed to settle 
on a metal electrode about 4 cm? in area. When the mixture had dried 
a second electrode was added. For most purposes this was of tin or 
aluminium foil which was stuck to the surface of the phosphor mixture 
with a little of the polystyrene solution or with silicone oil. When it 
was desired to make measurements with the phosphor illuminated, 
the second electrode was of glass rendered conducting with a film of tin 
oxide. The thicknesses of dielectric films prepared in this way ranged 
from about 0-01 cm upwards and were measured with an accuracy of 
about 0:0005cm by means of a dial gauge. These thin films were 
particularly useful for measurements at high field-strengths. 

For the measurement of the capacitance C and equivalent parallel 
resistance R of the condenser formed by the specimen and its two electrodes 
a bridge of the type shown in fig. 1 was used. Such bridges have been 
discussed by Clark and Vanderlyn (1949) and have the great advantage 
that they are virtually unaffected by stray capacitances to earth, so that 
no Wagner earthing device is needed. With the notation of the figure, 
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the balance conditions are : 


My C _ &,k,+R,R,+R,R,+ RR, 
Hise Oy a eo nnr See 


The coils n, and n, are wound on a high-permeability core and it has 
been shown by Oatley and Yates (1954) that high accuracy can be 
achieved with the bridge even when the turns ratio is 100 or more. Thus 
the capacitance and power factor of a very small condenser can readily 
be measured. The numbers of turns 7, and n, were chosen to suit the 
range of frequencies required and, with two or three pairs of coils, it was 
found possible to cover the range from 10 c/s to 50 ke/s, though relatively 
few measurements were made below 30c/s. The balance detector 
consisted of a sharply tuned amplifier and an oscilloscope, so that no 
trouble was caused by harmonics of the input voltage. 


Fig. 1 


Bridge circuit. 


' For measurements in the frequency range from 50 ke/s to 10 Mc/s 
a standard Q-meter was used. 

With either method of measurement the voltage applied to the specimen 
did not normally exceed 3 v so that there was no appreciable variation of 
dielectric properties throughout the cycle. When the non-linearity of 
the dielectric was being investigated, the voltage across the specimen 
could be increased to 500 v. 

It has previously been suggested that non-linearity of the phosphor 
dielectric is caused by traces of moisture, but we have not been able to 
confirm this. In general our specimens were kept in a vessel containing 
a little phosphorus pentoxide, while measurements were being made on 
them. Some measurements were made, however, after the specimens had 
been heated to 200°c in a vacuum and kept there for several hours. The 
results obtained under these conditions were practically the same as in 
air; in particular, the non-linearity was as marked as ever. 
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§ 3. GENERAL SURVEY OF RESULTS ON DIELECTRIC DISPERSION 


It is convenient to express the results of measurements of the permit- 
tivity and power factor of the dielectric in terms of the components 
c, and ¢, of the complex permittivity «,;—je, and fig. 2 shows the variation 
of these quantities with frequency for specimen A when irradiated with 
ultra-violet light. In the curve for «, the large values at very low 
frequencies followed by a fall and then a peak at about 15 ke/s suggest 
that two quite separate mechanisms are operating simultaneously, and 
this is confirmed by the two points of inflexion in the curve for ¢. 


Fig. 2 
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Frequency (¢/s) 
Dispersion curves for sample A irradiated with ultra-violet light. 


When the measurements were repeated with the specimen in the dark, 
the curves of fig. 3 were obtained and there is now no sign of any anomaly 
at 15 ke/s. If we assume that the two mechanisms postulated above 
operate independently of each other, and that the one centred on 15 ke/s 
is caused by irradiation with ultra-violet light, it is legitimate to subtract 
the values of fig. 3 from those of fig. 2 to give the curves shown in fig. 4, 
which are now characteristic of the single mechanism caused by irradiation. 

Returning to the curves of fig. 3, it is clearly of interest to know what 
behaviour would be encountered at still lower frequencies. With the 
present form of apparatus and with the relatively low capacitances 
(<100 pF) which were to be measured, it was not found possible to work 
at frequencies below 10 c/s. However, if the dispersion is of the Debye 
type, there should be a peak in the curve for e, and the frequency at 
which this peak occurs should be greater at high temperatures than at 
low. Some measurements were therefore made on sample B and on an 
electroluminescent phosphor in powder form at a temperature of 300°c, 
but even at this temperature the values of ¢, and ¢, were rising steadily 
at the lowest frequencies at which measurements could be made. 


Electrical Properties of Electroluminescent Phosphors 1085 


bo—o- Experimental curve 
Theoretical points 


(a) 
Frequency (c/s) 


»—e- Experimental curve 
+ Theoretical points 


(0) 
Frequency (c/s) 
Dispersion curves for sample A in the dark. 
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§ 4, Discussion OF RESULTS ON DISPERSION 


The curves shown in fig. 4 are clearly of the general type that can be 
represented by the Debye equations (Frohlich 1949). If the dispersion 
could be ascribed to a process governed by a single relaxation time, the 
maximum value of «, should be given by (¢«,—€.)/2, where «, and «,, 
are the values of ¢, at zero frequency and infinite frequency respectively. 
In fact the maximum value of ¢, is rather less than this, which suggests 
that there is a small spread of relaxation times. A more detailed analysis 
by means of a Cole—Cole plot (Cole and Cole 1941) confirms this suggestion. 


Fig. 4 


Results obtained by subtracting the values of fig. 3 from those of fig. 2. 


Results of the type shown in fig. 4 are well known for ordinary non- | 
electroluminescent phosphors when stimulated with ultra-violet light 
(Garlick and Gibson 1947) but in such cases the relaxation times have 
usually been an order of magnitude smaller than the mean value of 
10-5 sec found from fig. 4. These results will not be considered further 
since they do not appear to be connected with the phenomenon of 
electroluminescence. On the other hand ordinary phosphors do not 
exhibit any marked dispersion when they are not stimulated by radiation. 


Thus the remainder of this paper will be concerned principally with 
results obtained under these conditions, 


Electrical Properties of Electroluminescent Phosphors 1087 


Turning now to the curves of fig. 3, since no peak value has been found 
for «, in the range of frequencies investigated it is not immediately 
obvious whether the results can or cannot be represented by the Debye 
equations. To test this point, use was made of the equations 


€9°=(e,—€,)(€y—€ x), Pi Clee ee ee 9 
=(e€,—€,)/(€y—€~) «SMe ce ae Co) 
where 7 is the relaxation time and the values of ¢, and «, correspond to 
angular frequency w. These equations should hold for any substance 
to which the Debye theory applies and for which there is a single relaxation 
time (Frohlich 1949). Taking values from fig. 3, putting ¢,—=10 
and using eqn. (2), <«, was calculated for frequencies of 1 ke/s, 100 ¢/s and 


30 c/s. These values of «, were then substituted in eqn. (3) to give values 
of 7. The results are shown in the table and it is clear that there is no 


if E, T 
1 ke/s 88 3-4 x 10-4 sec 
100 c/s 137 Zk XL Os Bec 
30 c/s 307 5-6 10- sec 


consistency in the values of either «, or tr. Furthermore, the values of 7 
are incorrect since they would correspond to a peak of ¢, in the range of 
frequencies covered by fig. 3. It is safe to conclude therefore that Debye 
theory with a single relaxation time cannot account for the experimental 
results. It is less easy to be certain that agreement between theory and 
experiment could not be reached by postulating a suitable spread of 
relaxation times, but such an explanation does not seem particularly 
plausible. 

Many authors have attempted to explain the dielectric properties of 
inhomogeneous materials in terms of a system containing particles which 
behave as capacitances separated by material of high, but not infinite, 
resistivity. It is easy to show that a simple system of this type containing 
only one capacitance in series with one resistance gives rise to dispersion 
which can be represented by the Debye equations, so the comments 
made in the foregoing paragraph apply to this model also. To explain 
the curves of fig. 3 it would be necessary to postulate a very complex 
system in which a wide spread of relaxation times was operative. If, 
however, the capacitances are assumed to vary with frequency, it is 
possible to construct a much simpler model and this, in effect, is what is 
done below, the capacitances being ascribed to the existence of barriers 
either between crystallites or at the surfaces of the specimen adjacent to 
the electrodes. 

Whenever barriers are formed in this way and an alternating electric 
field is applied to a specimen which has appreciable conductivity the 
local accumulation of charge can contribute to the measured capacitance 
and power loss of the specimen and may swamp all other factors. 
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Conditions in the barrier region will be determined by the numbers of 
the different types of carrier, by their mobilities and rates of formation 
and recombination, by the ease with which they can pass to the electrodes 
and by many other factors. Progress towards a general solution of the 
problem has been made by Jaffé (1933, 1952) and by Macdonald (1953) 
but their formulae are very complex and involve assumptions which 
may or may not be valid in the present experiments. In the first instance, 
therefore, our experimental results will be compared with those predicted — 
by a simplified treatment due to Volger, Stevels and van Amerongen 
(1953). These authors assume that only one type of conducting carrier 
is present in the specimen and that carriers of this type cannot pass from _ 
the specimen to the electrodes. Under these conditions the type of 
barrier formed will depend on the ratio of the conductivity o to the angular 
frequency of the applied electric field. If w<4mo/e where « is the 
permittivity of the medium itself when no barriers are formed, the 
build-up of charge near the electrodes will be governed almost entirely 
by space-charge and the barrier will be of the well known Schottky 
type. Such barriers lead to a surface capacitance which is independent 
of frequency and cannot therefore account for our experimental results. 
If, on the other hand w>47a/e the quantity of charge which can build 
up near an electrode during the half-cycle before the direction of the 
applied field is reversed, is so small that space charge is negligible. The 
accumulation of charge is then opposed only by diffusion which is brought 
into play by the concentration gradient. If x is the distance from the 
electrode in a one-dimensional system and n the concentration of charge 
carriers, the relevant equation is 
on 07 
aoe a Uf 
where D is the diffusion coefficient. At the electrode the total current 
is by hypothesis, zero, so 


on 


Solving (4) and inserting the boundary condition (5) we find that, when 
the applied electric field is given by H=E, sin wt, the steady-state value 
of n is 


ok, (a0) F ves ] wW 
N=Ny+ eae exp | -« JG) sin | ot Teste Ja) . (6)% 


where 7, is the carrier concentration when there is no field. 

Volger, Stevels and Amerongen then show that, in addition to the 
normal capacitance C and equivalent parallel resistance R of the specimen, 
there is now a barrier capacitance C', and parallel resistance Ff, such that 


(a 
R= |(Z) [o= a lie a 
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where 6= 4/(D/2w) may be taken as a measure of the thickness of the 
barrier. 

Since a barrier of this type will be formed at each electrode, the 
equivalent circuit for the complete specimen will be as shown in fig. 5. 
Both C, and R, tend towards zero as w gets very large, so C and R can 
be determined from measurements at high frequencies, and C, and R, 
can then be calculated at any other frequency when the total impedance 
of the specimen is known. 

A comparison of this equivalent circuit with the experimental results 
shows that the model gives a good qualitative picture of the variation of 
€, and «, with frequency but is wrong in certain particulars. For 
example, at very low frequencies, eqns. (7) and (8) taken with the complete 
equivalent circuit suggest that «, should be greater than ¢, whereas it is 
actually only about half as great as «,. Furthermore, as will now be 
shown, when numerical values are inserted into the equations, it becomes 
clear that the model so far considered is too simple. 


Fig. 5 


Simple equivalent circuit of specimen. 


The variation of <«, and «, with frequency (fig. 3) indicates that the 
barrier must be of the type in which space-charge is negligible and the 
variation of charge is governed by diffusion. But this implies that 
4na/e<w even at the lowest frequency of measurement which was 10 c/s. 
Putting « equal to 10, this means that o must be less than 10 e.s.u. or, 
converting to practical units, that the resistivity of the material must 
exceed 10410hmcm. Such a value is several orders of magnitude 
greater than the value (about 10¢ohm cm) which can be deduced from 
the equivalent parallel resistance of the specimen at frequencies above 
100 ke/s. The high-frequency value of the equivalent parallel resistance 
might, of course, be influenced by other sources of dielectric loss and in 
that case it would bear no relation to the resistivity of the material. 
However, if the resistivity really were as great at 101 ohm cm, the bulk 
resistance of the specimen would be so great (R in the equivalent circuit 
of fig. 5) that it would be quite impossible to account for the measured 
values of impedance by any barrier effects. 
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Consider next the possibility that a thin layer of the specimen adjacent 
to each electrode (or at the surface of each crystallite) has a resistivity 
of the order of 10!1ohmecm, but that the bulk resistivity is about 
10°ohm cm. The equivalent circuit would then be as shown in fig. 6 (a) 
where C', and R, are the barrier constants as before, while C, and R, 
are the normal capacitance and resistance of the high-resistance layer. 


Fig. 6 
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Improved equivalent circuit of specimen. 


This is a much more promising model which represents the experimental 
results very satisfactorily provided R, makes a negligible contribution 
to the total impedance. Assuming for the moment that this is the case, 
the equivalent circuit is as shown in fig. 6 (b). 

As before, the values of C and R can be obtained from measurements 
at frequencies above 100 ke/s. Furthermore, if we assume the resistance 
of k, to be equal to the reactance of (, and each to be equal to \/(2D/w)/ce 
as theory indicates, the values of \/(2D)/o and C,, can be found from 
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measurements of «, and e, at one other frequency, which may con- 
veniently be low enough for the reactance of C to be very large in 
comparison with k. Then, assuming C’, to be independent of frequency, 
the magnitudes of C, and R, at any other frequency, and hence the 
overall values of «, and «, can be calculated. In figs. 3 (a) and 3 (b) 
a number of values of «, and ¢, calculated in this way are plotted. The 
values of \/(2D)/c and C, used were chosen to give the best overall 
agreement with the experimental curves, and are indicated in fig. 6 (0). 

The assumption that R, is negligibly small raises a difficulty which 
has not been completely resolved and which must now be discussed. 
Consider first the case where the high-resistance layer consists of homo- 
geneous material of composition different from that of the bulk. Since 
our specimens were activated by the diffusion of copper inwards from 


Hig. 7 


Schematic representation of barrier. 


the surface, it would not be surprising if the resistivity of the surface 
region were different from that of the bulk of the material. However, 
with resistivity p and thickness d, the resistance of such a layer of unit 
cross section would be d . p, while its capacitance would be </47d. It then 
follows that, at all relevant frequencies, the reactance of this capacitance 
would be small compared with the resistance of the layer. In terms 
of the equivalent circuit of fig. 6 (a) this means that, if R, is negligibly 
small in comparison with R,, the impedance of R, and (, will be effectively 
short-circuited by the reactance of C,, and this is contrary to the experi- 
mental results. 

Following Jaffé (1933, 1952) it may be assumed that both positive 
ions and electrons (or, alternatively, negative ions and holes) are present 
in the material but that, in the vicinity of an electrode, a rectifying 
barrier is set up. Using the conventional representation of such a barrier, 
as in fig. 7, there will be a region of the material of width A from which 
free electrons have been largely drained and this region can thus con- 
stitute the high-resistance layer which has previously been postulated. 
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As is well known, the capacitance of this barrier layer corresponds to a 
thickness of dielectric but, when calculating its resistance, allowance 
must be made for the fact that the resistivity increases rapidly from 
Ato B. This occurs because electrons from the bulk of the material can 
diffuse into the barrier region, but the concentration of such electrons 
decreases very rapidly with distance to the left of A. Assuming the 
conductivity at any point to be proportional to the electron concentration, 
a rough calculation suggests that the resistance of the layer would be 
about 0-O05ppA per cm2, where p, is the resistivity at the electrode. The 
ratio of resistance to capacitive reactance is less by a factor of twenty, 
therefore, than in the case of a homogeneous barrier previously considered. 
This factor is still inadequate to make R, (fig. 6 (b)) negligible ; 
nevertheless it seems probable that the picture presented by Jaffé’s 
theory is essentially correct. In view of the inhomogeneous nature of 
the specimen, it would not be difficult to introduce further plausible 
assumptions which would have the effect of reducing RF, to a sufficiently 
low value. In particular, one might postulate a non-uniform distribution 
of immobile ionized donors which would concentrate a large part of the 
voltage drop into the immediate vicinity of the electrode without altering 
the width of the layer from which free electrons had been drained. This 
would have the effect of reducing Rk, without appreciably changing C,. 
However, the present measurements throw no light on the mechanism 
which makes Ff, negligibly small, so further speculation is unprofitable. 

The representation of the behaviour of the specimen by means of an 
equivalent circuit cannot, of course, be more than a crude approximation 
to the truth. At the same time it does provide a rough physical picture 
of what is taking place and the separation of the various effects, which is 
inherent in it, makes comparison of theory and experiment much simpler 
than it otherwise would be. The measured values of the components 
of the equivalent circuit, which are tabulated in fig. 6 (6), will now be 
used to calculate various properties relating to the charge carriers. 
In view of the approximation introduced by the use of an equivalent 
circuit, the results will probably not be very accurate, but they should be 
correct as to order of magnitude. ; 


§ 5. CALCULATION OF CONSTANTS 


It has so far been assumed that only two barrier regions exist in the 
material but, with a polycrystalline specimen, this may not be the case. 
Let us assume that there are S pairs in series. Then the true magnitudes 
of C, and C, will be S times as large as the values given in fig. 6 (b), 
while the value given there for R,, will have to be divided by S. However, 
with these new constants, the overall variation of «, and ¢, with frequency 
will be the same as before. 

From the magnitude of C, and taking the permittivity of the material 
to be 10, the value of A the width of the layer from which electrons have 
been drained is found to be 10-2/S cm. This region will be referred to 
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as the electron barrier. Within it, at any rate near the electrode, the 
ionic conductivity must be less than 10-1! (ohm cm)-! if diffusion of 
ions is to produce an additional impedance with the correct frequency 
variation. Furthermore, electronic conductivity in this region must be 
still lower if it is not to short-circuit the effects produced by the ions. 

Substitution of this maximum value of the ionic conductivity and the 
measured value of C’, into eqn. (7) gives for the diffusion coefficient 


te Dero oC SCG lem eee. Be a (9) 


Similarly, from eqn. (8), the maximum value of 5, the width of the 
region where variation of ion density is appreciable, is found to be 
3X 10~4/S cm at 100 ¢/s and 3x 10~5/S cm at 10 ke/s. This region will 
be referred to as the ion barrier and it is to be noted that it occupies 
only a small fraction of the width of the electron barrier. This fact, 
together with the original assumption that variation of ion density does 
not produce appreciable space-charge, provides justification for the 
separation of the effects due to ions and electrons respectively and their 
representation by components in an equivalent circuit. 

If the ions in question were normal constituents of the crystal moving 
by way of lattice defects, it has been shown by Mott and Gurney (1940) 
that the ratio o/D should be of the order of 3x10, whereas our 
measurements give a value of about 10-®. It therefore seems more likely 
that we are concerned with the diffusion of interstitial ions. For these 
the Einstein relation 

Ce ee OE es. coe ALO) 


where v is the mobility, should hold and this leads to a maximum value 
of 4x 10~3/S? cm/sec per v/em for v. Finally if n) is the concentration of 
interstitial ions 
Ca Otero FY etd) 
and this gives 
Mee lO. CMe eee ss fe. Ce (TS) 


These values of D, v and n, are all reasonable for interstitial ions if S 
lies between 1 and 100. For a single crystal S would presumably be 
unity, but with polycrystalline specimens such as were used in the present 
measurements it is not improbable that sufficient of the activator copper 
would diffuse along the grain boundaries to set up barriers there. A rough 
determination of the crystallite size by means of x-rays showed that, 
in samples A and B, the number of grain boundaries in series between the 
electrodes was of the order of 10. 


§ 6. D.C. MEASUREMENTS AND THE VARIATION OF ELECTRICAL 
PROPERTIES WITH TEMPERATURE 


The d.c. resistances of the polycrystalline specimens were so high that 
reliable measurements of them could not be made without the use of 
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guard-ring electrodes. It was found possible to deposit suitable silver 

electrodes on sample B by vacuum evaporation and the results of d.c. 

measurements at room temperature are shown in fig. 8. The corres- 

ponding resistance per unit area for small applied voltages is about . 
21011 ohms, which is consistent with the upper limit for electronic 

conductivity deduced above. ) 


Fig. 8 


Current density (amp/cm?) 


d.c. volt 
Dark current/voltage curve for sample B. 


Measurements were also made of the variation of the low-voltage d.c. 
resistivity p with temperature and a semi-log plot of the results is shown 
in fig. 9. From this it appears that the change of d.c. resistance is to be 
associated with an activation energy of about 0-5 ev, and it is reasonable 
to identify this with the height of the rectifying electronic barrier. 

This height, Vp, is connected with the concentration n, of ionized 
donors within the electronic barrier and the width A of the barrier, by 
the well known expression 


A=[V pe/2anje|Y*: eo ee 
Taking for A the value 10-*/S cm found previously, this gives 
N= 2X10 Stems. 6 As 


Comparison of eqns. (12) and (14) shows that it would not be inconsistent 
to identify the ionized donors with the mobile interstitial ions. However, 
there is no reason why n, and », should be equal, since some of the donors 
may be immobile. Furthermore, if the distribution of donors is not 
uniform, their concentration in the ionic barrier may differ from the mean 
value throughout the electronic barrier. 
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Since the electronic barrier region occupies an appreciable part of the 
total volume of the specimen, it is reasonable to take n, as the concentra- 
tion of free electrons in the bulk of the material. 


The bulk conductivity 
o,, Will then be given by 


Ono = 6U 1 oer (15) 
where w is the electronic mobility. Taking 10-§ (ohm cm)! for Ou 
this gives 

u=30/S? cm/sec per v/em. (16) 


This is a reasonable value if S is equal to unity, but rather low if 8 is of 
the order of 10. 
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103/°K 
Temperature variation of d.c. resistance for sample B. 
From the diffusion theory of rectification at an electronic barrier, 


the limiting value p,, for the resistance per unit area in the reverse direction 
with large applied voltage should be given by 


1/po=(meu/A) exp (—eVp/kT). . . . » ~ (17) 
Inserting the values already found for u, n,, A and Vp, this gives 
Oe 62) Cen er ee ee ee ee LS) 


This is about one hundred times as great as the measured value, but 
such discrepancies are commonly found in investigations on rectifiers 
and have been attributed to field-emission of electrons from the metal 
into the semiconductor. In any case, bearing in mind the approximate 
nature of our calculations, the discrepancy is not serious. 
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Measurements have also been made of the variation of a.c. impedance 
with temperature for sample A and from these, using the equivalent 
circuit of fig. 7, the variation of C, with temperature has been calculated. 
This time the semi-log plot, shown in fig. 10, leads to an effective activation 
energy of 0-15 ev. According to eqn. (7) C, is proportional to o/,/D and 
both c and D may be expected to vary with temperature. If the number 
of interstitial carriers is assumed to be constant, o will change only because 


Fig. 10 
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103/°K 
Temperature dependence of barrier capacity (sample A). 


the mobility v changes and, since v is proportional to D, (eqn. 10), C, 
should then be proportional to \/D. In this case the activation energy 
found from fig. 10 will be half the height of the potential hill to be sur- 
mounted by an interstitial ion in going from one position to an adjacent 
one, viz. 0-30 ev. If d is the distance between these two positions and v 
the frequency with which the ion can vibrate, D should be given as to 
order of magnitude (Mott and Gurney 1940) by 


D=dVexp (sei) Pen a eee) 


eae U is the height of the potential hill. Putting »v—10!3 and 
d2—= 10-15 om. we find for the order of magnitude of D the value 
10~7 cm?/sec, which is consistent with the value given by eqn. (9). 
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Fig. 11 
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Dependence of complex permittivity on applied voltage. (Sample A in 
the dark.) 
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§ 7. Non-Linear EFFECTS 


It has previously been reported that the electrical behaviour of an 
electroluminescent powder is a function of the strength of the electric 
field applied to it and, in the present experiments, an investigation of 
these non-linear effects was carried out. As before the results were 
found to depend on whether the specimen was kept in the dark or 
irradiated with ultra-violet light, and in figs. 11 (a), 11 (6) and 12 are 
curves showing the variation of ¢, and «, in the two cases for sample A. 
Since it has already been shown that two quite separate mechanisms are 
involved, it is legitimate to subtract the ‘ dark’ from the ‘ light ’ values 
to give the curves of fig. 13 for the mechanism which depends on stimula- 
tion by ultra-violet light. 


Fig. 12 
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Variation of complex permittivity with applied voltage (sample A). 


Curves very similar to those of fig. 13 have been obtained with certain 
ordinary phosphors, which are not electroluminescent, when these were 
irradiated with ultra-violet light. There is no reason to suppose, therefore 
that this part of the non-linearity is connected with electroluminesoenee 
and it will only be discussed briefly in the present paper. 
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It has been shown in § 4 that the mechanism brought into play by 
light-stimulation conforms to the Debye equations. A model which is 
often invoked in such cases is one in which dipoles (or single charges) in 
the solid have two or more positions of equilibrium separated by potential 
hills. The discussion of this model (Fréhlich 1949) usually assumes that 
Ey <xT, where uw is the dipole moment and £ the applied field, but if this 
is not the case it can be shown that, to a second approximation, 


ea a ae) yay eet 
ae : Q+¢ oe ), eee ie )) 


leer" 4x27? 7224 1 
(€,—€,.)Tw (1+ pH? oS) | 


Aa 21 
“2 1+ 77w? An? T? 7241 ep 
If Eu is comparable with «7’, still higher terms will become important. 


Fig. 13 
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Result of subtracting the ‘ dark ’ from the ‘ light ’ values in fig. 12. 


The curves of fig. 13 were taken at a frequency of 1 ke/s, for which 
wt <1, so eqns. (20) and (21) predict that both ¢«, and e, should decrease 
with increasing field-strength, in agreement with the experimental results. 
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A more stringent test of the equations would be to see whether the 
non-linearity changes sign at sufficiently high frequencies. This was 
not possible with the existing apparatus since the Q-meter with which 
the high-frequency measurements were made did not provide a suffi- 
ciently high field-strength. The matter is therefore left for further 
investigation at a later date. 

Returning now to the curves of figs. 11 (a) and (6), consider the effect 
of applying relatively large alternating voltages to the specimen in the 
light of the mechanism which has been invoked to. explain the low- 
frequency dispersion in the dark. In the theory leading up to eqn. (6) 
it was assumed that the conductivity o of the high-resistance layer is 
a constant. However, this will no longer be true if the change in 
concentration of carriers in the barrier region, brought about by the 
alternating field, is an appreciable fraction of the initial concentration 7p. 
Thus, with sufficiently strong fields, c in eqn. (6) will vary throughout 
the cycle and the barrier impedance will become non-linear. An 
analytical solution in this case would be extremely complicated and will 
not be attempted. It is to be noted, however, that changes in o are in 
step with changes in n so that we might expect the variations of n to be 
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larger than they would have been if o had remained constant. This 
implies that C’, should increase with increasing applied voltage while R, 
should decrease, as is found to be the case. In this context C, and R, 
are mean values averaged throughout the cycle. Furthermore, eqn. (6) 
shows that the variations of n are inversely proportional to s/w so it 
might be expected that the non-linearity of C, and R » (measured by 
changes from the small-signal values) would vary with frequency in 


: 
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roughly the same manner. The experimental curves of fig. 11 show 
that this is the case. 

Substituting in eqn. (6) the value of o/,/D previously deduced, it is 
easy to show that the variation of ionic concentration throughout the 
cycle at, say, 100 c/s will be comparable with the mean concentration 
when the voltage applied to specimen A is a few tens of volts. Under 
these conditions non-linearity would be expected and is, in fact, observed. 

So far, only variations of C, and R, have been considered, but (Obs 
also will depend on the applied voltage. The width of the electronic 
barrier is governed by space charge, so will increase and C,, will decrease 
with increase in applied voltage. The changes in C,, will thus be in the 
opposite direction to those in C’, and to those observed experimentally. 
It must be assumed therefore that changes in ©, predominate. The 
fact that electronic barriers alone would lead to a variation of capacitance 
with applied voltage of the wrong sign, provides additional evidence for 
the existence of effects caused by mobile ions. 

Measurements have been made of the impedance of sample A using 
a small alternating voltage of about 1 v superimposed upon a d.c. bias 
which could be varied over a wide range. Some results of such measure- 
ments are shown in fig. 14, but this case has proved very intractable 
from a theoretical point of view. It can easily be shown that when 
sufficient bias is applied to produce an appreciable change of impedance, 
ionic space-charge must be of great importance. Under these 
circumstances eqn. (6) is no longer valid. 


§ 8. CONCLUSIONS 


A great many measurements have been made on electroluminescent 
powders of various kinds and in all cases the results were qualitatively 
similar to those described above. In general the powders exhibited less 
marked dispersion than the polycrystalline specimens, but the two 
separate regions of dispersion were nevertheless present. It has been 
shown that the high-frequency dispersion is very similar to that which 
occurs in phosphors which are not electroluminescent (Garlick 1949), 
though the relaxation time is usually about one-tenth as long in these 
latter cases. This difference may or may not be significant. 

The low-frequency dispersion has not been reported in phosphors 
which are not electroluminescent and it is therefore natural to attempt 
to correlate it with the phenomenon of electroluminescence. To explain 
this phenomenon other investigators have invoked the presence of 
rectifying electronic barriers (Zalm e¢ al. 1954, Piper and William 1954) 
and such barriers form an essential feature of the model which has been 
put forward to account for the present results. By themselves, however, 
they cannot readily be used to explain the observed low-frequency 
dispersion and we have therefore been led to postulate the existence of 
mobile ions which cannot pass to the electrodes. The resulting model 
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is selfeonsistent and provides a satisfactory explanation of a variety 
of experimental results. 

Since electrcluminescence is not quenched when a specimen is placed 
in a magnetic field of the order of 10° gauss (Ince 1954), it seems necessary 
to assume that electric fields present in certain parts of the phosphors 
are much greater than would be expected from the applied voltage and 
the thickness of the specimen. It is possible that the necessary iInhomo- 
geneity of field is brought about solely by the presence of electronic 
rectifying barriers. (The width of these barriers has already been discussed 
and the voltage drop across them at any frequency can be calculated from 
the equivalent circuit of fig. 6 (6). It appears that the field within the 
barrier might be large enough to account for electroluminescence. 

This leaves unanswered the question whether ionic mobility plays an 
important part in the phenomenon of electroluminescence. Evidence 
for the presence of mobile ions in electroluminescent phosphors has been 
given, but the effects ascribed to these ions would probably not have 
been observed if electronic barriers had not also been present to reduce 
the electronic conductivity to a very small value. For this reason ionic 
mobility may have escaped observation in phosphors which are not 
electroluminescent. Nevertheless, experiments such as those of 
Waymouth and Bitter (1954) in which an electroluminescent phosphor was 
subjected to pulsed fields, suggest that ionic mobility might provide an 
explanation of the very slow build-up and decay of light output which 
were there observed. 
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ABSTRACT 

The attenuation of ultrasonic waves by conduction electrons is analysed 
in terms of the free-electron model of a metal. At such low frequencies 
that the electronic free path is smaller than the wavelength, the result 
obtained agrees with that found by other authors using simpler concepts, 
but this agreement breaks down when the free path becomes comparable 
with or larger than the wavelength. At very high frequencies the attenua- 
tion of longitudinal waves increases in proportion to the frequency without 
limit, while that of transverse waves tends to a constant value. The 
drop in attenuation observed when a metal becomes superconducting is 
shown to present quite serious problems of interpretation, and tentative 
suggestions are made as to possible explanations. 


$1. INTRODUCTION 


RECENT experiments by Bommel (1954) and Mackinnon (1955) have 
revealed that ultrasonic waves propagated through pure metals suffer an 
attenuation which increases at the lowest temperatures towards a 
maximum at the absolute zero. If the metal becomes superconducting 
the attenuation falls sharply, but apparently not discontinuously, below 
the transition temperature. The effect in normal metals has been dis- 
cussed along similar lines by Mason (1955), Morse (1955) and Kittel 
(1955), the first two in particular giving a numerical analysis of Bommel’s 
results for lead and tin which shows that the explanation advanced by 
these authors is substantially correct. If the electrons in the metal are 
regarded as a degenerate Fermi gas enclosed in a box which is subjected 
to periodic distortion, then in the absence of collisions the Fermi surface 
will be adiabatically distorted, a spherical surface for example becoming 
ellipsoidal. Collisions of the electrons with the ionic lattice tend to restore 
the surface to its spherical shape, but the process is never complete, 
particularly if the relaxation time is long. In consequence the total 
electronic energy is on the average somewhat greater than its equilibrium 
value, and there is a steady dissipation of mechanical energy as heat 
through the influence of collisions. The rise of attenuation with lowering 
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temperature is then readily interpreted as due to an increase of relaxation 
time, as is directly indicated by the behaviour of the electrical resistance. 

This view of the process, though it leads to the correct answer at 
sufficiently low frequencies, is nevertheless deficient in some respects. 
In the first place it suggests that the attenuation should vary quadratically 
with frequency up to such frequencies that wr is comparable with unity, w 
being the angular frequency and 7 the relaxation time of the electrons. 
However, as was pointed out to the author by Dr. Morse in conversation, 
on account of the Fermi velocity being much greater (by a factor of 
about 300 in most metals) than the velocity of sound, long before wr 
reaches unity, the electronic free path has become comparable with the 
wavelength of the sound wave, and it is then no longer permissible to 
isolate for consideration a portion of the metal small in comparison with 
the wavelength. Secondly, the assumption that deformation of the ionic 
lattice leads automatically to a change in the distribution function of the 
electrons, as if they were enclosed in a box, does not accord with the true 
situation. In reality it is the combined effect of collisions and electric 
fields which brings about this change. The principal object of this paper 
is to give a solution of the problem which takes these points into account, 
for the ideal case when the electrons form a free degenerate gas with 
spherical Fermi surface. It will then be possible to see more clearly what 
problems are raised by the behaviour of superconductors, and suggestions 
will be made as to the nature of the new mechanisms operating under 
these circumstances. 


§2. LonerrupInaL WAVES 


Consider an infinite block of metal carrying a plane longitudinal wave 
in the z-direction, of frequency w and wavenumber k. All variables 
associated with the wave may then be understood to be multiplied by 
exp {i(wi—kz)', so that the differential operators take the form d/dz= 
o/ey=0, o/dz=—ik, d/et=iw. The wave velocity v, is of course w/k. 
In the undistorted metal let there be N free electrons per unit volume, 
occupying all velocity states up to v9, the Fermi velocity. During the 
passage of the wave, which may be characterized by a particle velocity w, 
the density of electrons will not necessarily remain constant but will have 
a periodic increment n. If the densities of electrons and ions do not 
keep exactly in step there will be space charges leading to a periodic 
electric field @ in the z-direction. We now calculate what distortions the 
Fermi surface suffers as a consequence of the combined effects of & and 
of collisions. For this purpose we adopt the simple kinetic device of 
following a single electron having velocity v, through the lattice ; as a 
result of & and collisions it acquires an additional velocity v9’ in the 
direction of its motion, and this value of v,’ is taken to give the displace- 
ment of the Fermi surface in this particular direction. 

An electron moving with velocity v, at an angle @ to the z-direction 
experiences an electric field which oscillates at a frequency w(1—v, cos 6/,), 
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and since in a field & its change of velocity is given by the equation 
dv, |dt=e& cos 6/m, it follows that due to & alone, collisions being neglected, 


iv[09'],(1—% cos O/v,)=eF cosB/m. . . . . - (1) 


To take account of collisions we assume that v9’ is continuously relaxing, 
with relaxation time 7, towards the value which it would have if the 
electron gas were locally in equilibrium. This equilibrium value of 
v9 is not zero for two reasons : first because the local electronic density is 
N-+n, not NV, and secondly because the lattice is in motion with velocity 
u. The velocity distribution in local equilibrium thus corresponds to a 
Fermi sphere of radius vp(1-+-4n/N) onto which is superposed a uniform 
velocity w. And the local equilibrium value of v9’ is thus given by the 
Teale U6 ‘oq ENN Fu. COS 5 a ee 
We may therefore introduce the effect of collisions by writing 

[dvp’/dt]on=%o[%p Jeou(1—U%p COs 4/v,)=(FNUp/N +u cosO—vq')/7. - (3) 
Hence, combining (1) and (3) 

Pee i any (e&r/m+w) cos 6+-4nv_/N 
Y% =[U ]-+[% Jeon= l+iwr—ikl cos 6 : . : (4) 

where / is the electronic mean free path, defined as vor. 


From (4) which gives the displaced form of the Fermi surface, the 
electronic current density J, may be immediately calculated, since 


5 : 3Nef Nevo f 
Jl =3 — ‘ef Th 0 
. Ne | sin 6 cos 6 dé ai (e&t/m+u) ae (5) 
where a is written for kl/(1+iw7) and 
q dx 
Teun 
1 Ogee (8) 


Now from J, the value of n may be found, since according to the equation 
of continuity J,=nev,. Hence from (5) and (4) 


_ 3Nefeés/m+u 


eo “aki 1+fl(iwr) > . . . . . . . (7) 
; ie cos 6—f/ —1 
and Vp =(e6 7/m+u) ee : o> lg, incae CL 


For the purpose of calculating dissipative effects we are interested 
not so much in v9’ as in the difference, 4vg’, between vp’ and the value 
U eq given by (2): 
At, = (e&'r/m+u) oe (Oe ee 
(l+-twr—tkl cos 0)[1+-f/(tw7)] 
Thus the departure of the Fermi surface from local equilibrium is deter- 
mined apart from the value of &. Provided the frequency is much lower 
than the plasma frequency w p(=[47NePc?/(me)]*/?) the space charge needed 


—ucosé. . (9) 
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to develop the required value of & is so small that no serious error results 
from assuming no space charge, i.e., putting J,= Neu (— Neu is the lattice 
current density). Then & is found from (7), and (9) takes the form 


Ain th. a? cos 6 a 9 
Ore 360 2iaco8 6) 314—COs O>. See LO) 
In general, however, if we do not assume J -—Neu, but instead write the 


total current density as J,—Neu and relate it to & by means of the 
equation (in electromagnetic units) 


Uwe 
J,—Neu=— Trea? 
then Av,’ takes the more complicated form 
Ave! uw (tot —w*/w,")[(1—f) cos 6—ifja]_ Ee sh 
{ (37fv,/av9) —w?/w ,?[ 1+ (fiw) }(1+iw7—ckhl cos 6) 


(11) 
The simpler form (10) is recovered by neglecting the terms in (11) which 
contain w?/w,”, and which are, in all practical applications to metals, 
many orders of magnitude smaller than the other terms. It is of interest 
to note, however, that the plasma terms appear in such a form as seems to 
allow no conditions under which they could contribute any marked 
effects, even if, as is perhaps possible in semiconductors, w becomes 
comparable with w,. The discussion which follows will be based on (10). 
Because of the departure from equilibrium represented by Av,’ in 
(10), there will be a continuous irreversible conversion of electronic energy 
into heat by collisions with the lattice. If the mean energy excess over 
that corresponding to local equilibrium is 4H per unit volume, the rate 
of production of heat will be 24#/7 per unit volume, the factor 2 arising 
from the definition of 7 as the relaxation time for the velocity, while the 
energy excess is proportional to | 4v,'|?. In fact if w represents the peak 
amplitude of the particle velocity in the wave, at any part of the Fermi 
surface the displacement oscillates between -+-| 4v,’|. If we consider an 
element dS of the Fermi surface, at the moment when the displacement is 
+ |Av,’'|. there are 3NdS|Av,'|/(47v) electrons displaced above the 
equilibrium level, each with an average excess energy 4mvp| Av |, so 
that the excess energy of this group of electrons oscillates between zero 
and 3NmdS|Av,'[?/(87v52), with a mean value 3NmdS|Av,'[?/(1670,"). 
Hence for the whole Fermi surface 


3Nm ([ 3Nm (* 
ce Av,’ |? dS= Av,’ |? sin 6 dé, 
oe eae Yo |? de 8 |, 0 | 
and the mean rate of heat production per unit volume 
3Nm (7 pies 
Va. ne ein Ooi) Pa ee wake emcee LZ) 


In general a and f are complex, and the evaluation of |42,‘ |? from (10) 
is very tedious. But in the interesting range of frequencies for which a 
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is not too much greater than unity, so that wr is still much less than unity, 
it is sufficiently accurate to take a as equal to kl, when f becomes real and 
equal to 1—tan-a/a. Then after some reduction (12) takes the form 


Nm 42 Jt —1lj, 
oz Nmu? Lk [2 baad = (13) 
2Q7 3 ki—tan-1kl 


A more searching analysis of | 4v,’ |? shows in fact that (13) is a very good 
approximation even when w7 is not negligible, and may be taken to apply 
at all frequencies. 

The attenuation constant of the wave may now be written down 
immediately, since 7, the attenuation of energy per unit length, is given 
by 2 Q/(v,pu?), p being the density of the metal. 


Nm (; k2]? tan—1 kl ) 


uae “ 3ki—tan-! kl 


= —— (14) 
pv .T 

At low frequencies, for which the mean free path is much smaller than the 

wavelength, that is kl<1, 


4Nmv,2w T 
l5pv < 


9 23 iA 1 
= — =~ (kl)?+ —~(kl)*—...). et 5 
x (1— es Fee...) (15) 
As the frequency is raised so that kl exceeds unity, the law of attenuation 
approaches the limiting form 
a7Nmvyw 


16 
S03 (16) 


AS 
The leading term of (15) is the same as is obtained by the simple treat- 
ments of Mason and Morse, which neglect free path effects. Morse’s 
analysis suggests that the attenuation continues to rise as w? until such a 
frequency as makes wr comparable with unity, when it flattens off to the 
constant value of 4Nmv,?/(15pv,37). The present analysis shows however 
that at the much lower frequency at which k/=1 the attenuation becomes 
proportional to w and independent of 7. It is perhaps worth remarking 
that this behaviour is reminiscent of the anomalous skin effect (Pippard 
1947, Reuter and Sondheimer 1948) which also becomes independent of 7 
when the mean free path is very long. Moreover, just as in the latter 
phenomenon most of the current is carried by a small group of electrons, 
so here, in ultrasonic propagation, when / is very long the deformation 
of the Fermi surface is mainly confined to the small region occupied by 
those electrons whose velocity component along the direction of wave 
propagation is equal to the velocity of sound. 


§3. TRANSVERSE WAVES 


In transverse waves there are no density changes and hence no electric 
fields resulting from space charge. On the other hand it is possible that 
the lattice and electronic currents may not cancel one another, in which 
case magnetic fields will be generated, and from these electric fields by 
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induction. Ifthe transverse wave is propagated in the z-direction, with a 
particle velocity w in the y-direction, the magnetic field will lie along the 
x-direction and the electric field along the y-direction. If the electronic 
current (in the y-direction) is J,, then according to Maxwell’s equations 


6 (0? =—h°E=4riw) =47iw(J,—Neu). . . . . (17) 
The analysis of the deformation of the Fermi surface proceeds exactly 


as before, leading to the result that for an electron moving at an angle 0 
to the z-direction and at an azimuthal angle ¢, 


= e6 r/m—+u 
~— |+iwr—ikl cos 6 


y 


VU sin @ cos ¢. Pk «8 


The value of vy’ corresponding to local equilibrium, that is, a uniform 
translation with velocity wu, is wu sin@ cosd. The electronic current density 
is found from (18) by integration over the Fermi sphere, 


J .=Neg(e&r/m+u), where g=$(1—f/a?—f), eee (L9) 
and f and a have the same significance as before. Hence from (17), 


(18) and (19) (ky2+k?)/ (gk? k2)u sin 6 cos d 


oe I-Fiar—thl cos 0 ee cw) 
ky? +k? 
‘= peel 
and Ag, = mR mss ae TOC8 0) 1} wsiné cos ¢, (21) 


in which k,?=47iwNe?t/m=4riwo, o being the conductivity. 

The detailed behaviour of 4v,’ is seen to depend on the relative magni- 
tudes of k and ky, that is of the sound wavelength and the classical skin 
depth of electromagnetic waves of frequency w. We shall not discuss all 
the cases that can arise, but confine our attention to that of most interest, 
in which / is long enough for kl to reach and exceed unity while 4? is still 
much less than /,2. Then so long as k? is also much less than gk,” we 
may write for Av,’ 


oh l/g aan 
Av, = (mes = 1) usin @ cos ¢, 


and, since g may be taken as real until such high frequencies that o7~1, 


(1/g—1)?+k/? cos? 6 
1-- kl? cos? 0 


[Avy |?= u? sin? @ cos? ¢. 2) 


The attenuation constant is now calculated exactly as before, yielding the 


result IN (23) 


? 


| 
where g=3/(2k2?)[(k2l?+ 1)/(Al) tan~"(kl)—1] and v, is the velocity of 
transverse waves. When kil<1, « may be expanded as a series, 


Ninvg?o*r Dv hus pacer oO ) one (24 
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As with longitudinal waves, the leading term in « is just what would be 
obtained by the simple relaxation theory, and comparison of (15) and (24) 
shows that the free path effect gives rise to very similar initial departures 
from this theory. 

As kl attains values greater than unity, while w7 is still very small, 
g tends to the value 37/(4kl), and « to the form 


AN mvyw 
—_ 


[SS 25 
a 37pv, ( ) 


which is similar in form to the corresponding result (16) for longitudinal 
waves. However, as k/ increases further it begins to be invalid to neglect 
k? in comparison with gk ? in (21). The critical frequency at which 
| gkp?| equals k? is given by (37?v,'30/1)!/2, which is independent of the value 
of r and equal to about 10° c/s for most metals. At about this frequency 
x begins to fall and eventually becomes independent of frequency. This 
last stage occurs either when k,)? becomes much smaller than k?, so that 
the skin depth exceeds the wavelength, or when wz becomes much greater 
than unity. Under either circumstance the electric field falls to a low 
value and the electrons are not appreciably affected by the passage of the 
wave, except for a small group, having v, cos 0=v,’, which move with the 
wave and acquire a drift velocity w in the y-direction. This group, 
however, in contradistinction to its importance at lower frequencies and 
in longitudinal waves, now plays a quite negligible role. The limiting form 
of « at the highest frequencies is given by the expression 

ee 

put 
being independent of frequency and inversely proportional to the relaxa- 
tion time. 
§4,. APPLICATION TO REAL METALS 


There is no need to repeat here the comparison which Mason and 
Morse have made between their theories and Bémmel’s experimental 
results, since the experiments refer to a regime in which ki is somewhat 
less than unity, and there is no significant difference between the results 
of the simpler theories and those obtained here. The satisfactory nature 
of the agreement obtained shows that the free-electron model gives a 
good semi-quantitative account of the real situation. If the electron gas 
is not free, but has band structure, the analysis given here corresponds, as 
Kittel’s (1955) analysis shows, to transfer of electrons from one part to 
another of the same Brillouin zone. It should be remembered, however, 
that if the electrons occupy more than one zone, the energy levels in two 
zones may be differently affected by uniform compression, and in the 
passage of longitudinal waves there may occur a periodic transfer of 
electrons between zones, with a relaxation time perhaps quite different 
from that used here. This effect would not occur with transverse waves, 
so that studies of the attenuation of both types of wave in the same 
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specimen would enable it to be detected. From a formal point of view 
this effect may be described as a second viscosity of the electron gas, 
leading to an imaginary term in the bulk modulus. The possibility of 
second viscosity is considered by Mason, who concludes from Bémmel’s 
measurements that it is much smaller than the first viscosity, the effect 
which we have analysed at length here. 

It is of interest to consider the attenuation to be expected at very high 
frequencies, of the order of 1011 c/s or more, such as occur in the thermal 
excitations of the ionic lattice. Transverse waves will not be highly 
attenuated, since the conditions will be such that (26) applies, but the 
attenuation of longitudinal waves will be considerable, as described by 
(16). There seems to be no reason why 1/« should not be regarded as 
defining a mean free path for phonons, provided the phonon wavelength 
is large enough for the model used to be valid. Putting in reasonable 
values of the parameters, and taking w equal to 5x 10! (corresponding 
to a temperature hw/k equal to 3°K), we find that the effective free path 
1/x is about 2x 10-%cm. This is of the same order of magnitude as the 
free path calculated by the quantum theory of electron-phonon inter- 
actions, and indeed what we have analysed here is the classical analogue 
of one of the principal scattering mechanisms, the interaction between 
electrons and the periodic potential associated with a longitudinal phonon. 

The magnitude of this potential may be simply derived from the electric 
field, since .kd=€, and from (7) we have, putting J,= Neu, 


mu {akl in : 


Now the equilibrium value of the Fermi velocity, and hence of the kinetic 
energy at the Fermi level, varies periodically with the wave, and so does 
the potential energy ef. The total energy varies with an amplitude 
ebd+mv,2u/(3v,), which from (27) may be written in the form 


U0) tae. 1—I1/f 3 is 
aa —— +57. - - + = (28 

AB toy 30, stor |e T Re ee 
It may be seen from this expression that at all frequencies 4F,,, is very 
much smaller than the kinetic energy variation mv,?u/(3v,) ; the field is 
always able to adjust itself to maintain the Fermi level practically constant. 


§5. SUPERCONDUCTORS 


Boémmel’s and Mackinnon’s observations that the attenuation falls 
rapidly as the temperature is lowered through the superconducting 
transition raises a number of interesting points. It might at first be 
regarded as a matter for some surprise that the behaviour is the same 
for both longitudinal and transverse waves, for in the former the accelera- 
tion of normal electrons is accomplished mainly by electric fields, while 
in the latter collisions are nearly sufficient to produce the required motion. 

‘Tf one is accustomed to think of a superconductor as a metal which cannot 
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support an electric field it is not obvious that the transition to super- 
conductivity should affect both wave types similarly. This difficulty is, 
however, readily resolved. It is true enough that in a uniform super- 
conductor an electric field gives rise to a steadily increasing electric current, 
but in a non-uniform superconductor it is not simply the electric field 
which is the accelerating force, but the gradient of the chemical potential 
yu, of which the electrical potential is only one part. This may be seen 
by considering a thermoelectric circuit consisting of two different super- 
conductors ; the known absence of a thermoelectric potential difference 
leads immediately to the consequence that is constant over all parts of a 
superconductor in which current is not changing, even if it is not at a 
uniform temperature. Now, as has just been shown in §4, the effect 
of the electric field in a longitudinal wave is to maintain the Fermi level, 
and hence », practically constant, so that there will be no tendency for 
any ‘superconducting electrons’ present to move in such a way as to 
eliminate the field. There is thus no significant difference between 
longitudinal and transverse waves in this respect, and the following 
discussion may conveniently be confined to the latter. 

If the two-fluid model of a superconductor be accepted it is under- 
standable that on passing the transition temperature the attenuation due 
to normal electrons should drop, since these become less in number as the 
superconducting electrons increase. Whether or not the latter give rise 
to additional attenuation is a question that cannot be resolved without a 
more detailed model, but if they do not, the most elementary view suggests 
that the observed drop in attenuation merely reflects the decrease in the 
number of normal electrons. Yet the rate of drop is very rapid, far more 
rapid than would be predicted by any existing model of the super- 
conducting state (e.g., Gorter and Casimir 1934), from which a decrease 
roughly as 7 might be expected. or this reason Morse (1955) suggested 
that the normal electrons are more tightly coupled to the lattice in the 
superconducting state, that is, there is a decrease in free path as the 
temperature is lowered past the transition. This explanation, however, 
cannot be accepted without question, since there is other evidence which 
appears to contradict it. For example, the thermal conductivity in tin 
does not begin to fall rapidly when the superconducting state is reached, 
as it would if there were a decrease in the free path ; indeed, in order to 
explain the temperature variation of thermal conductivity Heisenberg 
(1948) was led to suppose that the free path of normal electrons was 
actually greater in the superconducting than in the normal state. It 
seems necessary therefore to discover a mechanism which decreases the 
free path concerned in the attenuation of ultrasonics without decreasing 
that concerned in the conduction of heat. One can only speculate about 
what this mechanism may be, but at least it does not seem an impossible 
task if the additional coupling between the normal electrons and the ionic 
lattice is effected through the intermediary of the superconducting 
electrons. 
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Let us consider, for example, a probably unrealistic model of the 
superconducting state, in which the superconducting electrons, amounting 
to a fraction (l1—x) of the total number, occupy a sphere of about the 
same radius as the Fermi sphere, while the remaining fraction, x, of normal 
electrons occupy a similar sphere, so that the total density of states remains 
constant. Let us suppose the model has the following properties : 

(a) The superconducting assembly, unlike the normal assembly, is 
incapable of thermal excitation. This is consistent with the absence of 
thermoelectric effects in superconductors. 

(6) The superconducting assembly is tightly coupled with the lattice, 
so that when a transverse wave passes the sphere remains spherical, but 
acquires the local lattice velocity. 

(c) The normal assembly interacts with the lattice in the same way, 
and with the same free path, as do the electrons in the normal state of the 
metal. 

(d) The two assemblies interact strongly with one another. It is not 
unreasonable to postulate this since it is known (e.g., Pippard 1951) that 
the width of the transition region between normal and superconducting 
phases may be much smaller than the free path. 

First let us neglect the interaction (d). Then at not too high ultrasonic 
frequencies both normal and superconducting assemblies will drift with 
the transverse displacements so that only a very small net current is set 
up by the wave. But while the superconducting assembly remains 
spherical, the normal assembly will be deformed in a manner similar to that 
expressed by eqn. (21). This will lead to a similar dissipation of energy 
as heat in the same way as in the normal metal, but reduced by a factor x, 
There will, however, be a second consequence—at different parts of the 
Fermi surface there will be a different energy gap between the ° Fermi 
levels ’ of the two assemblies. Since the magnitude of this gap is deter- 
mined by the processes leading to the superconducting condensation and 
is not arbitrary, it seems likely that the interaction (d) will lead to a 
sharing of momentum between the two assemblies so that to a good 
approximation the gap stays constant at all points, and each surface has 
the same deformation from spherical form. In order that momentum 
may be locally conserved at all points of the Fermi surface, the effect of 
the interaction will be to reduce the deformation of the normal assembly 
by a factor a. Since therefore the attenuation of the sound wave is 
proportional to the number of normal electrons and to the square of the 
deformation we may expect « to vary as x*; according to the Gorter— 
Casimir theory x varies as 7", so that « may vary as 7’. This rapid 
variation is roughly what is needed to account for the observations, 
although the data published are not sufficiently detailed for any precise 
exponent to be fitted to them. 

It should be noted that the interaction between the normal and super- 
conducting assemblies need not alter in any way the thermal excitation of 
the former, since the latter is assumed incapable of excitation. Thus 
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this process does not lead to a diminution of the free path which is 
applicable to problems of thermal conduction. It need hardly be remarked 
that the model developed above is not intended as a serious contribution 
towards a theory of superconductivity. It is introduced simply to point 
out the sort of problem raised by the experimental observations and to 
suggest the lines along which a solution might be possible. There may be 
other solutions along quite other lines, but it seems that the behaviour of 
ultrasonic attenuation in superconductors is not a trivial matter by any 
means, and that these new observations may contribute valuable evidence 
leading to a better understanding of superconductivity in general. 
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CXXIIT. Thermodynamics of a Sheared Superconductor 
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RECENT measurements by Olsen (1955) have shown that the rigidity 
modulus n differs in the normal and superconducting states of tin by an 
amount which varies with temperature approximately as (l= StS); 
without a discontinuity at the transition temperature, 7',. This work was 
carried out by torsion measurements on polycrystalline samples, while 
Landauer (1955), working with single crystals of tin, found evidence for a 
discontinuity at 7’,; otherwise the temperature variation was similar to 
that found by Olsen. In this note we give a simple thermodynamic 
treatment of the problem which provides an explanation of Olsen’s result 
and reveals a difficulty in reconciling it with that of Landauer. 

Let us consider unit volume of a single crystal specimen in the form of a 
rectangular or nearly rectangular parallelepiped, and take the edges of the 
specimen to be the xyz-axes. Let it be subjected to a uniform shear stress 
+ applied along those faces which contain the z-axis, and let the corres- 
ponding shear strain be 6. Then it is readily shown by an extension of the 
usual thermodynamic treatment of the phase transition in superconductors 
that the critical magnetic field H, is determined by the equation 


PO ae ee es a8 


in which G,, and G, are the values in zero field, for the normal and super- 
conducting states respectively, of the potential U—7’S—2é@ which has 
the property that (dG/02’),;——6. Hence, from (1), 


Hef oi . 
Ad=0,—0,, = is (=), . 5 e > : ' 0 (2) 


The significance of this result is readily apparent. For a cubic crystal, ora 
crystal of lower symmetry cut so that the a- and y-axes are axes of 
crystal symmetry, positive and negative values of 2 are equivalent, so 
that (@H,/02),—0 and correspondingly #,=—0,. But in general (dH,/0~), 
need not vanish, and the non-vanishing value of 0,—@, corresponds to a 
spontaneous shearing of the crystal at the phase transition, even in the 
absence of a shear stress. There is no need to assume a change in crystal 
symmetry ; for example a spontaneous shearing could occur ina tetragonal 
crystal such as tin if the axial ratio were different in the superconducting 
and normal states. It is necessary in this case, however, that the tetrad 
EO hs 5 Na ere ED a i eer 
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axis shall not lie along any axis of the specimen, for then the specimen will 
only change the ratio of its sides without deformation of its angles. 

Differentiating (2) with respect to 2 we obtain an expression for the 
change in isothermal rigidity modulus n(=(02/06);) when the phase 
transition is initiated by a magnetic field, 


a Ly ene | (ce Oh 
An=n,—N, = —n?A (-) =— ie {(S5), +H, ( Ae ),} 5 : (3) 


It is here assumed, as is borne out in practice, that 4dn<n. Under such 
conditions as those considered above, when 46=0 for Y=0, Mn takes the 
form —(n2H,/47)(02H,/02?)7, and clearly vanishes at the transition 
temperature since there H,=0. On the basis of simple assumptions 
Olsen’s result may now be derived from this expression. These assump- 
tions are that H, varies parabolically with temperature and that the 
electronic specific heat is not appreciably altered by shearing. Then H, 
takes the form H,(1—T7?/7'.2) and H,/T7', is independent of X. It follows 
from (3) that, since we are assuming dH,/d2=0. 
nH, d*H, 


An Sg UTD es 


in agreement with Olsen’s results. 

If, on the other hand, we do not assume 4é@ to vanish when 2’=0, the 
first term in (3) also contributes to dn. In particular as the transition 
temperature is approached dn tends, not to zero, but to the value 


: A6\2 
(An)p = lim < —4rn? | — J: «ge Sa ie eee 
c H { 
T>Te c J 
Consequently if 4@ rises, like H,, linearly from zero as the temperature is 
reduced below 7'., there will be a discontinuity in 4n at the transition. 
It might seem, at first sight, as if the explanation of the discrepancy 
between Landauer’s and Olsen’s measurements might be resolved by the 
use of these results, but in fact no such easy solution seems possible. 
For according to (5) (dn), is essentially negative, that is, the rigidity 
cannot increase when a metal becomes superconducting ; if then there are 
ya Ty ene ee ety ; Baas : 3 
any crystal orientations for which (4 n)p, does not vanish these will show 
up when a polycrystalline sample is used, as in Olsen’s experiments. 
The fact that Olsen observed no discontinuity in rigidity implies therefore 
that there are no orientations for which the discontinuity is very great, 
in other words the axial ratio in tin does not suffer significant change 
bet ween the superconducting and normal states at any temperature. It 
is Important to note that one cannot argue from the effective isotropy of a 
polycrystalline sample that 4é in (2) vanishes and in consequence (4), 
ra re ; ‘ a4 1a 7 rahi ‘ r a > 2 
must also vanish; for it is ins alid to apply the arguments given above 
to any but a monocrystalline specimen. This may be made clear by 
considering the composite specimen shown in the figure, consisting of two 


crystallites with their tetrad axes symmetrically disposed. Each 
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crystallite may suffer a spontaneous deformation 4@ when the phase 
transition occurs in a field H,, so that according to (5) there may be a 
discontinuity in at the transition temperature. But the application of 
the same argument to the specimen as a whole would lead to the con- 
clusion that 46 was necessarily zero, and hence also (An)p,. The flaw in 
the latter argument lies in the assumption of a homogeneous sample. The 
composite specimen becomes thermodynamically inhomogeneous as soon 
as a stress is applied, since the stress tends to increase the value of H, 
in one crystallite and decrease it in the other; thus no critical field is 
definable in the presence of a stress, and the thermodynamical treatment 
breaks down. The polycrystalline specimen is merely a more elaborate 
example of the same type of system, and one can see that the discontinuity 
in nm must be calculated as the mean over all crystallites, and not found 
by applying a spurious thermodynamical argument to the specimen as a 
whole. 


Composite specimen containing two crystallites, with tetrad axes along the 
directions indicated by the lozenges. The dotted lines show the effect of 
a spontaneous shearing at the phase transition. 


For the same reason care must be exercized in considering a monocrystal- 
line wire subjected to a torsional couple J" and suffering thereby a twist ¢. 
If » and @ in the above argument were replaced by I and ¢ it would 
immediately follow, if the argument were sound, that no discontinuity in 
the torsion constant could occur at the transition temperature, since at 
no temperature, when J’=0, would the wire spontaneously twist when the 
phase transition was initiated by a magnetic field. The fallacy in such 
an argument is again that a twisted monocrystal is in general inhomo- 
geneous, the critical field being raised on one side and lowered on the other 
exactly as in the example shown in the figure. By consideration of each 
element of the wire separately, however, it follows that (An), does not 
necessarily vanish and that therefore a discontinuity in torsion constant 
may occur at the transition temperature, provided the crystal is not so 
oriented as to preclude a discontinuity in » for reasons of symmetry. 
Now in fact Landauer’s experiments were performed with a erystal of tin 
having the tetrad axis very nearly perpendicular to the axis of the wire, 
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and this is one of the special orientations for which (4m), should vanish. 
This is because the shear stresses produced by twisting lie in planes which 
contain the tetrad axis. In order that a discontinuity in 7 shall be 
observed it is necessary to choose a crystal orientation such that the wire 
may spontaneously shear at the phase transition and not merely change its 
cross section, as with the orientation used by Landauer. It is, of course, 
possible that Landauer’s observation of a discontinuity in” at 7’, is due 
to the fact that the tetrad axis lay a few degrees off the normal to the 
wire axis. But if this is so it is to be expected that there are other 
orientations for which (4dn)p, is even higher than that observed by 
Landauer, and it is then remarkable that Olsen found no discontinuity with 
a polycrystalline specimen. Itseems more likely that Olsen’s result is to be 
believed, that (4), is small for all orientations of tin, and that Landauer’s 
result arises from some instrumental error associated with the phase 
transition. 
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CXXIV. Formation of Etch Pits on Slip Lines in Aluminium 


By G. Wyon and J. M. Marcutn 
Centre de Recherches Metallurgiques de l’Ecole des Mines de Paris* 


[Revised MS. received May 19, 1955] 


ReEcENT work suggests that it is possible to find a close relation between 
the location of etch pits and the points at which dislocations emerge through 
the surface of slightly deformed metallic or ionic crystals. Examples are 
_ experiments on 35% Zn brass (Jacquet 1954), on ionic crystals (Amelinckx 
1953) and on pure single crystals of Zn (Abdou 1954). 

In the course of the present work, it has become clear that in the 
particular case of Al, great care must be taken in the interpretation of etch 
pits formed by dipping specimens into acid. In previous publications 
(Wyon and Marchin 1954, Wyon and Lacombe 1954) it has been shown 
that for Al the presence of impurities on dislocations appears to be the 
cause of the localization of etch pits on grain boundaries and sub- 
boundaries and even of the etch pits distributed at random on crystals’ 
surfaces. On the other hand the experiments of Amelinckx on etching of 
slip bands in Al seem to indicate that dislocations alone are responsible 
for the localization of etch pits on slip lines. 

In order to verify this very important point, these experiments have been 
repeated and the following method has been used: refined aluminium 
samples (0-004°, Fe and 0-005 $i) in sheets about 1 mm thick were recry- 
stallized with large grains by annealing at 620°c for eight days. The 
samples were electro-polished in a cooled Jacquet bath (666 cm? of 
(CH,CO),0 and 333 cm? ClO,H) and subsequently slightly deformed in 
tension to produce some slip lines. In this condition, the sheets were 
immediately etched, first in fluorated aqua regia (Wyon 1953), then in a 
reagent derived from the preceding one and including a wetting agent 
which emphasizes the grain and sub-grain boundaries, and which also 
produces micro-etch pits (Wyon 1953). This etching was preceded by an 
electrolytical re-polish after the deformation in order to determine the 
part played by the dislocations locked on the slip planes inside the crystals. 
On the basis of these experiments, the following observations could be 


made. 


ErcHING AFTER DEFORMATION WITHOUT REPOLISHING 


We have found that generally the slip lines apparent on the surface 
are only etched if the crystal has an orientation such that a (100) axis 1s 
approximately parallel to the surface. After etching, micro-etch pits 
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appear with indefinite shapes as observed by Amelinckx (fig. 1 (figs. 1-3, 
Plate 23)). In certain very rare cases, there occur for grains with any 
orientation, true micro-etch pits along the slip bands (fig. 2). The 
distribution of etch pits along the slip bands can vary widely. 

We have observed : 

(a) compact local alignments of etch pits, as seen by Amelinckx (fig. 1), 

(b) alignments similar to (a) but continuous along some slip lines (fig. 3), 

(c) clearly marked micro-etch pits occur well separated from each other, 

(d) micro-etch pits arranged in a broader and more diffuse cluster, 
overlapping in the region of crystal around the slip line. 

Generally several slip lines remain unetched along their length, and the 
following remarks can be made. 

1. When several slip planes are acting in the same grain : 

(a) one system out of two, or two systems out of three may be etched 
(fig. 3); all active systems may be attached but not all the slip lines ofa 
single system may show the same phenomenon, 

(b) etch pits may or may not be observed at the ends of cross-slip lines. 

2. When a slip line stops in front of a crystalline boundary, the micro- 
etch pits, when they occur, are spaced further and further apart while at 
the same time the line disappears (fig. 5 (figs. 4-6, Plate 24)). There is never 
a concentration of etch pits at the end of such slip lines. 

3. At the point at which a slip line meets a grain boundary, the micro- 
etch pits on those lines which show them suddenly stop at the boundary, 
even if the lines cross the boundary (fig. 6). But not all the slip lines have 
micro-etch pits up to the grain boundary (fig. 6). 


ETCHING ON SURFACE RE-POLISHED AFTER DEFORMATION 


Figure 7 (figs. 7, 8, Plate 25) shows the appearance of an electrolytically 
polished specimen slightly deformed in orderto develop some slip lines, then 
repolished in a Jacquet bath in order to smooth out any slips on the surface. 
The sample has also been etched and no alignment of etch pits is observed. 
A second deformation is then carried out, in order to show the direction of 
the active sliplines ; these are the last formedslip linesseen on fig. 7. Starting 
with the idea that dislocations by themselves cannot (in the case of Al) be 
revealed individually by an acid, we attempted to reveal the slip lines 
obliterated by polishing. For this purpose an ageing treatment was 
applied after deformation to enable impurities to diffuse to dislocations. 
The samples were then polished and etched again. In the case of refined 
Al, it was never found possible to reveal the slip markings by this method. 
In fact the impurity content is very small. The impurities (generally 
iron) have a low mobility at the low ageing temperatures which one is 
obliged to use to prevent other effects, e.g. * polygonization ’, from 
occurring (15 days at 100° for instance). The experiment gave a positive 
result, however, with an alloy of Mn (0-56%): fig. 8 shows the slip 
markings observed on a sample repolished and etched after the deforma- 
tion and subsequent ageing treatment. 
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CONCLUSIONS 


From the experiment described above it appears that dislocations alone 
are not sufficient to explain the localization of some etch-pits on the slip- 
lines of aluminium. If it were so, at least sometimes a concentration of 
etch-pits should be observed at the ends of the slip lines stopped at grain 
boundaries (Jacquet 1954). The etch pits should continue to be localized 
on the slip planes in the areas under the surface. On the other hand, the 
etch figures seem to occur only on the slip lines in crystals having a fairly 
particular orientation and even then some slip systems do not show these 
etch-pits. Nevertheless, the influence of some factors should not be 
disregarded, for instance the rupture of the surface oxide film at the slip 
lines or the * microgeometric ’ factor caused by the undulations of the 
surface at these slip lines, can modify profoundly the local reactivity of the 
metal. 

In the light of some previous results (Wyon and Lacombe 1954) and 
others to be published soon, it may be concluded that dislocations alone 
cannot modify the reactivity of Al with regard to the acid reagents 
sufficiently to give rise to etch-pits at the points of emergence of the 
dislocations. But in Al alloys with small amounts of the alloying element 
such as the above-mentioned Al 0-56°% Mn, the Al 3°% Mg alloy studied by 
Crussard and Aubertin (1951) and alloys of Zn (6—12°/%) (Berghezan 1952) 
it is possible to demonstrate without any doubt the alignment of etch-pits 
on slip-band markings, even in the areas of crystals lying under the surface. 
But this result is only obtainedifthesamples are aged after the deformation ; 
this ageing treatment promoting without doubt the combination of 
dislocations with impurities which can be detected by etching. From this 
point of view, Jacquet’s work on 65-35 brass also shows that impurities 
play their part in corrosion phenomena on the slip lines ; the reagent used 
by Jacquet is unable to reveal dislocations in the case of fairly pure copper. 
For this purpose, it is necessary that the Zn-content is sufficiently high 
and close to the content at which the beta-brass phase appears. 

It is known that in the copper lattice, the dislocations are easily split 
up into half-dislocations giving rise to ‘stacking faults’ (Heidenreich 
and Shockley 1948). Nevertheless these stacking faults alone do not 
cause dislocations to be revealed by etching with hyposulphite reagent. 
On the other hand, in the case of 65-35 brass, the presence of a high 
proportion of Zn atoms in the face-centred cubic matrix enables isolated 
dislocations to be revealed by etching. 

All these experimental results show that the problem of revealing the 
dislocations by micrography in metallic lattices is rarely simple and that it 
is necessary to perfect etching reagents which are more and more sensitive, 
reproducible and selective. 
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SUMMARY 


The electrical discharge theory of the variation in excitation observed 
in the atmospheres of the long-period variable stars offers an explanation 
for the combination spectra exhibited by many of these stars, which is 
in accord with several of the most outstanding changes in their spectra 
and magnitude. 


§ 1. INTRODUCTION 


Ty an earlier note (Bruce 1955 a) it was shown that the hypothesis of 
electrification due to freezing in low temperature stellar atmospheres 
would account for many of the phenomena observed in the spectra of the 
long-period variable stars. Subsequently it was suggested (Bruce 1955 b) 
that the simpler and better established phenomenon of electrification due 
to asymmetrical impacts will probably account satisfactorily for all the 
electrical phenomena observed in dust storms, volcanoes and terrestrial 
thunderstorms, as well as those already mentioned in the atmospheres of 
the long-period variable stars. 

It is the purpose of the present report to suggest that the electrical 
discharge theory may account for the puzzling combination spectra 
observed in some of these stars. The states of excitation indicated by 
these spectra are so diverse that they have been, and often still are, 
believed always to have their origin in a binary star and an associated 
nebula. However, a consideration of the nature of the long electrical 
discharges, already postulated to account for the general characteristics 
of long-period variable stars, suggests that they may sometimes account 
for both the nebular spectrum, and for the spectrum which has hitherto 
been attributed to a companion star. 


§ 2. GENERATION OF STELLAR ATMOSPHERIC ELECTRIC FIELDS 


It is well known that asymmetrical interactions between pieces of 
the same material can lead to separation of electric charge (Shaw 1926) ; 
the two pieces of material carry away the complementary charges of 
opposite polarity. Impacts between particles of different sizes and shapes 
eee eee Pelee ee ee See See 
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are an example of such asymmetrical interaction, and it is well known that 
electrification does result from such impacts (Shaw 1928). It is thus 
to be expected that the larger particles will on an average carry one charge 
and the smaller particles will carry the complementary charge of opposite 
polarity. When such impacts take place in a gravitational field, then 
extensive electric fields will be built up. In sand and dust storms and in 
the debris which is ejected from volcanoes, electric fields generated in 
this way can reach values which lead to electrical breakdown of air at 
atmospheric pressure, and to long electrical discharges. 

Though the electrification in terrestrial thunderstorms is often attributed 
to charge separation occurring in the freezing process itself, it would appear 
that the same process of electrification by asymmetrical impacts will 
account satisfactorily for the separation in that case also, as was indeed 
in effect suggested by Simpson and Scrase (1937), when they first showed 
that ice particles, rather than water drops, must be involved in the 
separation process. Shaw (1928) found that the electrification of cold 
dry ice particles seemed to be of the same order as that observed with sand 
particles in his experiments with particles blown by air blasts. 

The same conclusion applies regarding the electrification processes 
in stellar atmospheres which have been invoked to account for the 
phenomena observed in the atmospheres of long-period variable stars 
(Bruce 1955 a), though the actual freezing out of solid particles from the 
atmosphere remains of fundamental importance. The particles involved 
will be the oxides, carbides, hybrides of the metals, as well as carbon and 
other materials, and ultimately the metals themselves, which freeze at 
temperatures from about 4000°K downwards. Electric fields will only be 
built up, however, when the electrical conductivity of the atmosphere 
with electric fields up to that required for breakdown, has fallen below 
the rate of charge separation due to impacts. The breakdown field in an 
atmosphere which is largely composed of hydrogen will be given approxi- 
mately by (PX,/7'), in which P is the pressure in atmospheres, 7’ the 
absolute temperature, and X, is the potential gradient required for break- 
down in a uniform field at atmospheric pressure, i.e. it will be of the order 
of (104P)/7' volts per em, where P is in atmospheres. 

The rate of freezing and of charge generation will depend on the 
frequency of collisions between particles and will therefore probably 
increase with the square of the gas pressure, as well as with the gravita- 
tional force. Both will increase towards lower regions of the stellar 
atmosphere. The rate will also increase with the velocities of the 
particles. The breakdown field itself is proportional to the gas pressure, 
so that it is to be expected that the time required for the electric field 
to reach the breakdown value will vary inversely as the gravitational force 
and gas pressure and inversely too as the velocity of the particles or as 
some power of that velocity. 

In terrestrial thunderclouds the time required for freezing and field 
generation is of the order of minutes or tens of minutes while the electrical 
status quo is reconstituted after the occurrence of a lightning flash in a 
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few seconds. The values of the above three quantities are not known 
precisely in the atmospheres of the long-period variable stars. If we 
may assume that the pressure is of the order of 10-6 atmosphere, that 
gravity is of the order of 10-4 the value in the earth’s atmosphere, and 
that the velocities involved are of the order of about 103 times those 
observed in thunderclouds, then the time required for the process of 
freezing and charge generation in the stellar atmosphere would be of the 
order of 104 to 107 times that required in a terrestrial thunderstorm 
according as the electrification process increases linearly with, or as the 
second power of the velocity. Taking the critical time required in 
a thundercloud as 10? seconds, the time required in the stellar atmosphere 
will be of the order of 10° to 10° seconds, in general agreement with the 
observed values for the periods of those stars of the order of 107 seconds. 


§ 3. CHARACTERISTICS OF LonG ELEcTRicAL DISCHARGES 


In considering the effects which the extensive electrical discharges 
will have in the atmospheres of long-period variable stars it is necessary 
first to describe an important characteristic of long electrical discharges 
to which the writer first called attention in relation to the leader strokc 
of the lightning discharge (Bruce 1944). It had been generally believed 
(Loeb and Meek 1941) that, in order to produce electrical breakdown in a 
long gap, a field of the order of that required for breakdown in a homo- 
geneous gap under the existing conditions of temperature and pressure 
had to exist throughout the whole gap, and remnants of this belief still 
are expressed (Loeb 1954). For air at atmospheric pressure this field is 
of the order of 10 000 volts/em even with water drops present, and fields 
of the order of 3 500 volts/cm have actually been observed during aeroplane 
flights through thunderstorms. 

However, when the field at one point rises to such a value as to give 
rise to corona currents reaching some critical value, a value such that the 
temperature of the gas reaches that required for adequate thermal ioniza- 
tion to carry the discharge, then the discharge becomes an electric 
are (Bruce 1941). That portion of the channel already formed is thus no 
longer a column of corona-like discharge requiring for its maintenance a 
field of several thousand volts per cm, but one requiring a potential 
gradient of the order of only 10 volts/em. The propagation of the leader 
stroke thus represents a partial short circuit of the pre-discharge field, 
and the discharge becomes, so to speak, self-propagating. It also causes 
an increase of the electric field concentration at its tip, over that which 
originally initiated the discharge, and transfers that field concentration into 
regions of the atmosphere in which the electric field was initially quite 
small. Thus if X volts/em is the average value of the predischarge field 
throughout the length, J, in cm, traversed by the advancing tip of the 
leader stroke, and x volts/cm is the potential gradient required to maintain 
the column of thermally ionized are discharge, then the tip will be at a 
potential relative to the surrounding atmosphere, which is greater by the 
amount / (X —2) than the value at the initiation of the discharge. 


- 
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The theory explained how breakdown could occur between thunder- 
clouds and earth with average potential gradients of the order of only 
100 volts per cm. Even with such an average predischarge field of 
100 volts/em, the difference in potential between the discharge tip and its 
surroundings builds up at the rate of about 100 volts/em travelled. By 
the time the leader stroke travels the two kilometres or so to earth this 
potential difference has reached 20 x 16° volts, and potentials of half this 
value at least have been recorded on overhead transmission lines struck 
by lightning. 

The large differences in potential between the channel and its surround- 
ings have an important consequence (Bruce 1944), in that they cause 
corona currents to flow laterally from the leader stroke channel. These 
too are probably of astrophysical importance. 

This thermal ionization, or arc, theory of the leader stroke has been 
verified experimentally by Komelkov (1947) in the case of the laboratory 
spark discharge, while the complete theory of thermally ionized column 
surrounded by lateral corona currents has been found by Pierce (1955) 
to account for the observations made at Cambridge on the electric field 
changes caused by the leader strokes of lightning flashes. 

This property of long electrical discharges of concentrating pre-existing 
differences in potential across relatively small spaces, and of transferring 
these electric field concentrations into regions which were initially subject 
to only small electric stresses, is of considerable astrophysical importance. 
In very extensive atmospheres the result of the field concentration and 
transference will be enhanced when the discharge is propagated into 
regions of lower pressure, and also for the further reason that the thermal 
losses from the very extensive discharge channels become very small. 
In the small laboratory arc discharge, or the leader stroke of the lightning 
discharge, in which the current reaches hundreds, or in intense flashes 
possibly thousands of amperes, the temperature gradients are high 
(King 1954) and the thermal losses correspondingly large. Despite such 
relatively large losses the potential gradient required for its maintenance 
is as we have seen a very small fraction, of the order of 10-3, of that 
required for the initiation of the discharge. In the very large stellar 
discharge channels of fairly uniform temperature and therefore very 
small lateral temperature gradients, the only losses will ultimately be 
those due to radiation, and the potential gradient required for their 
maintenance will be correspondingly small. The process of electric field 
concentration and transference will therefore be effected by the long 
electrical discharges which occur in these stellar atmospheres, even more 
pronouncedly than it is by the leader stroke of the lightning discharge. 


§ 4. DURATION OF THE ELECTRICAL DISCHARGES 
Though the breakdown field will be reached first at some region relatively 


low down in the stellar atmosphere, considerable fields will in the same 
interval have been reached throughout much of that atmosphere. The 
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effect of the important characteristic of these long leader strokes discussed 
in the preceding section will be that, when once the discharge or dis- 
charges are started, they will be propagated throughout practically the 
whole stellar atmosphere. The velocity of propagation of such a discharge 
is not known accurately, but, since it depends on the product of the break- 
down field and the mean free path (Schonland 1938), and since the former 
of these varies directly, and the latter inversely with pressure there is no 
reason to believe that this velocity will differ greatly from the value 
observed at atmospheric pressure, i.e. of the order of 107 cm/sec. Since 
the dimensions involved in these atmospheres are of the order of 10! cm, 
the discharges will last for a period of the order of 10% sec or 100 days, 
in agreement with the times during which emission line spectra are 
observed in the long period variable stars. 


§ 5. COMBINATION SPECTRA 


During the propagation of these long discharges the excitation of the 
matter in and around the discharge channel will not remain uniform. 
As a result of the intensification and projection outwards of the field- 
concentration caused by the discharge, discussed in § 3, and because 
the discharge is being propagated outwards into regions of lower gas 
pressure, the excitation of the gas is increased. The spectrum of the 
discharge will therefore change during the outburst from one of high 
temperature at relatively high pressure, to one of high excitation at low 
pressure, with forbidden lines entering as the very low pressures of the 
outer regions of the stellar atmosphere are reached. It is suggested that 
the theory will thus account for the combination spectra in such stars as 
R Aquarii, Z Andromedae, BF Cygni and AX Persei. In these the 
initial bright line spectrum comprising lines of H, He I, Fe I, Ti II and 
Si IT gives place, after a period of the order of 100 days, to a spectrum of 
higher excitation, containing lines of He II, N IIT, C II, [O IIT], [Ne IT] 
and [Fe III]. The nature of this latter ‘nebular’ spectrum with its 
forbidden lines will be seen to be in accord with the suggestion that it 
originates in regions of very low pressure, far out in the star’s atmosphere, 
towards the completion of electrical neutralization. 

‘The stars showing these spectral combinations and variations thus 
probably differ from ordinary long-period variables in degree rather than 
in kind. They are stars in which the atmosphere outside the photosphere 
is in a more advanced stage of development. The inner atmospheric 
pressure is higher and the atmosphere is radially more extensive. The 
discharges in the inner regions then give rise to the spectrum, which has 
normally been attributed to a small hot ‘companion’ star, while the 
later stages of the discharge yield the ‘ nebular ° spectrum. 

In a study of some of the better known of these combination-spectra 
stars—AX Persei, RW Hydrae, CI Cygni, and Z Andromedae—Merrill 
(1944) concludes that the “‘ maxima and minima of the forbidden lines 
occur about one-fifth of the period after the corresponding phases of the 
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permitted lines’. The periods lie between about 650 and 900 days. The 
observed time intervals are thus in good agreement with the rough estimate 
given above for this interval, namely of the order of 100 days. The time 
intervals between the velocity-curves for the two sets of lines, between 
120 and 200 days are also in agreement with this interval. 

It may be pointed out that, though the belief that in all cases two stars 
and a nebula are required to explain the spectra is widespread, Merrill 
himself in his monograph (Merrill 1940) and papers has been careful to 
emphasize that in many cases, including R Aquarii itself, there is no 
positive evidence for the existence of the ‘companion’ star, as he has 
usually so written it, and that all might in fact come from one M-Type 
star and its atmosphere. 

Thus on one occasion Merrill (1935) wrote that “ the procedure (i.e. 
the ‘ analysis of the light curve of R Aquarii into two curves’) has a 
hypothetical basis: that is to say, it assumes the existence of two light 
sources. It does not, however, postulate the exact character of the 
source of the companion spectrum. This source may be a separate com- 
panion star, itself variable in light, or possibly small areas of high 
temperature within the outer limits of the tenuous red variable”. 

Recently, regarding the same star, he wrote: “* It is somewhat uncertain 
whether the bright lines . . . should be ascribed to the Me star or to 
the nebula. Perhaps the question is without meaning.” (Merrill 1950.) 
Summing up his discussion of this type of star in his monograph, Merrill 
concluded that “* it would be inadvisable at the moment to accept without 
reserve the hypothesis of actual duplicity for all combination spectra’’. 
(Merrill 1940, p. 105.) 


§ 6. CHARACTERISTICS OF STARS WITH COMBINATION SPECTRA 
6.1. Relation between M-Type and ‘ Companion’ Spectra 


Several characteristics of the variations of the light of these stars 
possessing ‘combination’ spectra are in accord with the proposed 
explanation of their origin, in addition to those already described. 

Perhaps the most important is the relationship between the intensities 
of the M-type and ‘ companion’ spectrum observed in R Aquarii by Merrill. 
He writes: “A remarkable fact is that during the years 1931-33, when 
the companion was especially bright, the Me variable had a succession 
of remarkably low maxima. Whether this was more than a coincidence 
cannot be ascertained until the physical nature of the companion and its 
relationship to the red variable are better understood.” (Merrill 1940, 
p. 84.) On the present theory the relationship referred to in that 
quotation is more than a coincidence. It was emphasized in an earlier 
note (Bruce 1955 a) that the electrical discharge theory is in accord with 
the suggestion put forward by Merrill to account for the marked falling 
off in the magnitude of M-type stars at their minimum. He suggested 
that this abnormal decrease is due to the formation of a cloud of particles 
or drops in their atmospheres. The denser this cloud of particles becomes, 


Combination Spectra in Long-Period Variable Stars 1129 


the lower is the magnitude of the M-star at minimum, and the greater 
will be the effect of the processes of static electrification occurring during 
the minimum period. The average electric fields generated throughout 
the atmosphere, during this period will also be greater, and therefore the 
intensity of—or electric currents in—the ensuing electrical discharges 
will also be greater. It is therefore to be expected that abnormally low 
minima of the M-type star will be associated with the abnormal increases 
in the intensity of the electrical discharges and hence of the ‘ companion ’ 
spectrum ; as was observed by Merrill in R Aquarii. 

Furthermore, the existence of these variations in magnitude are in 
accord with the general electrical discharge theory of the cause of this 
type of stellar variation, just as earlier, when applied to the lightning 
discharge (Bruce 1944). the theory of the initiation and propagation of the 
leader stroke accounted for the very marked variation observed in the 
currents in lightning flashes to earth. For the initiation of these long 
atmospheric electrical discharges there are two requisites: First, an 
extensive electric field which is adequate to maintain the discharge once 
it is initiated, and secondly an electric field concentration at one point in 
that general field of sufficient magnitude to initiate the discharge. The 
more uniform the atmospheric conditions in any one cycle, the longer will 
the formation of the region of intense initiatory electric field concentration 
be delayed, and the larger will be the general average electric field through- 
out the whole atmosphere, and hence the greater will be the electrical 
discharges, when they do occur. 


6.2. Effects of Atmospheric Absorption on the Spectra 


Another characteristic of the spectra of long period variable stars, 
which is in good accord with the discharge theory has been observed by 
Merrill in many of their spectra. The discharges originate far down inside 
the stars’ atmospheres and are propagated outwards mainly. As they 
progress, not only are the oxides and other particles vapourized and 
dissociated by the energy liberated in the discharge. so that the ‘ veiling ’ 
is diminished, and the stars increase in magnitude, but the discharges 
emerge from the clouds of TiO and ZrO in the lower atmospheres, and 
ultimately yield light which is little affected by the absorption bands of 
these oxides. This effect is, for example, well illustrated by the variation 
in the intensities of the individual lines of multiplet (2) of Fe I (Merrill 
1947, Table 4) of which Merrill writes: “In the last column, phase + 
162 days (ie. 162 days after maximum light) the relative intensities are 
nearly the same as in the laboratory. At earlier phases the intensities 
are modified, probably by TiO band absorption, as in R Leonis. The 
behaviour of this multiplet presents another example of the general 
tendency of bright lines to escape from the effects of the reversing layer 
as the phase advances.” 

The same explanation probably will account for the well-known varia- 
tion in the relative intensities of Hy : H8 at different phases of the bright- 
ness cycle of these stars (Merrill 1930). It is only some considerable 
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time after maximum light that this ratio approaches the value observed 
in the laboratory, when the lines are freed from the effects of atmospheric 
absorption as the discharge rises in the star’s atmosphere. 


6.3. Rapid Changes in Spectrum Line Intensities 


The occurrence of discharges in these extensive stellar atmospheres 
will result in their atmospheres possessing considerable variations in 
density. Along one discharge channel the current will flow for periods 
of several months, resulting in a considerable change in gas density 
relative to that in neighbouring regions of the atmosphere. When the 
next cycle of discharges occurs the leader strokes will be propagated through 
regions varying considerably in density. As each region of relatively high 
gas pressure is encountered there will result a considerable increase of the 
current density in that region and a consequent brightening of the light 
from the discharge. As distances of the order of one to ten million kilo- 
metres are covered in a day by the advancing head of the discharge 
column, the effect of changes in density can be considerable in times of 
this order. Rapid changes in line intensities are observed in many of 
these stars, and particularly well in BF Cygni (Merrill 1943). 


§ 7. CONCLUSIONS 


1. The combination spectra observed in several long period variable 
stars can be explained on the basis of the electrical discharge theory of 
their variability. 

2. On this theory the ‘companion’ spectrum is emitted from the 
electrical discharges at relatively low levels in the star’s atmosphere. 
The nebular spectrum containing forbidden lines is emitted when the 
discharge reaches the extremely low pressures of the outer regions of the 
star’s atmosphere. 

3. The theory explains the observed time interval between the maxima 
of these two spectra. 

4. It explains the observed relationship between the intensities of the 
M-type and *‘ companion’ spectra observed by Merrill in R Aquarii. 

5. The theory explains the variation in the relative intensities of the 
individual lines of multiplets as the outburst progresses. 
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CXXVI. CORRESPONDENCE 


The Decay Scheme of 1*La 


By C. F. CoLEMAN 
Atomic Energy Research Establishment, Harwell 


[Received July 4, 1955] 


AN investigation of the decay scheme of 14°La, prepared by slow neutron 
irradiation of Johnson Matthey Specpure La,O,, has been carried out 
with sodium iodide scintillation spectrometers. Gamma rays of energies — 
328, 486, 815, 893-+40, 1597, and 2570+150 kev were found, with 
relative intensities of 25-10, 48+8, 27+5, 1144, 100, and 4+1. No 
others were observed with energies greater than 50 kev, and in particular 
the intensity of the 2900 kev radiation reported by Wattenberg (1947) 
was found to be less than 1% of that of the 1597 kev radiation. The 
only previous report of a gamma ray with an energy of about 900 kev 
was made by Cork et al. (1951), who observed conversion electrons of a 
920 kev radiation. Since the ratio of the heights of the 815 and 893 kev 
photopeaks in the lanthanum spectrum changed by less than 20% 
over a period of three half-lives, the 893 kev radiation is assumed to belong 
to the 14°La decay, and is identified with Cork’s 920 kev radiation. 

Coincidence studies were then carried out by setting a single channel 
pulse analyser on each photopeak in turn, and observing the coincident 
spectrum in a second counter by means of a thirty channel pulse analyser. 
The 1597 kev line was found to be in coincidence with all the lines of 
lower energy, including the 920 kev line (fig. 1). The other coincidence 
spectra showed that coincidences also occurred between the 328 and 486 
kev gamma rays, and that, within the experimental error caused by 
the presence of the 1597 kev Compton distribution below the photopeaks 
of the lower energy gamma rays, both the 815 and the 920 kev lines 
were in coincidence with the 1597 kev line alone. These results, together 
with the beta-gamma coincidence measurements of Robinson and 
Madansky (1951) and the beta-ray energy measurements of Peacock 
et al. (1954), strongly support the decay scheme given in fig. 2, proposed 
by the latter group of workers. 

Values obtained for the angular correlation of the 815, 1597 kev 
gamma ray cascade have differed widely (Robinson and Madansky 1951, 
Bishop and Perez 1955), probably because the observed correlations 
include different contributions from the 920, 1597 kev cascade. The 
correlation of this latter cascade was isolated and measured by feeding 
the output of one detector to a single channel pulse height analyser 
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arranged to accept as much of the 920 kev photopeak as was consistent 
with complete rejection of the 815 kev photopeak. A rate controlled 
gain stabilizer was used to minimize the effect of amplifier and multiplier 
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gain drifts. The observed correlation, after correction for finite angular 
resolution, is given by 
W(0)=(1-000-40-010) + (0-200-+0-035)P,(cos 0)— (0-002-+.0-043)P,(cos 8) 

which can be fitted by the spin sequence for the 2517 and 1597 kev 
levels, 2{M1+-(0-1% to 1-0) E2}2{H2}0, where the M1 and E2 components 
in the first transition are in antiphase according to the convention of 
Biedenharn and Rose (1953). The spin shown in fig. 2 for the 2083 kev 
level is that which Bishop and Perez derive from their measurements on 
the angular correlation of the 486, 1597 kev cascade. 

From the above decay scheme over 95% of all 14°La disintegrations 
are accompanied by emission of the 1597 kev gamma ray, and fewer 
than 5°, by the emission of any other gamma radiations with energies 
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sreater than 920 kev. Thus it should be possible to estimate the absolute 
disintegration rates of 14°La sources to within 5% by beta-gamma 
coincidence measurements if a thallium activated sodium iodide crystal 
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is used as the gamma ray detector, and only those pulses are accepted 


which correspond to an energy release in the crystal greater than about 
1000 kev. 
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Electric Quadrupole Lifetimes in Even—Even Nuclei 


By C. F. Coneman 
Atomic Energy Research Establishment, Harwell 


[Received July 4, 1955] 


Two groups of coincidence measurements have been carried out in an 
attempt to determine the lifetimes of a number of E2 gamma ray tran- 
sitions. In the first group, values of mean lifetimes were obtained from 
the side slopes of coincidence curves with a full width at half height of 
2x10-* seconds, taken using T.P.B.D. phosphors, E.M.I. 6260 photo- 
multipliers, and 30 Mc/sec fast amplifiers in a parallel coincidence arrange- 
ment. With this equipment mean lifetimes = 6 x 10~1° seconds could be 
detected. In the second group, for which the corresponding limit was 
~3 x 107-1! seconds, the lifetimes were deduced from the centroid shifts 
(Bay 1950) between pairs of beta-gamma coincidence curves with effective 
full widths of 1-6 x 10~° seconds, taken under comparable experimental 
conditions, using T.P.B.D. phosphors, R.C.A. 5819 photomultipliers, 
and 200 Mc/sec fast amplifiers, again in a parallel coincidence arrange- 
ment. The differences (~5 x 10-1! seconds) between the flight times of 
beta rays and gamma rays were calculated from the energy bands selected 
by the slow channels (widths ~20°) and the geometry of the equipment. 
In most cases the paired curves were taken from a single source (‘ self- 
comparison ’), a beta absorber being placed before each detector in turn 
to exchange the functions of the beta and gamma counters, but the 
coincidence curves from 7As and 124Sb were paired with comparable 
curves from 24Na. The experimental technique used in these measure- 
ments will be described in detail in a later paper. The quoted uncertain- 
ties in the mean life-times for this group arise from the different distribu- 
tions within the phosphors of the points of origin of light photons produced 
/ by beta and by gamma irradiation. The statistical uncertainties are 
usually less than 10-1 seconds, going down to 2x 107 seconds for the 
24Mo measurement. Where two gamma rays are emitted in cascade* 


Nucleus Tobs E,(mev) Texp 7H2(S.P.M.) | M|2 Method 

24 gr (842) x10! sec_| 1-38* | (3:6 2°5) x 10-11 see | 1-6 10-11 sec | 20-25 Self Compn. 
ai <2 10-41 sec | 0:89* <4 10-1 sec | 5:°9X10-1! sec | >1°5 Self Compn. 
56He 0-85 <6 10-19 gee | 6:0X10-1! sec | >0-10 

1:26 <3 10-1 sec | 8:3X 10-12 sec | >0-02 
6°Ni Se N= Bere || lalla. <4x10-!gec | 1-0X10-' sec | >0-25 Self Compn. 
16Se (342) x 10-11 sec | 0:56* | (3342-2) x 10-11 see | 3:0 10-19 sec | =5 Compn. *4Na 
Be Sr 1-85 <6 x 10-19 gec | 6-9 x 10—-™ see | >0-001 
B24Cd 0:55 <6 x10-1% sec | 2:0X10-!° sec | >0:3 
B2*"T'e (20-2) x 10-11 see | 0-57* (22+3)x10-11 sec | 1-4x10-7° sec 0-64-£0-1| Self Compn. 
imTe “<2 10-11 sec | 0-60 <2 10-11 sée | 1:1X10-2%see | >5 Compn. *4Na 
BSCS 0-122 | (5-9-40-8) x 10-® sec — 


‘Tobs—See text. 

texp—Mean lifetime deduced from experiment. 

7k2(S.P.M.)—Theoretical mean lifetime for E2 radiation deduced from 
the Single Particle Model. 

|M|-2—Ratio texp/TH2(S.P.M.). 
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a value for the mean life-time of the softer transition has been calculated 
from the experimental average of the two lifetimes, 7,,,, by taking the 
lifetime of the harder transition as zero. The metastable level found in 
131('s is fed by the 0-494 Mev gamma ray from the decay of "Ba, con- 
firming the decay scheme proposed by Lu et al. (1955). 

Most of the data presented here are for transitions between the first 
excited states and ground states of even-even nuclei. It will be seen 
that several of the observed lifetimes are comparable with those estimated 
from the single particle model, or even smaller, in two cases by a factor of 
five, giving strong evidence for the existence in these nuclei of collective 
motion (Goldhaber and Weneser 1955). 
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The X-ray Spectra of u-mesonic Copper and Lead 


By F. D. 8. BuTEMENT 


Atomic Energy Research Establishment, Harwell 
[Received July 27, 1955} 


THE theoretical significance of measurements of the 2P-—IS radiative 
transition energies of .-mesonic atoms of the heavier elements (from which 
the nuclear radius may be calculated) has been discussed theoretically by 
Wheeler (1953), Cooper and Henley (1953), and Hill and Ford (1954). 
Experimental measurements for a number of elements have been made by 
Fitch and Rainwater (1953). This letter describes measurements made 
on copper and lead. 

The arrangement of the apparatus is shown in fig. 1. The external 
7 and ~~ meson beam of the University of Liverpool cyclotron (after 
passing through a channel in the 12 ft. thick iron loaded concrete shielding 
and deflection by an analysing magnet) was filtered through 48 mm of 
copper (1) in order to stop the ~ mesons. The emergent beam of 
# mesons and fast electrons in approximately equal proportions passed 
through a collimator in a wall of iron loaded concrete (2), then through 
scintillation counters with plastic scintillators 8 mm thick and 63 mm 
diameter (3 and 4) and fell on a slab of the target element 75 mm square 
and 10 g cm ? thick (5). The mesonic x-rays emitted were detected by a 
Nal(Tl) crystal 44mm in diameter and 47 mm thick (7) viewed by a 
type 6260 E.M.I. photomultiplier tube. Events were selected by coin- 
cidence between counters 3, 4 and 7, with a resolving time of 3x 1077 sec. 
Counter 6 with a plastic scintillator 6 mm thick and 8 cm square operated 
as an anticoincidence counter to reject events due to a charged particle 
traversing 5 and entering 7. About 400 ~ mesons per minute stopped in 
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the target. The pulse height distributions from counter 7 was analysed 
on a 100 channel analyser similar to that described by Wilkinson (1950). 
Energy calibrations were made with the 1:33 Mev y-ray of ®Co and the 
4-43 mev y-ray from a Po—Be source. The attenuator scale of the spectro- 
meter amplifier and the linearity of the pulse analyser were calibrated and 
found to be accurate to better than one part in a thousand. The main 
source of error was due to slight drifts in the gain of the photomultiplier, 
but numerous calibrations showed that this did not exceed 1%. 


Fig. 1 


& a? ley 


Arrangement of apparatus. 


The mesonic x-ray spectra observed were superimposed on a continuum 
of pulses extending to above 10 Mev, due to bremsstrahlung produced by 
high energy electrons. Nuclear y-rays following various modes of nuclear 
disintegration caused by eventual nuclear capture of the mesons must also 
have been recorded, but these y-rays would have individual intensities 
considerably less than those of the mesonic x-rays, and were not observed. 

The spectrum of copper (fig. 2) shows the (unresolved) 2P-1S transition 
as a photoelectric peak at 1-55 Mev together with the corresponding 
Compton peak at lower energy. Fitch and Rainwater (1953) also obtained 
a value of 1:55 Mev. 

Figure 3 shows the 2P-1S transition of lead. At high y-ray energies 
(2) there is a peak due to pair production at (H— 1-02) Mev, a peak of almost 
equal intensity at (7—0-51) Mev due to pair production with subsequent 
absorption of one of the annihilation quanta and a weak peak at the full 
energy due partly to photoelectric absorption and partly to pair produc- 
tion with absorption of both annihilation quanta. The resolution on a 
single peak was 8% (width at half maximum intensity) and it was there- 
fore not possible to resolve the 2P,).—-18 and 2Ps),-18 components, which 
should theoretically have a separation of about 3-6%, for the pair produc- 
tion peak (assuming that the , meson is a Dirac type particle of spin 3). 
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Fig. 2 


ENERGY MeV 


Spectrum of »-mesonic copper. 


ENERGY MeV 


Spectra of 4-mesonic lead and of 4:43 Mev y-ray. 
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Study of the peak profile for the 4-43 Mev y-ray from a Po-Be source, 
obtained with good counting statistics, indicated that the most reliable 
method of determining the energy of a statistically imperfect profile was to 
determine the centre of the composite (Z— 1-02) and (H—0-51) Mev peak at 
a level of 80°, of the maximum peak height. The value thus obtained 
was very close to the mean, i.e. (H#—0-77) Mev. Application of this method 
to the lead spectrum gave a value of 5-84 mev for the unresolved 2P-1S 
transition. Hill and Ford have shown that for any plausible nuclear 
model the 2P3), and 2P,). levels have a separation of between 0-176 and 
0-186 Mev and relative statistical weights of 2:1. With the resolution 
quoted above, the 2P3).-1S transition should then be approximately 
0-05 Mev higher than the unresolved peak, i.e. at 5-89 mev. Fitch and 
Rainwater found 6-02 Mev. 


The author wishes to thank Professor H. W. B. Skinner for providing the 
facilities of the University of Liverpool cyclotron, and Mr. D. West and 
Mr. E. F. Bradley for help with the operation of the apparatus. 
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A Note on the Ground State Spin of 1 Hafnium 


By B. Hartmann and T. WIEDLING 
Laboratory of Physics, University of Stockholm, Sweden 


[Received July 28, 1955] 


RASMUSSEN made spectroscopic fine structure measurements on peeett 
and reported a ground state spin of <3/2 (Rasmussen 1935). It is 
well known that angular correlation measurements are a method of 
determining spins not only of excited states but under favourable condi- 
tions of ground states as well. Some experiments of this kind were made 
on 177Hf some years ago in a trial to exclude one of the proposed spins 
given above (McGowan et al. 1951, Wiedling 1952). According to the 
results obtained it seemed possible to exclude spin 1/2. . 

We have now made a reinvestigation of the angular correlation of 
177Hf with the use of some different chemical compounds. The measure- 
ments were made on the two strong gamma lines in the well-known decay 
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scheme of 177Hf (Douglas 1949, Marmier and Boehm 1955). The investiga- 
tion was made with two scintillation spectrometers coupled in coinci- 
dence. One spectrometer accepted only the full energy photopeak of 
113 kev and the other spectrometer only the peak of 208 kev. 

Measurements were made on the following chemical compounds : 
Lu,O 3, Lu(NO,), as a powder, Lu(NOs), in water solution, and Lu(NOs)s 
in water solution mixed with glycerine. 

The two sources mentioned last were chosen because of the fact that the 
attenuation of an angular correlation is a function of the viscosity of the 
solution under some circumstances. The angular correlation measure- 
ments were made on several sources of each compound. Data were taken 
for every fifteen degrees. A least squares fit of the experimental points to 
the function W(@)=1+A4,P, (cos @) was made. The results are shown in 
the table. 


Source —A, 


Lu(NO3)5 0-142 
(powder) 

Lu(NO3)3 0-176 
(in water 
solution) 

Lu(NOs3)3 0-140 
(25°% water and The errors of A, 
75% glycerine) are about +4% 


Within the experimental errors the value of A, obtained here for Lu,O, 
is In agreement with earlier results. However, the less pronounced 
anisotropy obtained for Lu( NOs), as a powder and for Lu(NOg3), in water 
solution mixed with glycerine is remarkable as well as the stronger 
anisotropy found for a dilute solution of Lu(NOs), in water. This last 
mentioned anisotropy may be a maximum because of the fact that solu- 
tions often give maximum anisotropy. The fact that the attenuation of 
the angular correlation is less pronounced when the crystalline source is 
dissolved in water shows that the source effect may be an interaction of 
the electric quadrupole moment of the first excited level of the nucleus 
with an electric field gradient. 

As pointed out earlier Rasmussen proposed a ground state spin <3/2. 
The shell model prefers spin 1/2. If we suppose the spin sequences of 
the levels of '’’Hf to be 5/2-+7/2-+3/2 or 3/2+5/2—+1/2 and assume the 
emitted gamma-rays to have pure character the angular correlation theory 
gives A,=—0-153 and A,——0-200, respectively, compared with our 
maximum value —0-18. Thus our measurements show that at present 
it is not possible to draw any certain conclusions about the ground 
state spin of Hf from angular correlations measurements but that it is 
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necessary to use other experimental methods to get information about this. 
It seems, however, to be possible to get other valuable information about 
the ’”Hf nucleus from correlation measurements and for this reason our 
investigations continue. 
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Screw Dislocations and Polytypism of Cadmium Todide 


By Ricuarp 8. MitcHeLn 
Department of Geology, University of Virginia, Charlottesville 


[Received July 28, 1955] 


INTRODUCTION 


SINCE the appearance of a previous summary on polytypism of cadmium 
iodide (Mitchell 1954), many new forms of this substance have been identi- 
fied, and the stacking structures have been determined for several of 
these. In this note the writer plans to summarize this new information 
and then draw some conclusions concerning the relationship of screw 
dislocations to the formation of these polytypes. In order to avoid 
repetition of the previous summary of this subject no space groups or 
unit cell constants are given here for the different polytypes. All of the 
determined stacking structures are summarized, however, since these 
are essential to the discussion of the dislocation theory. Each structure 
is represented both by the classical A, B, C notation and the iodine atom 
zigzag sequence. Details on the determination of these structures will 
be published in the future. The crystals used in this investigation were 
grown from aqueous solutions at room temperature, and were picked at 
random from these solutions for study by the Weissenberg x-ray method. 


SUMMARY OF POLYTYPES 


The writer has encountered a few crystals of type 2H. Lach of these 

has the structure [(AyB)],, (11), which was first determined by Bozorth 
1922). 

: ane 4H is by far the most common modification of cadmium iodide. 
This type has the stacking structure [(AyB) (C«B)],, (22). Several 
earlier investigators reported cadmium iodide crystals with this structure. 

Two specimens of 6H were recently found. Both of these have a 
[(AyB) (CB) (AyB)],, (2211), structure. This stacking sequence differs 
from Pinsker’s (1941 a, b) 6H type which has the arrangement [(AyB) 
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(CBA) (CaB)],, (33). Theoretical considerations for the existence of 
these two types will be given later in this note. In the future the writer 
will designate the (2211) type as 6H,,) and the (33) type as 6H,,). This 
follows a scheme used by Ramsdell and Kohn (1952) for naming two 
polytypes of silicon carbide with equal numbers of layers but different 
stacking. 

Two crystals of type 8H have been noted. One of these, reported 
previously, has the structure [(AyB) (CB) (AyB) (AyB)],, (221111). The 
second example appears to have an identical structure, but, because it is 
syntactically coalesced with 4H, a good evaluation of the intensities of 
some of the x-ray reflections is impossible. 

Type 10H has been encountered twice. Both crystals have a structure 
represented by [(AyB) (CaB) (AyB) (CaB) (AyB)],, (222211). 

Three structural variations of type 12H have been identified. Type 
12H,,, has the structure [(AyB) (CaB) (AyB) (CaB) (AyB) (AfC)],, 
(222123). In this particular case the iodine zigzag sequence is not, 
by itself, a sufficient symbol for designating the complete structure. The 
zigzag arrangement does not fix the cadmium atom positions when there 
is no horizontal plane of symmetry for the iodine atoms, and, hence, the 
cadmium atoms may lie on one of two groups of positions. The A, B, C 
notation, however, gives an unambiguous designation. Type 12H,,) has 
the structure [(AyB) (AyB) (C«B) (Ca«B) (AyB) (ABC)],, (22121121). 
The structure of type 12H,,) appears to be [(AyB) (CaB) (AyB) (CaB) 
(AyB) (AyB)],, (22221111), but this is still uncertain because the crystal 
is syntactically coalesced with 4H which makes the intensity evaluation of 
certain x-ray reflections difficult. 

One crystal of 14H was identified. This has a structure represented by 
[(AyB) (CaB) (AyB) (CaB) (AyB) (CaB) (AyB)],,, (22222211). 

In addition to the aforementioned crystals, the writer has also observed 
types 16H, 18H, 20H, 24H ,); 24H,,), 24H, 24H), 26H, 28H(,,), 28H), 
32H), 32H,,), 40H,,, 40Hy), 40H, 44H, 50H, 56H, 62H and 64H, and 
is working on the structures of some of these at the present time. There 
is still some uncertainty about the identity of a rhombohedral type, 24R, 
and several other types consisting of many layers. 


EXPLANATION OF POLYTYPISM 


This study has not only proved Forty’s (1952) idea that cadmium iodide 
exhibits polytypism, but also now virtually proves his explanation for 
this phenomenon. He suggested that various slips in the [0001] direction 
of the ‘ideal structure’ 4H, [(AyB) (C«B)],, would produce different 
kinds of steps on the (0001) surface, and the crystal generated by spiral 
growth on this surface would have a structure controlled by the structure 
of these steps. In fig. 1 atoms below line ABCDEF are arranged according 
to the 4H stacking sequence (JKLM outlines a 4H unit cell). This struc- 
ture has been dislocated in the [0001] direction along BG and DH. BC 
represents a slip of the structure equal to a Burgers vector of unit strength 
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(equal to one sandwich of atoms or 6-835 4), and DE is a similar slip equal 
to a vector of strength three. The structure at step BC is (AyB) or 
(11). Ifthe crystal grows on this step alone a 2H type would be generated. 
A unit cell of this is outlined by NPQR in fig. 1. Step DE has the 
structure (AyB) (C«B) (AyB) or (2211). If this step existed by itself 
type 6H (STUV) would be generated. One should also notice how well the 
atoms fit together along the planes of dislocation. Following is a sum- 
mary of the structures that would be produced by single slips of different 
magnitudes of the type illustrated below. 


% Fig. 1 OA 
OER AIS 
LO AIT AACACEGACE GA 
dears ee 


Creve eee tamara 


A (1120) cross section of a dislocated 4H cadmium iodide structure. White 
circles represent iodine atoms while black circles represent cadmium 
atoms. The atoms and unit cell dimensions are drawn to scale relative 


to each other. 


It should be noted that every time the strength of the Burgers vector is 
equal to 2n the 4H structure results, and, therefore, a single dislocation of 
this type would never produce an 8H, 12H, 16H . . . 4nH polytype. 
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From the theory, one might also expect to encounter type 4H more 
frequently than other modifications. 

A close examination of the experimental results, outlined in the earlier 
part of this note, will reveal that experimental observation completely 
supports the above theory with a few modifications. Types 2H, 4H, 
6H,,), [OH and 14H have been found and have structures identical with 
those predicted by the theory. In addition type 4H is by far the most 
abundant polytype of this compound. 


Burgers ate a Bh . Polytype 
vector strength Bestisins seer so ug generated 
1 (AyB) 2H 
(11) 
Z (AyB)(CaB) 4H 
(22) 
3 (AyB)(CaB)(AyB) 6H 
(2211) 
4 Coy a oe a 4H 
5 Os 10H 
6 (AyB)(CaB)(AyB)(CxB)(AyB)(CaB) 4H 
(222222) 
7 (AyB)(CaB)(AyB)(CaB)(AyB)(CaB)(AyB) 14H 
(23222211) 
etc. 


Polytypes 8H and 12H,,, reported above, can easily be explained by 
assuming the presence of two simple screw dislocations of the same 
hand. In fig. 1 it was illustrated how 2H and 6H each originated by a 
single dislocation. Now if both dislocations occurred in the same crystal 
and close enough together to co-operate (Frank 1951), the final effect 
would be an addition of the two step structures, resulting in a step with 
the structure (AyB) (AyB) (C«B) (AyB), (112211). This is illustrated in 
fig. 1 where a unit cell of the structure is outlined by WXYZ. This is 
the structure that was determined experimentally for 8H. A step needed 
to generate the 12H,,) structure will result from a similar combination of 
two dislocations whose steps have the 2H and 10H structure. 

To explain types 6H,,), 12H,,) and 12H,,) one must consider a further 
complication of the theory. In all the cases up to this point one will 
notice that when atoms began to take their place on the newly formed 
step they fell into the position of atom g in fig. 2, which is the position to 
be expected by the control of the exposed step. However, if the first 
atom fell into position h certain variations in the growing crystal would 
result. Also it should be noted that the arrangement of neighbours 
around atom h is similar to arrangements in the dislocated region 
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along BG. A structure of this kind, built upon a dislocation with a Burgers 
vector of unit strength, would give an (ABC), (11), arrangement which is 
equivalent to the simple 2H structure. A unit cell of this type is outlined 
by NPQR in fig. 2. Likewise, a structure of this type built upon a dis- 
location with a vector strength of three would produce a 6H type with the 


structure [(AyB) (C«B) (ABC)],,, (33), which is equivalent to the 6H, 
structure reported by Pinsker (1941 a, b) 


Similar to fig. 1, but showing a different arrangement of atoms growing on the 
dislocation step BC. 


The 12H, structure [(AyB) (CaB) (AyB) (CaB) (AyB) (ABC)],, is equiva- 
lent to the simple 10H structure plus one (ABC) sandwich. <A step with an 
arrangement necessary to generate this stacking could be formed by the 
combination of two dislocations of the same hand. One of these would 
have a dislocation strength of one and would grow by adding the first 
atom in the h position. The other dislocation would involve a slip equal 
to a vector strength of five and would grow by adding atoms in the normal 
g position. Type 12H,,) could result from the combination of five small 
dislocations, the growth of three involving atoms in the g position and the 
growth of the other two involving atoms in the h position. This does not 
seem to be unreasonable. Forty (1952) has stated that if a crystal plate 
buckles fairly uniformly the resultant slip will produce several screw dis- 
location groups of about the same size. 

Because polytypes can be generated both from multiple dislocations, 
and growth in g and h positions, it is easily understood why one polytype 
may have several structures. The fact that type 26H, for example, fits 
in the series 2H, 6H, 10H, 14H, etc., but does not have the structure 
(22222222292211) of this series, does not invalidate the theory. Its 
structure, still undetermined, could be one resulting from multiple 
dislocations. 

It is interesting to note that Forty (1952) measured spiral step heights 
on cadmium iodide and found values ranging in magnitude corresponding 
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all the way from 16H to 96H. Values corresponding to types less than 
16H were not determined because of limitations in the method of measure- 
ment. Some of his steps no doubt occurred on pure 4H crystals which 
had one dislocation with a Burgers vector strength equal to 2n, while 
others were no doubt present on polytypes of many layers, yet to be 
discovered. He found steps equal in value to the writer’s types 16H, 
24H, 28H, 32H, 40H, 44H, 50H, 56H, 62H, and 64H. 

The writer feels that the observations presented in this note reasonably 
prove the theory that polytypism of cadmium iodide is a phenomenon 
that results directly from screw dislocations of a variety of magnitudes 
in the ‘ ideal structure’, 4H. In addition to explaining polytypism, this 
work also seems to validate the more general theory that crystals grow 
from screw dislocations. A more detailed discussion of this work will 
be published in the future. 
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Ionization at the Origin of Electron Pairs, and the Lifetime of the 
Neutral Pion 


By D. H. PERKins 
H. H. Wills Physical Laboratory, University of Bristol 


[Received August 12, 1955] 


Ir was first suggested by King (1950) that, as a result of the close proximity 
of the trajectories of the unlike charges of the electron and positron created 
in the materialization of a high-energy photon, the ionization near the 
origin of such a pair might be less than that of the two particles separately. 

During an investigation of the soft component associated with high 
energy * jets’ produced by cosmic rays in nuclear emulsions, preliminary 
results had suggested that the lifetime of the neutral pion, as measured by 
Anand (1953 a), might be an overestimate. Anand’s method consisted 
essentially in measuring the length of the gap between the centre of a 
star and the first clearly resolved grain in the track of a ‘ direct ’ electron 
pair originating from the alternative mode of decay 79-+e++e-+y. 
The gap distribution was compared with that of the tracks of protons 
of twice-minimum ionization projected from the same stars, and it was 
a possible point of objection that the method depended on the assumption 
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of a constant ionization along the path of the ‘ direct ’ pair. It was there- 
fore decided to make measurements of the ionization along a number of 
tracks of pairs of energy of the order of 100 Bev, for in such conditions the 
effect suggested by King, if indeed it existed at all, might be present for a 
considerable distance from the pair origin and go be accessible to measure- 
ment. 


o so 10 200 300 400 500 600 700 800 900 1000 


DISTANCE FROM PAIR ORIGIN, MICRONS — 


In the figure is displayed the result of measurements on seven energetic 
electron pairs, the coefficient g of the gap-length distribution (Fowler and 
Perkins 1955) being plotted against distance x from the pair origin. The 
full line in the graph refers to the mean value of g for x between 0:3 and 
1-0 mm, and the reduction in the value of g, and hence of ionization, for 
values of x<0-3 mm is at once apparent. The energies of the pairs 
were estimated from the relative scattering of electron and positron over 
the first radiation length or so, and from the features of the soft cascade 
at greater distances. They lay within the range 80—400 Bev, with a mean 
value of the order of 180 Bev. At such energies, the relative scattering of 
the two particles is quite small, and the tracks to which they give rise 
are unresolvable for a distance of several mm from the origin. 

Theoretically, the reduction in ionization near the origin must pre- 
sumably arise from the fact that for large impact parameters between the 
pair and the target electrons of the emulsion, the incident field is that of 
an electric dipole, whereas at greater distances from the origin, the electron 
and positron are well separated and behave as two independent monopoles. 
The ionization is expected to increase with x approximately as log d/dmax, 
where d is the lateral separation of the electrons, roughly proportional to 
x, and dmax is the maximum effective value of the impact parameter for 
momentum transfers contributing to the activation of the silver bromide 
crystals traversed. The dotted curve indicates the form of the expected 
variation of ionization with distance, calculated from the opening angle 
of the pairs and for an arbitrary value Ol ax 0 Oe Cle 
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As stated above, the mean energy of the pairs examined was 180 Bev. 
whereas that of the direct pairs in Anand’s experiment was about 200 Mev. 
Thus, assuming the relation between opening angle and energy for direct 
pairs to be similar to that of pairs produced by external conversion of 
photons (Anand 1953 b), the reduction in ionization in Anand’s experiment 
would only be serious within less than 0-5. of the pair origin. This 
distance is small compared with the mean gap-length (~2-5 ;) in a twice- 
minimum track. What is more important, the proximity of the parent 
star very often obscured the first micron or so of the path of the direct 
pair. Taking these two factors into account, it turns out that, although 
Anand’s assumption of a constant ionization along the path of the pair 
must lead to an overestimate of the lifetime of the neutral pion, the error 
introduced is very small in comparison with the statistical errors of his 
measurements. 

Recently, an independent estimate of the lifetime of the neutral pion, 
based on a modification of the intercept method of Carlson, Hooper and 
King (1950), has been carried out in this laboratory (Brisbout et al. 1955). 
The value so deduced falls with a 90% probability within the limits 1-2 
and 5x 10~1° sec, in fair agreement with Anand’s value of between 2-5 
and 14% 10-1 sec for the same confidence limits. 

If the effect described above can be considered as established, it offers 
a useful alternative method of estimating the energy of y-rays in the region 
above 101" ev, subject of course to large statistical fluctuations. At such 
energies, the conventional measurement of pair energy based on the 
relative Coulomb scattering of the two particles, has to be carried out 
over distances greater than or comparable with one radiation length 
(3 cm). Thus confusion arises because of the presence of bremsstrahlung 
pairs, and recourse must be had to estimation of the energy from the 
somewhat tedious analysis of the ensuing soft cascade. 
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On the Magneto-resistance Effect in Metals 


By Rosemary A. Conpwein-Horsraty and D, Ter Haar 
Department of Natural Philosophy, St. Salvator’s College, St. Andrews 


[Received August 15, 1955] 


Onz often finds in the literature the statement that the magneto-resistance 
effect in metals will be zero, if the one-band model is used. This statement 
is only true, if one understands it to mean that magneto-resistance is really 
a second-order effect in the same way as are the thermo-electric power and 
thermal conductivity. Indeed one finds in Seitz: (1940) a formula for the 
magneto-resistance derived on the basis of the one-band model which is 
stated to be correct in the general case of an energy-dependent mean free 
path, and which is quoted from an earlier paper by Sommerfeld and Frank 
(1931). However, Sommerfeld and Frank’s derivation is only correct in 
the case of a constant mean free path. It therefore seemed to us to be of 
some interest to derive a general formula for the magneto-resistance in a 
metal with an energy-dependent mean free path derived from the one-band 
model in the usual way. We have found, up to terms of higher order, 


[o(H)—o(0)]/o(0) = —DH?/(1+-CH?) (1) 
in which 
D=B(1—$ [d In Id Ine’},,_..)?, (2) 


where all symbols have the same meaning as in Seitz’s book.* 


From eqns. (1) and (2) it is clear (i) that the magneto-resistance calcu- 
lated in this way will never be of the same order of magnitude as_ the 
observed effects which is the reason why Jones (1936) and Sondheimer 
and Wilson (1947) developed the two-band model ; and (ii) that the effect 
always has the conventional sign and can thus not be used to explain the 
recently found effects of opposite sign in In-Sb (Lark-Horovitz, private 
communication) and Sb (Steele 1955). 
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exponent of z should be 2, and a factor ¢ is missing from the denominator of the 
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Slip in the Interior of Aluminium—Silver Single Crystals 


By N. A. McKinnon 


Aeronautical Research Laboratories, Australia 
[Received August 2, 1955] 


Tue mechanism of deformation in the interior of metal crystals has given 
rise to much speculation in recent years (Brown and Honeycombe 1951, 
Kuhlmann-Wilsdorf et al. 1952). However a reliable method for disclosing 
the distribution of slip bands in the interior of, for example, super purity 
aluminium single crystals has not yet been found. The author has been 
investigating a possible technique which involves the addition to super 
purity aluminium of a solute element, which, after deformation has been 
carried out in the solution treated and quenched condition, will pre- 
cipitate at a suitable temperature on any dislocations remaining in the 
crystal. Subsequent removal of the surface by electropolishing and 
conventional metallographic procedures enables the location of the pre- 
cipitate and hence the distribution of dislocations to be found. 

An alloy of aluminium with 3-8 wt. % of silver was used. Single 
crystals with a gauge length of 2-3in. and cross section 0-3 in. x 0-1 in. 
were grown by a modified strain anneal procedure, electropolished. 
solution treated, quenched at a rate not fast enough to cause detectable 
quenching strains, and then quickly repolished and strained in tension. 
Part of one surface of the single crystals used had been abraded before the 
final polish during a study of the effects of such abrasion on the nature of 
slip bands, but in view of later work to be described separately it is not 
believed that this abrasion had any effect on the results to be detailed. 

A typical surface micrograph after an elongation of 5% of erystal LG103 
whose orientation, fig. 1 (Plate 26), is towards the middle of the stereo- 
graphic triangle is shown in fig. 2 (Plate 26). Under oblique illumination 
it was found that some slip bands were more prominent than others. 
No slip bands on other systems could be detected. After ageing and 
repolishing with an accentuation of the structure by an anodising treat- 
ment, the interior micrograph of fig. 3 (Plate 26) was obtained. 

This shows heavily developed primary slip bands, obviously related to 
the prominent bands of fig. 2, (Plate 26) but no trace of the other more 
closely spaced surface slip bands between the prominent ones. Whether 
closely spaced bands exist in the interior is the subject of further investi- 
gation. For the moment, the important feature of fig. 3 (Plate 26) is 
the indication in this crystal of ‘ middle’ orientation of short slip bands 
on the conjugate (Rosi and Mathewson 1950) system lying across the 
primary system. The slip bands on the conjugate system became con- 
tinuous as the edge of the crystal was approached and here evidence was 
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also found of slip on the cross slip system (fig. 4 (Plate 27)). This effect 
at the specimen edges has not been found for other crystals—it may be 
related to the fact that the glide direction was towards a corner of the 
test specimen and not, as in the other cases, towards the surface being 
examined. 

Crystal LG81, of orientation near the [100]-[110] side of the stereo- 
graphic triangle, showed only primary slip bands on the surface after an 
extension of 2°%, but in the interior short slip bands were again found on 
the conjugate system crossing the primary one (fig. 5 (Plate 27)). 

For crystal LG80 on the [111]-[{110] side of the stereographic triangle, 
primary slip only was observed on the surface after an extension of 7%. 
Tn the interior no evidence was found of any aggregation of dislocations on 
other than the primary system (fig. 6 (Plate 27)). 

None of the 3-8°% silver crystals so far tested have shown any deforma- 
tion bands of the S-type either on the surface or in the interior but bands 
of secondary slip of the step ladder type (Calnan 1952) are sometimes 
observed. The stress strain curves for the particular crystals described 
were not determined, but in more recent experiments, the shape of the 
curves has been consistent with linear hardening, the slope being approxi- 
mately the same as that for super pure aluminium crystals of similar 
orientations. 

Part of the explanation of these results is believed to be that at the 
applied stress which is necessary for several per cent elongation on the 
primary system the total resolved stress on other less favourably oriented 
systems due to the applied stress plus that set up by groups of dislocations 
on the primary system, is large enough to cause the generation on such 
a system of a limited number of dislocations. These dislocations are 
apparently not sufficient in number to produce a detectable step on the 
surface. Experimental evidence for the generation of slip bands in the 
interior prior to the applied stress exceeding the conventionally determined 
yield stress has been reported in sodium chloride by Pratt (1952). 

For crystal LG81 where short slip bands on the conjugate system were 
observed the second most favoured slip system is the conjugate system. 
For crystal LG80, on the other hand, where only the one set of slip bands 
was found in the interior the slip system with the second highest resolved 
shear stress is a different direction in the same slip plane and this of course 
cannot be distinguished from the primary system. Crystal LG103. 
though further from the [111]{110] boundary than LG80, shows a 
similar order of favoured slip systems. However the interior metallo- 
graphy results differ from those for LG80 in that slip on the conjugate 
system which is the third favoured one was actually found in the interior 
while at the edges of the crystal some slip on the cross slip system was 
also detected. The specimen shape and the orientation of glide plane 
and direction towards that shape as they affect the rotation of the gauge 
length relative to the shoulders during plastic extension is being further 


investigated. 
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The detection of less favourably oriented slip systems in the interior of 
crystals for which only primary slip was observed on the surface provides 
experimental justification for the model used by Haasen and Leibfried 
(1951-52) to explain the orientation dependence of strain hardening 
curves. 

The interactions between dislocations on secondary systems with the 
primary system clearly provide a source of work hardening although, if 
the preliminary stress strain curves are typical of all orientations of alloy 
crystals, then these results provide experimental evidence for a mechanism 
of linear hardening only. 
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ERRATA 


Photometric Measurement of Lonization in Cloud Chambers, by 1. BurtEr- 
WORTH, 1955, Phil. Mag., 46, 884. 


In the table of mass measurements, p. 892, the figures in the column 
Mass from visual estimation’? are upper limits, no lower limit being 
assignable. The limits given do not apply. The same applies to the 
value 1220 in the photometric measurements. 

In the abscissae of fig. 2 the small divisions represent 200 electron masses. 
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CXXVIT. REVIEWS OF BOOKS 


Annual Review of Nuclear Science, Volume 4, 1954. Editor James G. BECKERLEY. 
Published by Annual Reviews Inc., Stanford, California in co-operation with 
the National Research Council of the National Academy of Sciences. [Pp.483.] 
Price (outside the U.S.A.) $7.50 postpaid. 

Taz contents of the fourth volume of the Annual Review of Nuclear Science 

cover an exceptionally wide range as will be seen from the following list : 
Recent developments in proton synchrotrons by J. P. Blewett. 

Theories of photonuclear reactions by J. 8S. Levinger. 

Radiofrequency and microwave spectroscopy of nuclei by G. E. Pake. 
Standardization of radioactive sources by G. G. Manov. 

Fission radiochemistry (low-energy fission) by L. E. Glendenin and B. P. 

Steinberg. 

Stable isotope dilution as an analytical tool by M. G. Inghram. 

Nuclear particle detection (fast electronics) by R. E. Bell. 

Nuclear particle detection (characteristics of scintillators) by R .K. Swank. 
Nuclear particle detection (characteristics of Cerenkov counters) by J. Marshall. 
Alpha radioactivity by I. Perlman and F. Asaro. 

Posttronium by S. DeBenedetti and H. ©. Corben. 

The interactions between pi-mesons and nucleons by M. Gellmann and 

K. M. Watson. 

Heavy mesons by C. Dilworth, G. P. 8S. Occhialini and L. Scarsi. 
Penetration of heavy charged particles in matter by E. A. Uehling. 
Biochemical effects of radiation by K. P. DuBois and D. F. Petersen. 
Vertebrate radiobioloyg (lethal actions and associated effects) by J. F. Thomson. 
Radioactivity in geology and cosmology by T. P. Kohman and N. Saito. 

Ne ESM. 


Radioisotope Conference 1954. Sponsored by the Atomic Energy Research 
Establishment, Harwell. Proceedings of the Second Conference, Oxford, 
19-23 July. Edited by J. E. Jounston, with the assistance of R. A. Fatres 
and R. J. Minuerr. Vol. 2, Physical Sciences and Industrial Applications. 
[Pp. ix+223.] (London: Butterworths Scientific Publications, 1954.) 
Price 45s. 

THE second volume of papers presented at the Radioisotope Conference at 
Oxford in 1954 presents a fascinating account of the use of radioisotopes in 
Physical Sciences and Industrial Applications. The twenty-four papers 
presented cover four main fields: Chemistry, especially exchange reactions 
and polymerization ; Metallurgy, including segregation, purification and self- 
diffusion in metals and alloys; Physics, including some new technical devices 
and Industrial Applications. To physicists these papers appear to form a 
more united group than the medical, agricultural and biological papers con- 
tained in volume 1. Nevertheless, it is already clear that radioisotope methods, 
as a tool of scientific research has reached roughly the same position as electron 
microscopy. The basic techniques and underlying principles are well estab- 
lished. Modifications and improvements of general importance may from 
time to time be reported. But on the whole individual applications of the 
technique to a particular field form an integral part of the field itself rather 
than another facet of radioisotopy. With the growing application of this 
technique to numerous branches of scientific research it will become increasingly 
difficult to convene further Radioisotope Conferences of this type. 

The volume is well produced, the papers are brief and clear and the 
discussions lively—where necessary persistent—and stimulating. The reviewer 
found particular pleasure in reading Dr. D. 8. Taylor’s account of a novel use 
of radioactive phosphorus to study the movement of wool fibres during the 


process of drafting. D. Tasor. 
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Linearized Theory of Steady High-Speed Flow. By G.N. Warv. (Cambridge : 
University Press.) [Pp. xv-+243.] Price 30s. 


For a decade, the linearized theory of compressible fluid motion has been 
intensively pursued by a large body of aerodynamicists throughout the world. 
The pioneering exploration seems now to be mainly completed, and it is proper 
that summaries are beginning to appear in book form. High-speed aerodynamic 
theory in general has been the subject of several recent treatises of an encyclo- 
paedic nature. Professor Ward has chosen instead to produce a single 
volume of moderate size by confining attention to steady motions and to the 
linearized theory only. 

The linearized equation of fluid motion is simply Laplace’s equation for 
subsonic flows and the wave equation for supersonic flows ; all the complications 
arise in the bondary conditions. It is therefore appropriate that the first part 
of the monograph, devoted to general theory, concludes with a careful discussion 
of boundary conditions and the various subtle points connected with singularities 
at wing and wake edges, shock waves, and conditions far from the body. The 
second part treats special methods for thin wings, conical fields, and nearly 
cylindrical shapes. The last part is devoted to slender-body theory, a useful 
further approximation beyond linearization that is valid only near the surfaces 
of smooth slender shapes. A selected bibliography lists 200 original papers. 

The monograph is based upon an Adams Prize essay, and the author conse- 
quently admits to a strong personal bias in the selection of material. However, 
conciseness demands some definite viewpoint, and the author’s own well-known 
researches do serve as a suitable basis for a connected account. The mathe- 
matical apparatus includes vector notation and Hadamard’s finite part of a 
divergent integral, but the physically more intuitive concept of sources, sinks, 
and vortices is also discussed. Mo DoV,.D; 


BOOK NOTICES 


Harmonic Analysis and the Theory of Probability. By S. Bocuner. (Berkeley 
and Los Angeles: University of California Press.) [Pp. viii+-176.] Price 
358, 


Degradation of Vinyl Polymers. By H. H. G. Jenninex. Volume ILL of a 
series of monographs on physical chemistry edited by Erie Hutchinson of 
Stanford University. (London: Academic Books Ltd.) [Pp. 329.] Price 
$8.50. 


Applied Mass Spectrometry. The report of a conference organized by the 
Mass Spectrometry Panel of the Institute of Petroleum and held in London 
from 29th—-3lst October, 1953. Published by the Institute of Petroleum. 
[Pp. vii+-333.] Price not stated. 


Science and its Background. By H. D. Anruony. 2nd edition. (London : 
Macmillan and Co, Ltd.) [Pp. ix +337.] Price 20s. 


Bessel Functions for Engineers. By N.W.McLacutan. [Pp.xii |+-239.] (Oxford: 
Geoffrey Cumberlege, Clarendon Press.) 2nd Edition. Price 35s. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Same as fig. 5, after intense electron bombardment. 
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(Figures 2—6 on formvar pellicle base.) 


Electron micrograph of 100 4 gold film, between two 100 A bismuth oxide 
films after annealing at 450°c (x10 000). (On carbon pellicle base.) 
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Preferential alignment of ‘ microetchpits’ on the slip lines (Al etched 
immediately after deformation). 


_- 


Fig. 3 x 1200 | 


Microetchpits aligned on a single slip system. (Notice the difference between 
etchpits with a geometrical outline and the smaller and rounder micro- 
etchpits.) 


G. WYON and J. M. MARCHIN Phil. Mag. Ser. 7, Vol. 46, Pl. 24. 


x 450 


Arrival of slip lines at a grain boundary A, B. One of the slip. lines shows no 
etchpits except at the point of intersection. 
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‘Fig. % 600 


Electrolytically polished specimen, weakly stretched, polished again and etched. 
A new deformation after etching shows the direction of the slip lines. 


*600 


AI-Mn_ 0-56%, alloy homogenized, weakly stretched, polished and aged. 
Microetchpits are aligned along the direction of the previously active 
slip lines. 
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Conjugate slip continuous 


and two traces on the cross. slip 


crystal. 


As fig. 3 but towards the edge of the 
system. 
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